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Vibrations of Rotating Beams with Tip Mass 
By Wit1aM E. Boyce and Grorce H. Hanpetman, Troy, N. Y., USA?) 


: 1. Introduction 


_ The problem of determining the influence of a tip mass placed at the free 
end of a uniform bar rotating at a constant speed about an axis at the other 
_ end and vibrating transversely to the plane of rotation (see Figure 1) has been 
_ considered previously [1]?). It was shown, through the use of Rayleigh-Ritz 
_and Southwell methods, that the presence of the tip mass depressed the lowest 
_ frequency. For this case, then, the effect of the tip mass on the kinetic energy 
_ outweighs the contribution to the potential energy produced by centrifugal 
force. Dr. KURT HOHENEMSER, of McDonnell Aircraft Co., pointed out in an 
informal oral discussion that this result agreed with calculations made by his 
group for the lowest frequency. It appeared to them, however, that the second 
and higher frequencies behaved in the opposite fashion with an increase in 
tip mass raising these frequencies. Furthermore, he suggested that frequencies 
higher than the first are of more direct concern to the designer of helicopter 
blades. 


Figure 1 
Rotating beam with tip mass. 


The purpose of this investigation is an attempt to analyze the effect ofa 
_ tip mass on frequencies higher than the first. It will be shown that asymptotic 
or singular perturbation methods can be employed to give information even in 
the case of variable flexural rigidity and variable mass distribution for high 
speeds of rotation. The uniform beam can be studied more explicitly through 
"these methods as well as by ordinary perturbation techniques, upper and lower 
bound estimates based on the Rayleigh-Ritz approach, and an extended 


1) Department of Mathematics, Rensselaer Polytechnic Institute. 
2) Numbers in brackets refer to references page 391. 
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370 Wiiu1am E. Boyce and Grorcr H. HANDELMAN ZAM 
Southwell procedure [2]. The general beam will be considered first and the 
detailed analysis for the beam of uniform section will be discussed later. It 
will be seen that these methods allow one partially to verify Dr. HOHEN- 
EMSER’S remarks. 

The asymptotic or singular perturbation method is partly an extension and 
partly a special case of a general theory developed earlier by Moser [3]. An 
attempt to provide a broader extension of MosEr’s work on general linear 
differential equations is now being carried out by Professor JosEPH B. KELLER 
of New York University and one of the authors of the present paper. The 
chief difference between the problem considered here and the general one of 
[3] lies in the presence of the eigenvalue and expansion parameter in one of 
the boundary conditions. This boundary condition is due to the tip mass. 


2. The equations of motion 


The deflections of the beam are assumed to be sufficiently small that the 
ordinary, linear, Euler-Bernoulli theory can be applied; in addition, shear 
deformation and rotary inertia are neglected. Let us consider a coordinate 
system which rotates with the beam at the constant angular velocity Q (see 
Figure 1). The origin 0is at the clamped end of the beam, the x-axis sweeps out 
the plane of rotation, and (x, t) measures the deflection out of the plane of 
rotation. Let EJ represent the flexural rigidity of the beam, o the mass per 
unit volume, A the cross-sectional area, and M the mass of the bob. All but 
the last can be functions of position x in 0 < x </, where / is the length of 
the beam. 


q+Aq 


é . m-+dm iy 


— 
p 


fc ane 


Figure 2 


Stress resultants on a typical beam section. 


The equations of motion can be obtained by considering an elementary 
section as shown in Figure 2 running from x to x + Ax. Let p(x, #) represent 
the tensile force at x, g(x, t) the shear force, and m(x, #) the bending moment. 
The equation of motion in the y-direction thus becomes 


2 
br 8 one (2.1) 


3 


Vol. XII, 1961 Vibrations of Rotating Beams with Tip Mass 371 
and the moment equation about one end of the element is 

om oy 
_ Differentiation of Equation (2.2) and substitution of Equation (2.1) yields 


02m 0 oy\ 0? 
Ox® Ox (65>) —@A On? * 


a 
| From the Euler-Bernoulli law, m = EI 0?y/0x?; the last pause becomes 


ga (ETS) - £ (0%) --04 2. ae) 


- The total tensile force at the section x is 


1 


p(x) = [| 08) A®) QEdE+ MIE. (2.4) 


x 


4 The beam is clamped at x = 0 and hence satisfies the boundary conditions 


At the opposite end, the bending moment must vanish, or 


Ory 
j BE op 0s 
a 

% 

| 


The final boundary condition is found from the equation of motion of the 
mass M in the vertical direction; namely 


—g=M s at x=7. 


However, g may be eliminated in favor of m by means of Equation (2.2) with 
_ the result 


The Euler-Bernoulli law and the fact that (J, t) = M 1 Q? reduce this equation 
to the boundary condition 


o (EIS oy) — M 1p SL = Mo ag oie 


Ox 
These equations can be simplified considerably by introducing harmonic 
oscillations (x, t) = w(x) exp(i4é) and the new non-dimensional variables 
~s=x/l and v(s) = w(x)/2. In addition, non-dimensional coefficients and para- 
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meters can be defined by 


1 
7 1 o(x) A(x) dx 
a= nO | 
0 


1 


asa f ole) A(x) dx 
Pst Bo (rere 
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where the maximum is taken over the range 0 < x </. It should be noted : 
that 0 <e(s) <1 on OXs <1. The quantity f? is the non-dimensional © 
eigenvalue which is being sought, «? is the non-dimensional speed of rotation, 


and y? is the mass-ratio whose influence is to be examined. 
The final eigenvalue problem can now be stated as 


s 


(ev")" — € \/ v(0) oda + y? 


v(0) = 0'(0) = 0, 
o"(1)=0, 
{e(s) v'} — yato’ = —72 fv, at s=1, 


) <preye, Oats ale 


(2.5) 


(2.6) 
(2.7) 


(2.8) 


where primes denote differentiation with respect to s. The problem is thus to 
find the eigenvalues f? and, in particular, to examine their dependence on the 
parameters « and y?. A closed form solution is, in general, almost impossible 
to find. For the. general case, then, an asymptotic or singular perturbation 
method, appropriate for large values of «, that is, high speeds of rotation, 


will be used. 
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3. Application of a singular perturbation method 


It is clear that a straight-forward application of a perturbation procedure 


_ in terms of 1/«? is not possible, for the lowest order term arising from Equation 


(2.5) must then satisfy a differential equation of the form 


Oe (« ole) ae 72 


) =prreye, 


which is of second order in the dependent variable v. On the other hand, 
since all four boundary conditions remain, the problem is over prescribed, and 
has no solution. Similar questions have been discussed rather generally by 
MosErR [3] using an asymptotic solution. His analysis does not include the 
present case in which the eigenvalue and expansion parameter appear in the 
boundary conditions. Nevertheless, his method can be extended to the case in 
question at the expense of a little more algebraic manipulation. 
Let us first set ¢ = 1/a? and uw = f?/a?; Equation (2.5) becomes 


é(ev")” — (¥ i/ r(0) odo + 7) = ris) 0. 


_ The term reduced equation will refer to the second order equation which results 
_ when ¢ is set equal to zero. The appropriate parameter for expansion is not ¢ 


itself, but one which depends on the difference in the orders of the original 
differential equation and the reduced equation; namely, 4 = e"?. 
The basic differential equation becomes 


Eee)" = (v i/ v(0) odo + 1) =pr(s)v; (3.1) 


s 


boundary conditions (2.6) and (2.7) remain unchanged ; whereas, Equation (2.8) 
is transformed into 


nf (ev) —yv' =—y uv at s=1. (3.2) 


One seeks a solution of the form 
v(s) = B (s; n) exp{y7* h(s)} | 
= »/ B,(s) nr exp {ni h(s)},  By(s) + 0. 
n=0 
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Substitution of this function into Equation (3.1) and division by exp {771 h(s)} 
yields the following equation arranged in terms of ascending powers of 7: 


fe|ew eo (ve) 2ae+ 7) B| 


s 


-+- 


L[oewar p+ tena s2e09Bersh'B 
n 


2 (/ (a) odo] (h” B+ 2h! m| 


(3.4) 
an [een ar a sem ns zene B'+6eh2B"+6eW h"B 


1 
roewe reapers e—( Je) adar) Bl 


s 


+nfeh"" B+4eh" B'+6ch" B"+4eh' BY +2e'h" B 
+6e'h" B’+6e h' BY +2e"h' B’+e"h" B 
+ 7 le BY’ +26 BY +e" BY )—purB=0: 


Since the coefficient of each power of 7 must vanish independently, and 
B,(s) + 0, it follows from the first term that 


1 
cH J ro) a4a-+8) =0. 


Ss 


Thus there are four solutions for h’; namely, 


E (/ \)" 
Fe bee v(0) odo + y? ; 


which can be integrated directly when the mass distribution is known. It is | 
important to observe that i(s) < 0 and h(s) 20 for) Ss = 1. | 
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zi Corresponding to each value of h is a solution v(s) which will be written in 
_ the form 

g Uy = My9 + 9 My + 7? My + + 7H Ry, (~My) , 

e Ug = Ugg + 1 May + 4? Ug + ++ + 7° Ro, (~ hy) , (3.5) 
eV exp (ths) = exp (7% hy) (ot % +o +2 S.) (mwa) | 
BW exp (— y-1hs) =exp (— 471hq) (wy + nwt ty’ T,) (~h). 


o- 
Moser [3] has shown that U,, U,,V, W will form a fundamental system of 
solutions of (3.1). Furthermore, the remainder terms R,,, R.,, S,, and T, are 
"functions of s, 7, and uw which can be estimated by 

j a! 
; | “Os! + 7 ioe 
-for0 <s <1and/=0,..., 3. This holds if 7 is sufficiently small and y is in 
, the interval | «| < u* where u* > 1 is any fixed number. The constant C, is 
independent of s,7, and w and the fundamental system is infinitely differen- 
_tiable with respect to s in 0 <s <1 and analytic in uw for | w| < p*. 

_ The equations satisfied by w;,;,v,;,#; can be found by substituting into 
Equation (3.4) and equating coefficients of like powers of 7. It is readily seen 


Ri, 


I 


0! 
+ | $57 Rao! 


, 
that for k= 1,2, m=0,1, and ” =2 

: t : 

- (/ | 

y ial [re o)odo+y? Un m = TU n» 

. | 

: (7 

| [ 10) o4a+08) A 
: s 

; 

x 


Equations satisfied by v, and w, can be found in similar fashion. 
The general solution is 


v= AU, + BU, + CV exp (y-1 ha) + DW exp (— 4-1), (3.6) 


‘where A, B, C, D are constants to be determined from the boundary conditions. 
Inserting Equation (3.6) into the boundary conditions (2.6), (2.7), and (3.2) 
‘leads to a set of four linear, homogeneous, algebraic equations in the unknown 
coefficients A, B, C, and D. The subsequent fourth order determinant guaran- 
teeing non-trivial solutions yields the secular equation from which the eigen- 
values can be found. 

| In order to simplify the structure of this determinant, for the present, at 
least, let us set 
| [VY exp (7-1 hy)] = V, exp (72 hs) , 


[W exp (— -1hg)]™ = W, exp (—n-th;), n=1,2,3 


me ali 
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where 


/ 1 i] , / if 
Vian tmV+V', Vea Phev +— ig V+ 2a V) +0", 


Va = n-2h2V + y-? (3 hy hg V4 3h V") 


+92 (WV +3MV' 43K) +", 
Waa oh Wi We, (3.7) 


Wo = 1-22 W — n-1 (a WH 2K W) + W", 
W, = — 23 W + y-2 (BM W + 3h2 W’) 
— nt (3KW" 43h. W' +h” W) +". 


The determinantal equation thus becomes 


U,(0) U;(0) 
U;(0) U,(0) 
oe V(0) exp {y~7 h(0)} V,(0) exp {7-7 h(0)} 
W(0) exp {— 77" h(0)}  W,(0) exp {— 971 hs(0)} 
m [e(1) Uy (1) — y? U(1) + y? w Uy (1) U7 (1) 
WP [e(1) UZ (1)]’ — y? U3 (1) + y? w U,(1) U3 (1) 4 
T,(1) V,(1) exp {77 hg(1)} 
T;(1) W,(1) exp {—-*hs(1)} 
where 


T(1) = [? {e(1) Va(1) + (1) Vo(1)} — y? V1) + y? w V(1)] exp {y-? hg(1)}, 
and 
T,(1) = [? {e(1) Wa(1) + e’(1) We(1)}— y? W,(1) + y? w W(1)] exp {— 9-1 h(1)}. 


The factor exp 4-1 {hs(0) — h(1)} can be taken outside this determinant 
by multiplying the third row by exp {— 9-1h,(0)} and the fourth row by 
exp {7-1 hg(1)}. The first two elements in the third row no longer contain an 
exponential factor; whereas the third and fourth elements contain the factor 
exp 7~1 {hg(1) — h3(0)}. Similarly, the first two elements of the fourth row 
contain this exponential factor; whereas, the last two elements have no eXxpo- 
nential factor. Since h,(s) is a decreasing function, the terms containing 
exp 9~1 {hg(1) — h5(0)} will vanish exponentially while the other terms will 
either grow or decay algebraically. Consequently, the terms with exponential 


p 
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_ The determinant J’ can be most conveniently expanded by the third row 
_ yielding 


“factors may be dropped and the determinant A replaced by A’ containing 
Bzero in each of these entries. 


70) exp 971 {Itg(0) — hg(1)} x 

q U;(0) yf? [e(1) UL (1)]'’-—y? U;(1) + y? uw U,(1) U, (1) 

|| U2l0) oP Let) UY — 9? Usd) + 7? w Ualt) Uy(1) 

S|} 0 WP fell) We(2) + el) MAY] — MA) + 2 WC) Wal) 

=[V’(0) + 4-2 45(0) V(0)] exp -* {hs(0) — h(1)} x 

5] Fu)? Lela) UL — 92 OY) + 9% aw Tal) Uy (1) 

; U,(0) m2 [e(1) USA)’ — 2 U4) + w UAL) UN) |=0. 

B} 9) Fle) WA) + 01) WO] — 7 MA) + 2m WL) Weld) 

Upon substituting Equations (3.5), and (3.7), one observes that the dominant 


~ term in A’ is of the order 7~*. Hence multiplication of A’ by 7° will leave each 
of the determinants with positive powers of 7 only. 

If the common factor exp 4-1 {h (0) — hs(1)} is dropped, the determinants 
may be expanded to give the following eigenvalue equation 


Riu) + Alo nts + hw rt ne Ron mw =9, (3.9) 


=r 


Tor, Se  e 


where one sees from MoseEr’s result on the behavior of the fundamental 
system of solutions that there exists a constant C, such that | R,(y, “) | <C, 
for all y> 0 sufficiently small and |u| < y*. Furthermore, each of the 


functions F,(u) ,... , #(u) is an analytic function of 7 in this region. 


ee a me Mia 


4. Evaluation of the zero and first order terms 


. Despite the generality of the coefficients, it is possible to make some 
' explicit statements concerning the eigenvalues, at least through zero and first 
- order terms in the asymptotic expansion. Since the functions f,.55. yt, are 
analytic in w and, as will be seen, Ff, has only simple zeros, Equation (3.8) can 


be inverted in the form 
B= pot nt + Me + te» 
_ where, again, the remainder term is bounded. It will be shown that sufficient 


information can be obtained about jug and jz, to shed some light on Dr. HOHEN- 
EMSER’S statement concerning the frequency behavior. 
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The differential equations satisfied by the functions ™,, have already 
been stated. In addition, the initial values of the functions and their first 
derivatives can be prescribed provided only that these conditions yield inde- 
pendent functions. For this purpose, one may set 


U4,(0) = 0, U4 (0) =1; 4, (0) =1, U59(0) es 
UsO=0, 4, O)=15 “a (O=15 t,,(0) = 0. 


One then finds, upon expanding Equation (3.8), that 
Fo(us) = hg(0) Vo(0) 5°(1) Wo(1) {= 7? myo(1) + y? w Mro(1)} - 


Since h,(0) + 0, 43(1) + 0 and V,(0) and W,(1) can be chosen so as not to 
vanish, Fj(u») = 0 becomes 


jo(1) +f Mo(1) = 0. (4.1) 
In addition, 


“(fr oto A | = 17 Osea (4.2) 


and 
(0) = 0. (4.3) 


Equations (4.1), (4.2), and (4.3) form the eigenvalue problem which describes 
the motion of a heavy flexible string, carrying a mass at one end, rotating 
about the other end in a horizontal plane, and vibrating transversely to the 
plane of rotation. The differential equation is precisely the reduced equation, 
and the boundary conditions are those which would arise from retaining the 
first of Equations (2.6) and (3.2) with 7 = 0. Thus the singular perturbation 
method has selected from among the original boundary conditions those which 
are appropriate to the zero-order problem. 

The n eigenvalue of the problem defined by Equations (4.1), (4.2), and 
(4.3) will be denoted by «\”) . Inspection shows that Uy9 = S is one of the 
eigenfunctions of this problem corresponding to the eigenvalue u = 1. This is 
indeed the lowest eigenvalue since U9 = S vanishes only at the boundary s = 1 
(see [4]). Therefore 


My’ = 1 corresponding to uw) =s. 


Unfortunately, the zero order term for the lowest mode yields no information 
about the dependence on y since it is independent of this parameter. It is 
therefore necessary to examine ju!) 

The solutions for higher modes are unfortunately not as obvious as the 
first; however, it is possible to obtain the behavior of u\”) for large order n. 
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(pu’)’+veu=0, o<s<n 


n> 2 
ue) ~ x Pen (4.4) 


[i yee | fr dere] ds 


s 


‘Equation (4.4) shows that as y? increases, jj (*) will also increase. Therefore, for 
sufficiently high speeds of rotation and sufficiently high modes of vibration, 
‘the natural frequencies will be increased if the ratio of the tip mass to the 
“beam mass is increased. 

_ Let us now return to the computation of the first order term 4. Carrying 
‘through the necessary algebraic calculations in Equation (3.8), the function 
_F,(u) can be found. The work can be considerably simplified, however, by 
“observing that a solution of the form yu = [fo + [, is desired and the functions 
Ff, and F, are analytic in y. Thus a solution of 


Fy (Uo +: 77 Pa) + 1 (Uo + 1 fa) = 9 


is to be found. Applying a Taylor-series expansion and keeping only the first 
two terms yields 


‘Since Fo(1uo) = 0, 


ee she 

4 iy OFM) 
3 Ou 
- order to find j,, one need only evaluate 
; Alu Vo(0) 52(1) Wo(t) [— 7* Usy (1, Mo) + Y? Mo M20 (1, Lo) ] 
j is 1s(0) hg2(1) Vo(0) WolH) (4.5) 
| eee PUE0. (1, pho) + 9 tag (Ls to) + po “sa (1, yo) | 

The function E(t) has been found by expanding A’ and making use of 


Equations (4.1) and (4.2). 


ae 


| 
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It is not necessary to determine the functions u) and Ou9/0u explicitly 

‘since Equation (4.5) can be rewritten in terms of “4 ) alone. The derivatives 
with respect to the eigenvalue ~ can be found by considering a perturbation: 
procedure. If one writes | 

b= ft 84, 
then 

Uyo (S, ) = Myo (S, Mo) + Ou(s) , 
and 
_ Otto (S, Mo) 

u(s) = a : 
Applying the usual perturbation expansion to Equations (4.2) and (4.3), one 
obtains 


if t 
(/ rodo+ *) Ur9(S, ra + My 7 Uy9(S, My) = 95 Uy9(9, My) = 0, 


(/ rodo+ | u’(s) 


Multiplying the second of Equations (4.6) by uyo(s, uo), integrating from 0 
to 1, applying integration by parts, and the first of Equations (4.6), we find 
that 


(4.6) 


, 


t for US) + 7 Uyo(S, Mo) = 0, u(0)=0. 


1 
Uyo(1, Mo) y2 w'(L) — W1(1, Mg) v2 w(1) + / r 024(8, flo) ds = 0. 
0 


Now %40(S, (4) also satisfies the boundary condition 


uo(1, Lo) = Mo Uy9(1, Lo) 


Thus 
x 
Uyo(1, Mo) y? [ee (1) — fy u(1)] + / r Uio(s, [4) ds = 0 
and hence : 
OF . 2 
poe ah 7 19 Y Uio(S, Mo) aS 
Ou (Ho) h;(0) Vo(0) h?(1) Wo(1) | ernie a | a a M49(1, Lo) 


Similarly the terms in F,(u9) involving U(1, Mg) can also be written in terms 
Of U49(S, Mo), for to satisfies 


1 / 
i(/ vy o do =P | Us9(S, | Se Lo Uo9(S, Ho) =0. 


ae. 
oy 
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Pistiptying by %19(s, #9) and integrating from 0 to 1, one finds 


a U49(S, [Lo) ( / vrodo + *) Uso (S, | ds 


Ss 


1 
se Ho [7 U19(S, Lo) Mao(S, fo) 4S = O. 


0 


If this expression is integrated twice by parts and the boundary conditions 
(0, fo) = Ma9(0, Mo) = 1 applied, the desired expression for #9 is found to be 


1 
fr odo+ yy 
4" Uso(L, Mo) + Mo Y? Uao(1, Mo) = - 


4 Uyo(1, Mo) 
Now 
fh 


1/2 
h,(0) = — [e(0)]~1? \/ rodo+ r| : 


0 


1 1/2 
Eo {f+ odo + || 
My = 1 - “I ce 7) 
[1 whols, Ho) ds + 7? Who(L, Ho) 
0 


Thus the linear term can be expressed entirely in terms of the eigenfunction 
of the reduced problem which in itself was required for the zero order problem. 
It is interesting to note that the denominator in Equation (4.7) is essentially 
to total kinetic energy of the string and added mass which is the physical 
‘counter-part of the reduced problem. 

For the lowest eigenvalue 


1 1/2 
lao {/ yodo+ || 


0 


(4.8) 


b= 1 
frsdst? 
0 


1 1 
[reds= | rstds 
0 0 


it is easily seen that du,/dy? < 0 and hence py is a decreasing function of y?. 


Observing that 


ee ee ee eS Oe me eee 


—— so 
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From the two term asymptotic expansion for yu it follows that for suffi- 
ciently high speeds of rotation the lowest natural frequency decreases as the: 
tip mass is increased. Furthermore, the limits obtained in Equations (4.4) and! 
(4.7) when y? approaches zero are finite. Thus the problem for the beam without! 
tip mass can be studied as an interior limit problem. It cannot be handled! 
directly by the singular perturbation method since the underlying theorem: 
requires that the coefficient of the leading term of the reduced equation shall] 
not vanish in the closed intervalO <s <1. 

Finally, the Southwell lower bound technique [6] yields as a lower bound! 
for the lowest eigenvalue, 

ee ee 


where y is the lowest eigenvalue of the problem 
yey) =yry, a See te 
v0) =v (0)=0, vf) =07 . 
y? (ce 0")’ —y? v' = —y* vv at s=1. 


Thus a one term asymptotic expansion will be a lower bound for the eigenvalue.. 
On the other hand, the second term in the asymptotic expansion contains a. 
factor which is linear in 1/7 whereas the Southwell bound does not. Hence for- 
small 7, the asymptotic result will exceed this bound. 


5. The uniform beam. Singular perturbation methods 


More specific results can be obtained for the case of the uniform beam. 
since the reduced problem is then capable of explicit solution. In this case: 
e(s) = r(s) = 1 and the boundary value problem becomes : 


wee ‘a 1 d 
nv — |» tr —s) +74) ee ee | 
v(0) = v'(0) = v"(1) = 0, (5.1) 
yo” —yv' =—yuv at s=1, | . 


The reduced problem now takes the form 


— [fF 0-9 47h] =n om, Pie Gicke bp.. 
(5.2) 
U9(0) =0=-— us, ALS + U9(1) . 


In order to solve the differential equation in problem (5.2), let us rewrite. 
it in the form 


{(L— 3? + 2%) wo} + 2 myo = 0 
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4 and set x= s/V1 + 22. + 2y?. Under this transformation the differential equation 
_ becomes 


d 
2fa-» Te} + 2 ty = 0, O<x< ere (5. 2a) 
or LEGENDRE’S equation. Setting 
Rk+1)=2p (5. 2b) 


the solution satisfying 1)(0) = 0, du,)(0)/dx = 1 may be written in terms of 
the Legendre function P,(x) and the associated Legendre function Q,(x) as 


io = sin ** dorcel JG ear, 
MRS he io 3) Die” 
Se ka hea 


- r(>+1)r(4) ON eae ea 


The eigenvalues jy can now be found by applying the second boundary con- 
dition in Equations (5.2). It should be noted that even this is somewhat 
awkward since the Legendre functions are not of integral order unless y? = 0. 
There are, however, several techniques which lead to very close approximations. 

Upper bounds for the eigenvalues can be found from the Rayleigh-Ritz 
technique. The standard methods of proof [7] will show that the lowest eigen- 


value is given by 
. i 


[Q—s*+2y%) w’? ds 


2 4 = min 


1 
[w* ds + y® w?(1) 
0 


where the only requirements on w are piecewise continuity of the first deri- 
vatives and w(0) = 0. The second boundary condition in Equation (5.2) will 
appear as a natural boundary condition. In fact, the resulting natural boundary 


condition is actually 


ee eg el) — pel) = 0. (5.4) 


Se SE ee Ne NI ee Te ee ee nC ee ee geen Leen Cen rn? 
ACA * ies ts as vv fee * 


“Tf. 


The limit term is of no importance except in the case y? = 0. Furthermore, if 
one derives the basic equations of the original rotating beam problem from a 
minimum principle, a similar limiting condition, guaranteeing no logarithmic 
singularity at s = 1, appears in the natural boundary conditions. 


; 

7 

| 

384 Witi1am E. Boyce and Grorcre H. HANDELMANN zAMe 

Since the solution of the eigenvalue problem given by the differential 

equation and first boundary condition in (5.2) with y?=0 and the second 

boundary condition given by Equation (5.4) is an odd Legendre polynomial, 
a trial function for the minimum principle has been taken as 


w=a P,(s) + b P,(s) + ¢ P;(s) . 


Here P,, P;, P; are Legendre polynomials and the coefficients a, 5, c are 
to be found from the minimizing condition. The lowest root of the resulting 
cubic equation is the lowest eigenvalue. The other two roots are upper bounds 
for the second and third eigenvalues. The upper bound for the second root as 
a function of y? is labeled RAYLEIGH-RITz in Figure 3. 


10:0 [ 
9-0/ High order 
Perturbation in y? Poe 
sf an 
80 
fe Asymiotic. in oe 

70; 
y? 
ae 

60 iy 

1 eee 1 L } 
01 02 03 0-4 05 
jag 
Figure 3 


Approximative eigenvalues of the reduced equation as a function of tip mass. 


A second method of approximation makes use of a truncated power series. 
It is more convenient to transform the equation into the form (5.2a) with k 
defined by (5.2b). The solution which satisfies 7, (0) = 0 and U;,(0) = 1 can be 
written as the hypergeometric series : 


i ao ATM ym MO 2) 2) Sa 2 
m=1 


(2m +1)! 


Differentiating this series and substituting the results into the second boundary 
condition yields an equation involving k and even powers of the quantity 6 
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“defined by bee watiint 


(1 + 2 O?)102 3 


Terms of order 6® have been neglected and the resulting equation solved for 
ee olution. The final results are marked Truncation in Figure 3 where it should 
"be noted that the exact solution for y? = 0 is 6.0. 

A third method is appropriate for small values of the added mass ratio and 
_ consists of a perturbation scheme in y?. For this case, expansions of the form 


yy = Dw, (2 9)", w= Dy, (2 9)?" 
n=0 n=0 


are taken. In order to take the case of the lowest order term, corresponding 
to y = 0, it is necessary to consider the boundary condition at s = 1 in the 
form given by Equation (5.4). Save for this remark, the procedure is com- 
ad standard and will therefore not be discussed in detail. 

It is found that the coefficient of the y? term vanishes and the solution for 
Ehe second eigenvalue up to terms of order y* is 


w= [6 + 420 y* + O(y*)] . 


As indicated in Figure 3, the results are only good for small values of y?. 
The same procedure has been applied to the m’* eigenvalue where it has again 
‘been found that the y? term in missing, and the coefficient of y* is positive 
provided that = 2. Consequently, for small values of y? the n‘ eigenvalue 
Ho is an increasing function of y? for n = 2. Recalling that up is the first term 
of the asymptotic expansion, we find a somewhat stronger result than that 
found for the beam with variable section; namely, for high speeds of rotation 
and small ratios of tip mass to beam mass, the second and higher frequencies 
_ increase with the tip mass. 

Asymptotic solutions for large y* can also be obtained by setting 


1 be 
oe R= 
22’ ye 


Dl nat ad 


and seeking solution in the form 


= 2 oe 2 .2n 
tee) w,e", A =A, e?, 
n=0 n=0 


valid for small values of ¢. The resulting differential equations have constant 
coefficients and the first two terms in the perturbation series can be readily 


found. The n‘# eigenvalue (n = 2) is given by 


oy 


pip = 2 (n —1)2aty® + [F + 4 (n — 128] +0(—r). 


ZAMP XII/25 


ew ee ee eS ee eevee ye eT 
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It should be noted once again that for 7 > 2, the eigenvalue is an increasing 
function of y?. Thus our conclusion previously found to hold for small ratios 
of tip mass to beam mass is now seen to be valid for very large ratios. The 
results for = 2 are also plotted in Figure 3 and are seen to be reasonably 
good even for moderate values of y?. 

Finally, for the sake of completeness, Equation (4.4) for high order modes 
has been evaluated for this special case for 7 = 2. One should not expect good 
agreement for the approximation in this case, but Figure 3 indicates that the 
results are not completely unreasonable. 

Once jy has been determined, the function 9 is fixed by Equation (5.3), 
and the second coefficient 4, in the asymptotic expansion can be found by 
evaluating Equation (4.7). This calculation leads to the results shown in Table 1 | 
on this page. 

Table 1 


First mode | Second mode 
y | oe Basy | Buy | y | o Basy Buy 
0 Z 2:87 4-10 0 5 14-34 25:34 
-f 4:95 5.45 10 26:67 33-32 
6 6-98 7:16 15 38-95 43-34 
ii 7:99 8-10 0:25 a 17-21 VAM fs 
0-25 2 2:64 3:25 10 31-33 33-59 
4 4-68 4-85 15 45-35 46-23 
6 6-70 6:70 0-50 5 20-40 23:49 
0-50 ; 2:53 2:93 10 36-50 36:74 
4 4:56 4-60 15 52-43 51-56 
6 6:57 6:46 1-0 5 26-44 25-46 
0:75 2 2:46 2°65 10 46-32 42:92 
“4 4-49 4-40 15 65:87 61-63 
6 6°50 6:30 


In this table B4sy is calculated from the two-term asymptotic series, and 
1 1/2 
Buy = 3 (8; ae pi.) 


where #7 and £7, are lower and upper bounds, respectively, for 6, found by the 
methods of Section 6. An inspection of the table reveals that the asymptotic 
results are valuable even for fairly small values of «. In general, the larger the 
tip mass, and the lower the mode of vibration, the better will be the asymptotic 
approximation. 
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6. The uniform beam. Upper and lower bounds 
on the second eigenvalue 


oa 
to 


a 


In the case of the uniform beam, a method previously used in [2] can also 

7 be applied to yield upper and lower bounds for the second eigenvalue for this 

case in which there is a tip mass. Since the technique has been explained in 
_ detail in the earlier reference, a brief outline only will be given here. 

The Rayleigh quotient for the eigenvalue problem set by Equations (2.5) 


4 (2.8) for the case of the uniform beam is 

e pe Dee 26 Dio, 2) 

.. RW) = Gy) — Hee, v) 

_ where 

a 1 1 

a 6 i] MW ‘il , 

y D(u, 0) = fo" ds+ Go | vw (L—s*+ 29) ds, 
3 1 

j H(u, v) = / vu ds + y® 0(1) u(1) . 

i 0 


4 


_ The following minimum principle can be proved for the lowest eigenvalue; 
' namely, 
4 

4 

~ where u has a piecewise continuous second derivative and satisfies boundary 
- conditions (2.6). The minimizing function is the first eigenfunction and satis- 


- fies Equations (2.7) and (2.8) as natural boundary conditions. Corresponding 
e principles can be proved for the n eigenvalue ; and, on the basis of these, the 
z 


P= min R(u) 


~ 


; Rayleigh-Ritz method can be applied. 
Thus, using the admissible function 


; 
: 
: 
, 


A u=as*+bs?+cs* 


- in the minimizing principle with a, b,c as free parameters leads to a cubic 
equation in the approximate frequencies. Each of these roots, in turn, is an 
_ upper bound to the three lowest eigenvalues. 

Lower bounds can be found by a generalization of SOUTHWELL’s method 
based on a minimum principle for the sum of the first ” eigenvalues [8] which 
states that if the eigenvalues are ordered according to increasing magnitude, 
ep? = 6? = --- By, then 


ae Ue een 


oe ae [D(d1, $1) <a Dbn $n) 
ar (es Ore 


ee ee TS NS 


388 WiutaM E. Boyce and Grorce H. HANDELMAN _ ZAMPT 


| 
where the 4, are admissible functions for the minimum principle also satis- | 
fying 


a | 
H($;, $j) = 64 = ‘ahh | 
lo, t+] 
Consequently, 
Daeis: = joeuhal >) Dadi: ¢;) + min D> Dali: W;) (6. 1) 
t=1 dy +++) On faq Vr ++ Un faq 
where 


1 1 
D,(u, u) = I “u"* ds, Dz(4, 4) = = a / wu’? (1 — s? + 2 y?) ds 
0 


0 


and the y; satisfy only the single boundary condition p;(0) = 0, but are other- 
wise admissible functions fulfilling H(y;, y;) = 0;;. 
It is easily shown that 


min SD, ($;, $3) = 2K 
where ju; are the eigenvalues of the system 
ie se ws Ones See 
w(0) = w’(0) = 0, 
w'(1) =O, wo" (1) + ty? ww) = 0. 
This problem corresponds to the non-rotating beam with tip mass; and, 


having constant coefficients, can be solved by direct methods. Clearly the 


eigenvalues w; depend on y but not on «. On the other hand, it can be shown 
that 


min SDs (y;, y;) = ye 
vi <4 f=1 


where y; are the eigenvalues of the reduced problem, Equation (5.2), written 
in the form 


[Udy aS 2 2) w' +2% w=0, 
w(0)=0, lim (1 — s?4+ 2?) w’ — 242 vw(l) =0. 


Equation (6.1) then yields the lower bounds 


Beat, B+ Ree + wet i+y 


+ 
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é Bo = wy + oe + + 2 — By 
where f} represents an upper bound to fy. 

d The results of evaluating these bounds for the second eigenvalue are shown 
in Figure 4. It should be noted that the bounds are rather close. Indeed, the 
Tatio of the difference between the upper and the lower bounds for f, to the 
“mean value is less than 6% for all values of « and y shown. The dependence of 
the eigenvalue on 7” is clearly exhibited in these bounds. For the non-rotating 
beam, «? = 0, the higher the tip mass the lower the frequency. As «? is increased, 
f however, the situation is reversed. This cross-over in behavior takes place first 
for the higher values of y?. 


cane Fe 


3 
7s 
J 
‘ 
c 
> ‘ = i L | L if | J 
a 0 50 100 150 200 
= a? —_— 
Figure 4 
“a Frequency versus speed of rotation for several values of typ mass. 


4 It should be noted that for y? > 0, there is a value of « beyond which the 
‘jower bound curve rises above the upper bound curve for ‘* = 0." Hence this 
calculation proves that for sufficiently high rotational speeds the second 
eigenvalue is an increasing function of the tip mass to beam mass ratio. 


y 7. Additional effects 


-- 


Certain additional effects can be traced without too much difficulty. I, 
for example, the tip mass is of finite size, its rotary inertia can be analyzed in 
“the following manner. Let x be the radius of gyration of the tip mass about 
its own center of mass. Then the boundary condition which leads to Equation 


m 


~ 
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(2.7) must be replaced by } 
| 
0? ea a2e POM. _ | 

EIS, +M25(3)=0, x=. 


If the same substitutions are used as before with the additional parameter 
k =x/l, Equation (2.7) should be replaced by . 


rf e(1) v"(1) — 72 Bw v'(1) = 0. (7.1) 


The problem which is then to be examined by an asymptotic expansion consists 
of the differential equation (3.1) and boundary conditions (2.6), (3.2), and. 
(7.1). The change in boundary conditions effects the analysis only in the last 
column of the eigenvalue determinant. The last column in the two determinants 
appearing in the expansion of A’, Equation (3.8) becomes 


PS oft wg (1) e(1) — 92 Be Sp ui, (1) , 
WP Dm" U2,(1) e(1) diem (1) 
e(1) Wi2(1) W(1) — e(L) {Ag (1) W(L) + 2 AL(1) W’(L)} + 9? e(1) WA) 


~~ By {- ; ni(1) W(l) + Sei 
The dominant term in the determinants is now of order 7~? rather than n~* as} 
in the previous case. Consequently multiplication of A’ by 7? will change the: 
determinant to one in which the entries contain non-negative powers of 7. The 
first two entries of the third column remain the same, whereas the last element 
becomes 


m e(1) hs*(1) W(L) — e(1) 9? {ag (1) W(1) + 2 hg(1) W'(1)} + 8 (1) W"(1) 
+ y? RP whg(1) W(1) — ny? Rw W'(1) . 


In the previous expansion of A’ through the first two terms, the only 
contribution from the third column came from the term /{2(1) W)(1) appearing 
as a factor. This is still true in the modified expansion except that this factor 


is replaced by y? k? whj(1) W(1). Thus the functions A(u) and F(u) in 
Equation (3.9) become | 


Fo(u) = hy (0) Vo(0) y? k? wo hg (1) Wo(1) {— y? wig(1, uw) + y? wu uyo(1, w)}, 
Fi (uw) = Vo(0) y? k? w h,(1) Wo(1) {— y? U59(1, Mo) + Y? fo Uao(1, Mo)} - 


Since 44 = 0 is a trivial solution only, as can be seen from the Rayleigh quotient 
for example, fy is unchanged. Furthermore the presence of u in the emee | 
cative factor in Aj(u) does not add an additional term to 0F(11») /Ou since 
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A 7? Uio(1, Mo) + Y? Mo Mr0(1, Mo) } = 0. Consequently, 


aa | aA 


Om 


mu 
tS 


will also be unchanged. 


— 


_ The foregoing analysis thus shows that the rotational inertia of the tip 
_mass does not affect the first two terms in the asymptotic expansion and con- 
sequently will not change the conclusions regarding general behavior drawn 
from these two terms. 


8. Conclusion 


Several methods for obtaining approximate solutions to the eigenvalue 
problem posed by the transverse vibrations of a rotating beam carrying a tip 
_mass have been considered. These include asymptotic representations in terms 
_ of the rotational speed and upper and lower bound methods based on minimum 
principles. In the first case, explicit formula are given for the first two terms 
and the lowest eigenvalue is determined to within quadratures which depend 
- only on the section properties. Higher frequencies have been explicitly calcu- 
lated for the uniform beam as well as upper and lower bounds on the second 
frequency. 

These results show the following behavior with regard to the question 
“raised by Dr. Kurt HOHENEMSER on the influence of tip mass. For high 
_ rotational speeds and general section properties, the lowest frequency decreases 
as the tip mass is increased ; whereas, for sufficiently high modes, the frequency 
"increases with the tip mass. In the case of the uniform beam, upper and lower 
_ bound analyses show that the latter statement is true even for the second mode. 
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Zusammenfassung ~ 


Es werden die zur Rotationsebene normalen Biegeschwingungen eines rotie-- 
renden Stabes mit Endmasse studiert; insbesondere wird im Anschluss an eine: 
friihere Arbeit, die sich auf die Grundschwingung beschrankte, der Einfluss der: 
Endmasse bei grésseren Winkelgeschwindigkeiten auf die héheren Eigenschwin-- 
gungen untersucht. Dabei werden St6rungsmethoden und asymptotische Nahe-- 
rungen in Verbindung mit den klassischen Verfahren von RAYLEIGH-RITz und! 
SOUTHWELL angewendet. Fiir Stabe veranderlichen Querschnittes wird nachge-- 
wiesen, dass fiir Eigenschwingungen geniigend hoher Ordnung und hinreichend! 
grosse Winkelgeschwindigkeiten die Eigenfrequenzen durch die Endmasse erhoht: 
werden. Fiir Stabe konstanten Querschnittes werden numerische Schranken an-- 
gegeben, welche dieses Resultat bereits fiir die erste Oberschwingung beweisen.. 
Dieses Verhalten steht im Gegensatz zu demjenigen der Grundschwingung, deren} 
Eigenfrequenz durch eine Endmasse stets erniedrigt wird. 


(Received: Dezember 9, 1960.) 


Meéthodes de calcul des milieux hétérogénes multiplicateurs de: 
neutrons. — Application a un systéme de 5 sources de neutrons; 
au plutonium-béryllium dans un modérateur de graphite’) 


Par PauLt Gavin, CHARLES MANDRIN et BERNARD VitToz, Lausanne?) 


1. Introduction 


L’assimilation d’un milieu hétérogéne A un milieu homogéne ne peut se: 
faire que sous les conditions suivantes: 
a) Le milieu hétérogéne est finement divisé en cellules identiques. 
b) Le libre parcours moyen des neutrons rapides et thermiques est au moins: 
du méme ordre de grandeur que les dimensions de la cellule. 


1) Subside du Fonds National Suisse de la Recherche Scientifique, N° A 183. 
*) Laboratoire de Génie Atomique de l’Ecole Polytechnique de l’ Université de Lausanne. 
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_ Examinons un milieu constitué d’un petit nombre de parcelles fissiles dis- 
_ posées dans un modérateur homogéne. Nous appelons «parcelle» un corps homo- 
| géne, fissile, ralentisseur, capturant et de volume relativement faible par 
_ rapport au modérateur. L’hypothése a) (division fine en cellules) n’est pas 
_ Satisfaite. On peut chercher a appliquer la théorie élementaire de la diffusion 
en considérant que chacune des parcelles et le modérateur externe forment 
_ autant de milieux homogénes. Ces milieux homogénes doivent étre couplés, et 
les conditions aux limites auxquelles on aboutit sont inextricables, sauf dans 
~ quelques cas trés simples. 

-_ FEINBERG [1]3) évite ces difficultés en négligeant le volume des parcelles 
 fissiles, et en les assimilant 4 des sources de neutrons ponctuelles, linéaires, ou 
planes. L’intensité de ces sources peut étre positive ou négative, le signe 
 dépend du caractére dominant (fissile ou capturant) de la parcelle. Le flux de 
" neutrons engendré par de telles sources est facile 4 calculer en un point quel- 


~ conque d’un modérateur infini; en utilisant le principe de superposition, on 
_ peut connaitre le flux de neutrons total en tout point de ce modérateur. 

Nous proposons [2] un autre procédé de simplification. Le milieu est 
constitué de n parcelles fissiles de volumes non nuls. Isolons une parcelle 
désignée par l’indice k dans le modérateur infini. La théorie élémentaire de la 
diffusion s’y applique sans peine si la parcelle a une forme simple. L’adjonction 


Ty ES 


_ des m — 1 autres parcelles fissiles modifient le niveau de flux et le taux des 
 fissions A l’intérieur de la parcelle k. Nous supposons que le systéme est en 
régime permanent. II est possible de représenter Vinfluence des n — 1 autres 
_parcelles sur & par des «sources auxiliaires de couplage» placées a proximité de 
k. En appliquant la théorie élémentaire de la diffusion 4 la parcelle k isolée, 
mais entourée de ses sources auxiliaires de couplage, on obtient une carte de 
flux approximative, qui est suffisamment exacte a lintérieur de la parcelle 
fissile, et A proximité immédiate de celle-ci. 

Les deux méthodes ont trouvé une application lors de I’étude d’un socle de 
graphite destiné a l’assemblage sous-critique de l’Ecole Polytechnique de 
VUniversité de Lausanne (EPUL). Ce socle contient 5 sources de neutrons 
Pu-Be, et il était essentiel de montrer qu’un tel ensemble est bien sous-critique. 

Les calculs numériques montrent que les deux méthodes sont en bon accord 
si l’écartement des sources est grand. Par contre, pour des positions rapprochées, 
la seconde méthode néglige certains effets d’ombres, et surestime de ce fait le 


taux de fissions. 


ie POTN NEP UNS 


er, Sra Te 


2. Méthode de Feinberg 
2.1 Méthode de couplage 


La méthode de FEINBERG [1] néglige le volume des parcelles fissiles. Ceci 
4 des points, des lignes ou des plans. Ces éléments 


Ts fF wee x. ¥ “™ * -_ 
Pe ey ee ee eee ee 


revient a les assimiler 


3) Les chiffres entre crochets renvoient 4 la Bibliographie, page 414, 


ete eee ee) ee ee 


| 
: 
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géométriques ont la propriété d’émettre des neutrons rapides, et d’absorber : 

des neutrons thermiques et de résonance. Le ralentissement des neutrons 

rapides dans la parcelle est négligé. : 
La parcelle fissile désignée par l’indice k émet: 

dx: neutrons rapides par seconde par suite de fissions ; 

Guz: neutrons rapides par seconde par suite de réactions indépendantes du 

flux de neutrons environnant (fissions spontanées, réactions («; 1), etc.) ; 
et absorbe: 
1,: neutrons thermiques et de résonance par seconde. 


Nous avons: 
Irk = Nk te» War=cte, = Vpk + Iak + (1) 


nombre de neutrons de fission émis par la parcelle k. 
“Ik nombre de neutrons absorbés par cette méme parcelle. 
y,: nombre de neutrons émis a la suite d’une fission. 
Ne = My (Upx!Zax)- 
Les parcelles d’indice 0,1,...k,... sont disposées dans un modérateur 
homogéne et infini. Le flux thermique en un point quelconque du modérateur 
est donné par: 


®,,(1) = 5 Get =i) — ill t —T 


); (2) 


F (| r — r, |) est le flux thermique produit en r par une source unité de neutrons 
rapides placée en r,: 


F(t + ty |) = fan (\e* = te) Ml = em) a, (3) 


Qu (|r — 1, |) est la densité de ralentissement thermique produit en r* par 
une source unité de neutrons rapides, placée en r,, et f (| r — r* |) est le flux 
thermique engendré en r par une source unité de neutrons thermiques, placée 
en r*. Notons en passant que la dépression due aux absorptions est représentée 
ici sous la forme d’une source négative de neutrons thermiques, et que ce 
formalisme suppose 

— que le modérateur n’est pas absorbant, 

— que les captures de résonance sont négligeables. 

Ces deux restrictions ont été levées par FEINBERG. Nous les laissons de cété, 
pour examiner comment ce modéle de calcul trés simple peut s’appliquer a un 
systeme réel, ot! le volume des parcelles n’est pas nul. 

Les couches périphériques de la parcelle fissile exercent un effet d’ombre 
sur les zones internes, ot les neutrons thermiques ont plus de peine a pénétrer. 
De ce fait, le taux d’absorption i, ne dépend pas seulement du flux thermique 
qui existe dans la région de la parcelle, mais aussi des dimensions et de la 


SEA Nae 


AN 


Vol. XII, 1961 Méthodes de calcul des milieux hétérogénes multiplicateurs de neutrons 395 


- nature de celle-ci: 


(Ry) 
Vk 


(4) 


$(R,): flux thermique dans la région de la parcelle. Ce flux est supposé égal 
a celui qui régne a la surface de la parcelle réelle de volume non nul, de 
rayon R,. 

V,: coefficient d’absorption, a déterminer pour un élément isolé dans un 
modérateur infini. Le calcul de y peut se faire, par exemple, par une 
théorie élémentaire de la diffusion, ou en se basant sur des résultats 
expérimentaux. 


= 


Cette méthode rend trés bien compte non seulement des effets d’ombre dans 
une méme parcelle, mais aussi des effets d’ombre réciproques entre parcelles. 
Toutefois, ceci suppose que celles-ci sont effectivement trés petites par rapport 
aux distances qui les séparent, et par rapport aux libres parcours moyens des 
neutrons. En effet, les masses capturantes ont été artificiellement concentrées 
en un point (ou une ligne, ou un plan) et la correction introduite par l’utili- 
sation de sources négatives ne suffit pas pour rétablir parfaitement 1’allure 


_ réelle de la courbe de dépression. 


D’autre part, il a été implicitement admis que le flux thermique 4 la surface 
d’une parcelle est uniforme. Cette hypothése ne se justifie que si elle est de 
dimensions relativement petites. 

On admet également que tout l’espace a les propriétés de ralentissement du 


_ modérateur proprement dit. En fait, les parcelles fissiles, qui modeérent les 


neutrons d’une maniére différente, occupent un volume non nul. Ici encore, il 


faut les supposer petites. 
Selon les relations (1), (2) et (4): 


St {n: Fir; ee 1 (%;%) 2 Oj xh = —S Fak U(";%) » (5) 
ou 
1 te ys 177 = Ay 


6;, est le symbole de KRONECKER: 


j,=t-s jak, et 4,=90 s1 G#hk. 
L’équation (5) est valable pour la parcelle n° j; nous pouvons écrire n — 1 
autres équations pour les m — 1 autres parcelles, ce qui forme un systéme 
linéaire de  équations aux  inconnues 7,, proportionnelles au flux qui régne 
4 la surface des parcelles correspondantes. En résolvant, on trouve: 

— Le taux d’absorption de tous les éléments fissiles. 

— Le taux d’émission de neutrons rapides. 


| 
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— En régime permanent, le facteur de criticalité, que l’on peut définir par: 


3 I rk 


— (6) 
S Trp + Caz 


‘a 


— En régime permanent le flux rapide, la densité de ralentissement thermique 
et le flux thermique en tout point. 


En posant que le régime permanent peut étre établi en l’absence des sources 
Jy, on obtient la condition de criticalité: 


| ie P15) al (tne) — Va One ==) (7) 


Dans ce cas, N, est négligeable devant N, et, en faisant tendre le rapport 
N,/N,; vers 0 dans l’équation (6), on voit que / tend vers 1. 

En régime surcritique, la théorie élémentaire de la diffusion, qui exprime 
un bilan stationnaire des neutrons, ne s’applique plus; en pratique, on obtient 
dans ce cas des flux négatifs, qui n’ont naturellement aucun sens. Par corollaire, 
le fait d’obtenir des flux positifs indique bien que le systéme est en régime 
permanent. 

Le systéme linéaire 4 m équations et inconnues se simplifie notablement si, 
parmi les m parcelles fissiles, # sont identiques, et sont disposées d’une maniére 
symétrique. En effet, il y a alors p équations identiques, et le systéme se réduit 
a lordren— p+ 1. 


2.2, Calcul du facteur y, 
Nous avons: 


O(R 
tg od oa) , (8) 
et 
t = | Eo. Po, av , (9) 
Ve 


ou l’indice 2 se rapporte au groupe de neutrons thermiques. 

y; ne dépend que des propriétés de la parcelle fissile ;- On peut donc déter- 
miner ce coefficient en appliquant la théorie A deux groupes a la parcelle k, 
isolée dans le modérateur. 


2.3 Calcul de la fonction f(r; ,) 


Dans l’équation (5), le terme — i, f(r;,) représente l’effet en r; d’une source 
négative de neutrons thermiques en r,, d’intensité totale — 1, n/sec. 

La fonction /(7;,) se calcule en appliquant la théorie de la diffusion A une 
source de neutrons thermiques isolée dans le modérateur. Pour fixer les idées, 


Z 
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_ supposons que les parcelles fissiles soient des sphéres trés absorbantes; dans ce 
_ cas, les absorptions ont lieu principalement sur la surface de la sphére, et il 
_ convient de représenter la source négative de neutrons thermiques sous la 
_ forme d’une coquille sphérique mince, de rayon R, et centrée en r,. Dans ce 
7 cas particulier, on trouve en accord avec GLASSTONE et EDLUND [3] (p. 129, 
» tableau 5.95) : 

tod wean fo Pl [PF 


29 29 


_ Dz, et Lz, sont le coefficient et la longueur de diffusion du modérateur externe ; 
_ pour simplifier le probléme, nous admettons que ces caractéristiques nucléaires 
_ sont également valables a l’intérieur de la sphére, ot se trouve le mélange 
- fissile. 

2.4. Calcul de la fonction F(r,,) par la théorie de lage (figure 7) 


Reprenons |’équation (3), en posant r, = 0. 


--—~ 


- Sp! 


FA A “ 
7 \, Sik Fareelle 
{ 1 
/ |= 
\ parcellek jal=? 
4 \ =0 / 
. Sey Weel =r [= Poot 
Figure 1 


Signification des symboles employés. 


On obtient alors la relation simplifiée: 


F(”) = faalr*) le — 9) ar (11) 


Considérons a nouveau le cas particulier d’une parcelle fissile ponctuelle. Selon 
GLASSTONE et EpLuND [3], p. 180 [6. 138. 2]. 


exp (— 7*?/4 t) 
Ql7*) = r ne 


t = Age des neutrons thermiques dans le modérateur externe. 


Il y a symétrie sphérique de centre &, donc: 


art = 47 dr*. 


D’autre part, en se référant a la relation (10) on a: 


fit =) = safe or [- 4 - er [- a], 


398 Paut GAVIN, CHARLES MANDRIN et BERNARD VITTOZ ZAMP | 


ae P exp (— 12/42) Lag f Spek aad OE ian Se" | ar®. 13 | 
al / (4 a x)8"" 27 Dog (ex? | Ly | exp| L2q | te " 


Cette derniére intégrale se résout a l’aide de l’identité: 


i= 
[vex — (ax? + bx +0) du = — 7 exp Hees 
a 


— bx 
x io eae a erf(y)| (14) 
avec y = Va(x + b/a) et erf(y) = fonction d’erreur. 
Aprés substitutions, on obtient: 


F(r) = exp (— t/L3,) 


8 a Dog 7 
+7 y es | e ") y Vr | y 
f f —— 2sh nts 
. lexp($) er F yr Log meee Log oe 2 yr Log : Log ( ) 


2.5. Utilisation directe des résultats de la théorie a deux groupes pour 
déterminer F(r) et f (r— r*) . 


Nous avons vu qu’il fallait de toute facon appliquer la théorie a deux | 
groupes a la parcelle k isolée pour déterminer le facteur y,. Ce calcul permet de 
trouver directement F(r) ou f(r). En posant r, = 0 dans |’équation (2), et en 
considérant une seule source, il vient: 


PD, (") = % Flr) —%& Kr), (16) 
dou: 
sty) — Perlr) + te f(7) 
Fi) ti (17) 
y\ — Gk P(r) — Pal) 
tae ee ae x, (18) 


Nx» % sont calculés a l’aide de la théorie 4 deux groupes, appliquée a une seule 
source, de méme que le flux thermique ®,,(7). Le coefficient g,, est donné. 
(Production imposée de neutrons, par réaction («; ) par exemple.) 

En utilisant la relation (17), on peut trouver F(r), A condition de calculer 
f(r) 4 Vaide de (10). De méme, |’équation (18) permet de déterminer f(r), en 
supposant que F(r) est connu par la relation (15). On préférera naturellement 
éviter les applications numériques laborieuses imposées par la relation (15), ce 
qui impliquera le calcul de f(r) par (10). 
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: 3. Méthode des sources auxiliaires de couplage 

: 3.1. Emploi d'une source auxiliaire de couplage 

__ La méthode proposée ici considére des parcelles fissiles de forme simple: 


br 


— sphéres, 
— cylindres de longueur infinie, 
— plaques d’étendue infinie. 


Leena ent 


Soit le cas simple de deux parcelles fissiles, isolées dans un modérateur 
“infini (modéle I). Pour fixer les idées, supposons qu'il s’agit de deux sphéres 


identiques (figure 2); les raisonnements seraient analogues pour des cylindres 


a, 


yO AtN e 


u des plaques disposées paraléllement. 


milieu modérateur intini 


ch) Sao 8) 


> he we ee 


mat. 
fissile 


ah ad tla 


Figure 2 
Modéle I. 


lh ee 


es Delile |. el 


Nous utilisons une méthode approximative de substitutions, qui permet 
_d’exprimer les conditions aux limites sous une forme simple. Pour cela, nous 
avons recours 4 un systéme fictif représenté par le modéle II, constitué par une 
‘seule sphere fissile, isolée dans le graphite et centrée en C, (figure 3). En plus 


‘ milieu modérateur infini 


\ 


Figure 3 
Modeéle II. 


de la matiére fissile, et afin de tenir compte de l’effet de la parcelle 1 dans le 
modéle réel, cette sphére comprend: 


/ 
~~ o  e 
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— Une source fictive de neutrons thermiques, infiniment mince et uniformément 
répartie sur la surface de la sphére. Les neutrons sont émis selon une distri-_ 
bution angulaire isotrope, a raison de S,, n/cm*sec. 

— Une source fictive de neutrons thermiques, uniformément répartie dans tout — 
le volume de la sphére 2, isotrope également, et émettant S, n/cm*sec. 


Ces deux sources forment ensemble la source auxiliaire de couplage (Source 
ADC). 

Les intensités S,, et S, devraient étre telles que le flux du modéle II soit 
égal en tout point au flux du modéle I. Cette exigence est impossible a satis- 
faire, parce que le modéle II présente une symétrie sphérique centrée sur Cy, _ 
ce qui n’est pas le cas du modele I. 

A lintérieur de la sphére 2, au voisinage immédiat de celle-ci, et dans les 
régions ou la sphére 1 n’exerce que peu d’influence, il est possible de faire © 
coincider approximativement les deux répartitions de flux. Appelons «zone 
d’influence» de la sphére 2, ce domaine ow la substitution par le modéle II est — 
correcte en premiére approximation. 

De méme, dans la zone d’influence de la sphére 1, une carte de flux appro- 
ximative s’obtient a l’aide du modeéle III, qui découle du modele II par per- 
mutation des centres C, et C, (figure 4). 


milieu modérateur infini 


Hi fy , 
frYrruyffffféh 


Figure 4 
Modeéle IIT. 


Le taux des fissions du modéle II doit étre égal 4 celui de l’une des sphéres 
du modeéle I, ceci résulte de l’équivalence des flux. 


3.2. Détermination des intensités S, et S,, 


Dans le modéle I l’influence de la sphére 1 sur la sphére 2 se manifeste de 
deux facons: 


~ En tant que sphére absorbante de neutrons thermiques, la sphére 1 diminue 
le niveau général du flux dans la sphére 2. Cet effet sera négligé; nous revien- 
drons sur ce point sous chiffre 4. 

- En tant qu’émettrice de neutrons, la sphére 1 provoque une augmentation 
du flux et du taux des fissions dans la sphére 2. 


x 
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a Dans le modéle II, en premiére approximation, la source superficielle 

dintensité S,, se substitue seule a la sphére 1. La source volumique d’intensité 

'S, correspond a une correction que nous calculerons plus loin. 

_ Appelons: 

. 32 le flux thermique supplémentaire induit par la sphére 1 dans la sphére 2 

a (modéle I). Ce flux se superpose au flux 3%, que l’on rencontrerait 

: dans la sphére 2 si elle était isolée dans le modérateur: B,, = B3° + B37; 

@}2* le flux thermique supplémentaire induit par la source pipeenecie 

; d’intentisé S,, dans la sphére 2 (modéle II). Ce flux se superpose a 

B GY): BF, = OE + Oe: 

iy, le flux thermique induit par la sphére 1 du modeéle III dans les environs 

3 dé.Cg; 

gx le flux thermique induit par la source superficielle du modéle IV dans les 
environs de Cy. 


| Le modéle IV (figure 5) s’obtient 4 partir du modéle II, en remplacant la 
sphere fissile par une sphére modératrice, de méme nature que le modérateur 
externe, et en posant S, = 0. 


v 


milieu modérateur infini 


YH yy 
Yy 


/, / / 
Wb iibie 


Figure 5 
Modéle IV. 


L’équivalence des modéles I et II s’exprime sous la forme: 


OPO, tou 0272 23, 
dans toute la zone d’influence de la sphére 2. Donc pour la détermination des 
taux de fissions, nous appliquerons les modéles II et III. 

La coquille-source d’intensité S,,, centrée sur C,, représente l’apport de 
‘neutrons en provenance de la sphére 1, étant entendu que l’intensité d’émission 
de cette derniére est renforcée par la présence de la sphére 2. Donc le taux des 
fissions dans la sphére 1 doit étre déterminé en considérant la sphére 1 isolée, 
“mais entourée d’une coquille-source d’intensité S,, (modéle III). Ainsi, l’apport 
‘de neutrons en provenance de la sphére 1 se mesure par le flux —,, calculé selon 

le modéle III. Ce méme apport se mesure aussi par le flux p¥ détermine par 


le modéle IV. 


_ZAMP XII/26 
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La condition d’équivalence de ces deux représentations s’écrit: 
px =, dans toute la sphere 2. (19) 


Les fonctions py, et y* étant différentes, il est impossible de satisfaire a 
légalité (19) en tout point de la sphere 2. Cette difficulté peut étre écartée en 
faisant la remarque suivante: par la forte absorption des parcelles fissiles, 
Vapport de neutrons fait sentir ses effets essentiellement a la surface de la 
sphére 2. C’est pourquoi il convient de remplacer (19) par 


px(R) =9,(C2). > , (19 bis) | 


Le premier membre de cette équation représente bien l’apport de neutrons 
de la sphére 1 a la surface de la sphére 2 (modéle IV); le second membre est 
une moyenne de cet apport a la sphére 2 (modéle III). 

Dans ce qui précéde, nous avons tenu compte en particulier des neutrons 
créés dans la sphére 1, qui parviennent a l'état rapide dans la sphére 2 ot ils 
se ralentissent. I] a été implicitement admis que les propriétés modératrices du 
mélange fissile sont celles du modérateur externe. Si tel n’est pas le cas, il en 
résulte un supplément (positif ou négatif) de neutrons thermiques crées a 
l'intérieur de la sphére 2. Nous représentons ce supplément par une source 
volumique uniforme de S, neutrons thermiques/cmsec, S, pouvant étre positif 
ou négatif. Conformément a la théorie de diffusion 4 deux groupes de neutrons, 


cette intensité s’écrit: : 

5, = Dy, Cie 2aige (20) 

Avec: 

@;7, flux rapide supplémentaire induit par la sphére 1 dans la sphére 2. Mais S, 
conduit a une correction relativement faible, ce qui autorise la détermina- 
tion de ce flux par le modéle III; 

2’, section efficace macroscopique de ralentissement de la sphere fissile 2; 

~}, section efficace macroscopique de ralentissement du modérateur externe. 


3.3. Application de la théorie a deux groupes aux modéles II, IIT et IV 


Pour le modéle III, les équations différentielles A l’intérieur de la sphére 1 
s'écrivent dans la théorie 4 deux groupes: 


Dey V2@,,— 24 ¢P,, + k 2,9, + S,=0, (21) 


Dg, V?®,, cz 22, Po, a 21.P, + S,=0. (22) 
Indice 1: neutrons rapides, 
Indice 2: neutrons thermiques. 
S,: intensité volumique (en neutrons/cm*sec) de la source de neutrons rapides 


qui est indépendante du flux environnant. Cette source est supposée 
homogéne et isotrope dans tout le volume de la sphere fissile. 
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_ Les autres notations sont celles qu’utilisent GLASSTONE et EDLUND [3]. 
L’indice c se rapporte a la sphere fissile 1. 
| Les solutions du systéme formé par (21) et (22) sont: 


= Ag, sinpr Bo, say Se + S$, 

Be ea at Seer Zi (23) 
sin sh vy Soe 

®, ,(7) ise A», paola e + Bo, = a (k ae 1) (24) 


Nous utilisons les notations suivantes: 
sinus hyperbolique 
cosinus hyperbolique 
tangente hyperbolique 
coth cotangente hyperbolique. 
- Les quantités u,v, S,,,S2, sont définies par GLASSTONE et EDLUND [3] 
(page 242), A,, et B,, sont des constantes d'intégration. En se référant encore 


» 


4 GLASSTONE et EDLUND, il est facile de trouver l’expression des flux dans le 


Dem D2) age El (25) 


Y 


D,, = BY (r) = Acg exp (= 1/Lou) (26) 


\ 


P, (R) = PD, ,(R) , (27) 
P, (R) a P, (R) , (28) 
aD, (r)]  _ dP, 4”) 
D ae +], — Dy,|—38 le (29) 
d®,,(7) a®,,(”7)] _<¢ 
—D,,| "|, Ie D, |" |, = Sn: (30) 


Z xX, Fie = Sa + Sy : 

e-° sa + Bac Sire | #agHie, = 24, (k = 1)" | 

. : Si +S; 

Ap, pert Beg VreiAig Seg Z1e ~ Ag, Zag = Sy, (k — 1)’ (31) 


BPs. Xie HBae eet Din Sig’ 


2c ees ly a 
Le 2¢ 


.. Dz, Xi, + Be. Dac Wie Aig Dy, S3¢ Zi6 — Azg Dog 234 = Sth: 
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sin R sh» R | 

avec yo Sa yee | 
le R , le R 


i = 
ie Zag = [Z| | 
Les indices se rapportent aux dérivées de ces expressions pour 7 = R. 

Appliquons ensuite la théorie de diffusion au modéle IV. Ce systéme ns 
comporte qu’une source de neutrons thermiques; il suffit de considérer I 
groupe thermique. Soit py(r) le flux thermique a l’extérieur de la coquille 
source, et px,(r) celui qui régne a l’intérieur de cette coquille. On obtien: 
facilement les relations: 


Les conditions aux limites pour ry = R s’écrivent: 


pr(R) = %,.(R), (34) 
es dye agF a 
Diep Dees )_ = Su- (Ge 


En combinant les équations (34) et (35), puis en procédant a quelques 
substitutions, on aboutit a la relation: 


Su=9H(R) z— Q. (36) 


avec 


= Dag Zag {gh + AN + Dyy Zar {8% (Rae) _ sh (Rag) 
Q 2g Pe a as Do, X¥, \ Ris. : , R2 a [2 


Mais en vertu de (26) 


Y,(Co) = DY (Cop) eyed bane 2 (— Ro/Ly,) ae, exp (— R,/Ly,4) 


0 Ry 


R, est la distance séparant les centres C, et Cy. En substituant p, a px dans la 
relation (36) on trouve: a 


= the 
Sin = S3¢ 41, Q Z,, Pie Aig OP, (37) 


- P, _ ©xp - oe , et B= exp = pace! : 
1g ‘to 2g to 


Rappelons d’autre part la relation (20): 
S, - OY" (C3) (21. “e X12) : 


ais, en vertu de (25): 


A —R 
ig(C2) = kz, exp(z_*) Aye Lag ba: 


ona Ai, Ly, Py (21. — »1,) : (38) 


Les relations (37) et (38) permettent d’éliminer S, et S,, dans le systéme (31), 
qui se transforme ainsi en un systéme linéaire de 4 équations a 4 inconnues; 
en résolvant ce systéme, on trouve les 4 coefficients d’intégration A,,, By,, Ai, 

A..: a V’aide de ces résultats, on peut déterminer la répartition des flux 


a 2g? 
Jans le modéle III, et par symétrie, celle du modeéle II. 


3.4. Evaluation du taux des fissions dans une sphere fissile du modele I 


En vertu des hypothéses faites, nous avons: 


Ny=h[EoPre AV. = | Dae Pre AV, (39) 
volume d’une sphére volume de la sphére I 
du modéle I. dans le modele III. 


N;: nombre de neutrons rapides émis par unité de temps par suite de fissions 


dans une sphere fissile. 


aa 
f= 
Pag 
£; 
rs 


La seconde intégrale est facile 4 calculer 4 l’aide des résultats dont nous 


isposons. 


THA 


a 3.5. Généralisation de la méthode. Application a des systemes plus compliqués 


QAP 


sé en détail la maniére de traiter le cas de deux sphéres 
elques modifica- 


Se 


_ Nous avons expo 
identiques, placées dans un modérateur infini. Il suffit de qu 
tions pour adapter ces raisonnements au cas de » sphéres quelconques, disposées 


dans un modérateur infini. 

at e . , . A , e 

~ Pour la sphére &, les intensités Si, et S* doivent étre déterminées en tenant 
y arrive en remplacant 


compte des influences des » — 1 autres sphéres. On 
Péquation (19bis) par: 
S Pje(Cox) = pe(R); j,4=1...0. (40) 


Se wee 


arr a 
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De méme |’équation (20) se généralise: ; 


‘ 
‘l 


bs a (Ze 7 a1) - oi" e (4 


| 
| 
| 
Les relations (21) a (31) restent valables, mais toutes les grandeurs sont - 
munir d’un indice 7, les propriétés nucléaires et géométriques pouvant varie 
d’une sphére a l’autre. Si les sphéres sont toutes différentes, et si on les dispos. 
sans aucune symétrie, le systéme linéaire (31) se transforme en un systém 
linéaire d’ordre 4 . L’élimination des S/ et des Si, se fait d’une maniére analogu 
a celle que nous avons exposée. . 

Si, au lieu d’un ensemble de sphéres, on considére des cylindres infinimen: 
longs ou des plaques infiniment étendues, tout ce qui précéde reste valable, 
l'exception des fonctions X,,, Y,,, 2,,Z,, et de leurs dérivées. Dans le premie 
cas, nous aurons des fonctions de BEssEL; dans le deuxiéme, des fonction 
trigonométriques simples. 

Enfin, il est facile de voir que si parmi les  parcelles fissiles, # sont iden 
tiques et disposées symétriquement, l’ordre du systéme linéaire (31) se réduit : 
4(n—p+1). 


4. Applications numériques 
Milieu constitué de 5 sources de plutonium-béryllium et de graphite 
Comparaison des simplifications utilisées et de leurs conséquences 


4.1. Description du systéme étudié, choix des constantes nucléaires, étude d’un 
seule source 


Le socle de graphite de l’assemblage sous-critique de l’EPUL contien) 
5 sources de neutrons au plutonium—béryllium; c’est un prisme d’axe vertica 
et de section octogonale. Les 5 sources sont logées dans des canaux horizontaux 
la position radiale de chacune d’elles est réglable A volonté (voir figure 6). 


fash 
2 ote meen 
aie 2 
NSAWZZZRSSSBSSSNVZZANSSY 
took 2 Sh e Hee] —j—-——}- 


SS 


Figure 6 


Disposition des sources Pu—Be dans le socle de graphité. Elévation et plan. 
g: ecartement des sources. En hachuré: coupe dans le graphite. En noir: coupe dans les sources 
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4 Les caractéristiques des sources de neutrons sont les suivantes: 

j— Intensité: 5 Curie «, 8- 10° neutrons rapides émis par seconde 4 la suite des 
| réactions («; 1). 

3 Composition : alliage métallique de 80 g de plutonium et de 60 g de béryllium. 
= Forme et dimensions externes: cylindres, ® 28,06 mm, longueur 51,30 mm. 

— Gainage: acier inoxydable et tantale minces. L’influence de ces enveloppes 
est négligée au point de vue neutronique. 


Pour le plutonium métallique massif, nous avons adopté les valeurs sui- 


0,4,= 1104 barns, of, = /8/barns, v= 2,88. 


Les constantes utilisées pour le béryllium métallique pur sont: 


t (Age thermique de Ferm) = Li, = 102,7 cm?, Lé= Oise 
: Zabs th = 9,03 10-* cm-?, 4,=1,46cm. 
Les notations sont celles de GLASSTONE et EDLUND [3]. Pour le graphite, les 
valeurs suivantes ont été choisies: 


t (age thermique de Fermi) = L}, = 365 cm?, x,é = 0,06lcem=*, 
. ie 51 cm, 2, th = 0,327 -10-* emer. 


Le plutonium étant fissile, il convient de s’assurer que le socle ne constitue 
‘pas en lui-méme un systéme critique. Il s’agit d’un calcul de sécurité, aussi 
les hypothéses simplificatrices seront choisies dans le sens d’une exagération 
du niveau du flux thermique, donc aussi du taux des fissions. Dans cet esprit, 
‘on peut supposer que le graphite occupe un volume infini, alors que ses dimen- 
‘sions sont de l’ordre du métre. 

_ Ilest difficile d’appliquer la théorie de la diffusion élémentaire 4 un cylindre 
fissile placé dans un milieu modérateur infini. C’est pourquoi il faut commencer 
par remplacer fictivement les sources cylindriques réelles par des sources sphé- 
‘riques équivalentes. Les conditions d’équivalence devront surestimer d’une 
“maniére certaine le taux des fissions. Nous étudions d’abord le comportement 
‘dune seule source dans le modérateur graphite infini. 

_ Envisageons deux cas extrémes: 

4 a) Le flux thermique est uniforme dans tout le volume d’une source. Le 
“taux des fissions est proportionnel au volume de celle-ci, et la sphére équivalente 
doit avoir le méme volume que le cylindre réel. 

__ b) Le flux thermique est grand prés de la surface de la source, et subit une 
forte dépression a l’intérieur. Le taux des fissions est pratiquement proportion- 
‘nel a la surface de la source, et la sphére équivalente doit avoir la méme surface 


“que le cylindre réel. 


wt 


Yee wee 


408 Paut GAvVIN, CHARLES MANDRIN et BERNARD VITTOZ ZA 


Dans le cas a), la concentration en plutonium de la sphére équivalente est Te 
méme que dans le cylindre réel. Dans le cas b), le volume de la sphére équi- 
valente est plus grand que celui du cylindre réel; nous maintenons la masse 
totale de plutonium 4a 80 g, si bien que la concentration de plutonium dans la 
sphére de surface équivalente est inférieure 4 la concentration réelle. On sim- 
plifie bien dans le sens de la sécurité, car on se trouve au-dessus du rapport 
optimum Noyaux de Pu/Noyaux de Be. 

Les constantes nucléaires utilisées sont indiquées au tableau 1. 


Tableau 1: Constantes nucléaires 


Sphére de volume Sphére de surface 
équivalent équivalente 
R = 41,9642 cm R = 2,1412 cnx 
V, = 31,745 cm? Ve, = 41,122 cm$ 
Li =a 
D,, = 0,830 
LZ, = 0,0694 cm? IZ, = 0,0898 cm? 
2i,, = 7,008 em! 2., = 5,41 cm—? 
Ry ==) = 23053 
€ =p=fel 
Bm? = 0,010242 cm? Bm? = 0,010237 cm? 
fi, = 365 cm* 
Di = L183 em 


L3, = 2601 cm? 
2a, = 0,327 + 10-8 


Les équations différentielles de la théorie A deux groupes se résolvent 


aisément pour une sphére de Pu-Be dans un modérateur graphite. Nous avons 
défini: 


I 


N; = intensité des neutrons émis par fissions = h | 2i3,P2,dV, nis. 


Ve 


I 


N, = intensité des neutrons émis par réaction (x, m) = 8-108 njs. 


N; 


= facteur de criticalité = — : 
/ Negro, 


A = facteur d’amplification due aux fissions — NT. 


a 


Les résultats numériques obtenus sont indiqués au tableau 2 pour les deux 
équivalences étudiées: sphére de volume ou de surface équivalente. 
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Tableau 2: Une seule source 


Sphére de volume Sphére de surface 


fe équivalent équivalente 
N; (108 n/s) 1,166 1,280 
., Mi 0,127 0,138 

a 1,146 1,160 


_ Pour les raisons de sécurité invoquées plus haut, nous ne considérerons par 
‘la suite que la sphére de surface équivalente en remplacement de la source 
cylindrique réelle, le taux de fissions dans ce remplacement étant supérieur au 
taux de fissions de la sphére de méme volume. 


; 4.2. Cing sources, méthode de Feinberg 


Les cylindres étant remplacés par des sphéres équivalentes, il est possible 
d’appliquer sans autre les méthodes décrites sous chiffre 2. 
Nous avons considéré le cas ot 4 sources périphériques sont disposées d’une 
facon symétrique par rapport a la source centrale, et nous avons fait varier le 
“rayon @ de cette couronne externe pour déterminer la position la plus dange- 
reuse (voir figure 6). ; 
La figure 7 représente la répartition du flux thermique obtenu en appliquant 
‘la théorie 4 deux groupes a une sphére Pu—Be isolée, de surface équivalente. 
Ce méme calcul fournit: 
i si @,, 3, dV, = 6,233, 10°n/s; B,,(R) = 6743 n[cm?s , 


Ve 


iy ae = 0,0108 cm?. 
i 


_ Nous avons résumé au tableau 3 les principaux résultats de la suite des 
-calculs. 

Les calculs numériques ont été exécutés conformément a la méthode exposée 
“sous chiffre 3. On constate que l’utilisation de la théorie de l’Age (éq. 15) pour 
‘le calcul de F (0) donne des valeurs sensiblement différentes des résultats 
 obtenus par l’équation (17) (théoriea deux groupes). Ces différences se retrouvent 
“Jorsqu’on calcule les facteurs de criticalité. Il convient de remarquer que la 
théorie de |’Age néglige la présence du béryllium, tandis que la théorie 4 deux 
groupes en tient compte; nous avons ainsi mis en relief le réle trés important 
“que joue le béryllium au point de vue de la modération des neutrons rapides. 
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19140 n/om* s 


5-70° 
Ont 
Be © 
{n/cm’ s] 2 
| 10° 
Pe oh tes 70 15 20 25 
a eleinl 
Figure 7 
Carte de flux thermique 
R= 2,141 1 source (80 g Pu/Be 8 - 10° n/s) K = 2,053 


Par la suite, nous ne considérons que le second cas, F(e) calculé par la 
théorie a deux groupes, a l’aide de (17). 

Ces résultats sont en effet plus réalistes, et surtout tendent a surestimer le 
niveau des flux et le taux des fissions. 


Tableau 3: Cing sources, méthode de FEINBERG 


Théorie de l’Age Théorie 4 2 groupes 


t(@) ee Pe Roe TTL 
Facteur de Facteur de 
@ d’apré F(@) se F(e) \ 
prés (10) ‘ eer | 
eae 10-2 em=—2 daprés (15) creas: d’aprés (17) pee 


[F(@) calculé [F(@) calculé 


10-4 cm-* selon (15)] 10-4 cm? selon (17)] | 
praises erp ge ee ee 
74 1,146 14,74 0,155 27,59 0,321 
14,2 0,499 13,98 0,219 21,38 0,327 i 
21,3 0,290 12,74 0,280 16,40 0,298 . 
28,4 0,189 11,45 0,277 13,16 0,266 | 
35,5 0,131 10,69 0,247 10,41 0,238 
42,6 0,0953 8,26 0,196 8,26 0,215 
49,7 0,0711 6,97 0,165 6,60 0,197 
oo 0,000 0,00 0,0588 0,00 0,138 ; 
. 
| 
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4.3. Méthode des sources auxtliaires de couplage 


__ Pour des raisons de sécurité déja invoquées nous adoptons ici les résultats 
obtenus pour des sphéres de surface équivalente. 


1 
! 


05-10° 
i 
. 
| 3 
| 
03 | 
2 
| 
| 
| 
j g2 I 1 
| 
if 
YD 
n/cm’s g1 ! | 
| | 
| 
ip 
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—>— [cps 


Figure 8 


_ Répartition des flux thermiques @,, en fonction de la distance au centre d’une sphére de plutonium- 
béryllium. R Rayon delasphére. 1 Sphére isolée dans un milieu infini en graphite. Théorie élé- 
mentaire de la diffusion. 2 Sphére excentrée d’un systéme de 5 sphéres disposées symétriquement 
dans un milieu infini en graphite: 9g = 7,1 cm (voir figure 6). Méthode des sources ADC. 3 Sphére 

centrale du méme systéme. 


La figure 8 représente la répartition du flux thermique dans la zone d’in- 


fluence d’une source, pour @ = 7,1 cm. 
Le tableau 4 présente quelques résultats de calcul pour la source centrale et 


_pour les 4 sources excentrées, pour @ = 7,1 cm également. 


Tableau 4: Cing sources, @ = 7,1 cm 


Source centrale Sources excentrées 


[ Gx, Fac WV, = 4,31 + 10° njs [Gaq Ea¢ WV, = 4,01 + 108 n/s 
Ve 


V, 
N, = 8,85- 10% n/s N, = 8,23 - 108 n/s 
N,+ Na = 16,85 108 n/s NetNg == 10)25 —L0Sns 
Pour l’ensemble des 5 sources: 
N5 = 81,76 - 108 n/s 
fe 
A, = 2,04 


ee eT. 
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tableau 5). 


a5 
04 
03 
pbc 
ee ee es eS asymptote 
ca 70 20 30 40 50 
oe, 
Figure 9 


Variation du facteur de criticalité f en fonction de I’écartement o@ des sources (pour la définition 


de 9, voir figure 6). 
1 Méthode de FEINBERG. 
2 Méthode des sources ADC. 
3 Méthode de FEINBERG simplifiée. 


Tableau 5: Cinq sources Pu—Be dans modérateur graphite 


E Facteur de Facteur de Facteur de 
cartement des criticalité f 
criticalité f, criticalité f, ‘ 
sources @ ‘ F méthode de 
ah méthode de méthode des FEINBERG 
FEINBERG sources auxiliaires simplifiée 
Qa 0,321 0,510 0,490 
14,2 0327 0,401 0,381 
21;3 0,298 0,330 0,324 
28,4 0,266 0,280 0,280 
35,5 0,238 0,245 0,247 
42,6 0,216 0,219 0,221 
49,7 0,197 0,200 0,199 


co 0,138 0,138 0,138 
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On constate que la méthode des sources auxiliaires et la méthode de 
FEINBERG conduisent a des résultats en bon accord tant que l’écartement des 
sources est supérieur 4 20 cm environ. Pour de plus petites valeurs, les hypo- 
-théses des deux méthodes se trouvent plus ou moins, en défaut, et on constate 
“que la méthode de FEINBERG conduit a des valeurs de 7 qui sont systématique- 
ment plus faibles que celles calculées par la méthode des sources auxiliaires. 
Ceci s’explique facilement si l’on se rappelle que la méthode de FEINBERG 
ne considére pas seulement l’effet de dépression créée par une parcelle k dans 
‘sa propre zone, mais également la dépression exercée par les n — 1 autres 
_parcelles dans la z6ne de k. La méthode des sources auxiliaires ne tient compte 
que du premier de ces deux effets; si l’écartement des sources est faible, le 
second prend de l’importance, et entraine une diminution du niveau des flux 
et du taux des fissions. 
Pour justifier cette remarque, nous avons appliqué 4 nouveau la méthode 
de FEINBERG, mais en négligeant cette fois le deuxiéme effet de dépression. Il 
suffit de modifier l’équation (5) comme suit: 


S Wee ae) Oren Wn Onkt = —S Iak Fr) (42) 


Cette méthode simplifiée conduit a des résultats qui sont en meilleur accord 
avec la méthode des sources auxiliaires, les écarts diminuent surtout pour de 
faibles distances entre sources. Un désaccord de quelques pour-cents subsiste 
_cependant. Nous l’attribuons 4 des effets de volume; la méthode des sources 
 auxiliaires tient compte du fait que les absorbants sont répartis dans un 
certain espace, tandis que la méthode de FEINBERG les concentre en un point; 
en ce sens, la’ premiére est plus exacte. 
| Le tableau 5 réunit les résultats de calculs du facteur de criticalité f en 
- appliquant les différentes méthodes exposées. 
D’autre part, la méthode des sources auxiliaires permet de mieux tenir 
compte des propriétés modeératrices du béryllium qui sont nettement meilleures 
que celles du graphite. Cet effet a pour conséquence d’augmenter le niveau 
’ du flux et le taux des fissions; il n’est sensible que si les sources sont trés 
_ proches les unes des autres, 4 cause du volume relativement faible du béryllium. 
Ace point de vue la méthode des sources auxiliaires est plus exacte que celle 
de FEINBERG. 
Toutefois, dans ce cas particulier, l’effet des dépressions secondaires semble 
- nettement plus important que ensemble des effets de volume et de ralentis- 
 sement. Par conséquent, la méthode de FEINBERG conduit certainement a un 
- résultat plus proche de la réalité que la méthode des sources auxiliaires. Par 
contre, il n’est pas certain qu'elle surestime systématiquement le taux des 
" fissions, a cause des effets de modération en particulier. Pour cette raison, la 
~ méthode des sources auxiliaires est 4 préférer pour un calcul de sécurité. 
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Ces conclusions ne sont valables que dans le cas particulier de 5 sources de_ 
Pu-Be dans le graphite. Toutefois, on peut dire, d’une maniére générale, que — 
les deux méthodes en présence sont complémentaires, et il est intéressant de les : 
appliquer simultanément pour ]’étude des milieux hétérogénes a petit nombre — 
de cellules. 
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Summary 
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In the field of nuclear engineering, heterogeneous systems are generally preferred — 
to homogeneous ones. Therefore, difficulties arise while computing critical sizes. 
It is possible to avoid them by considering the heterogeneous media as homogeneous 
ones, and by properly defining the characteristics of the latter. Such a method 
cannot be used in the case of a small number of fissionable pieces, placed in an 
infinite homogeneous moderator. We compare two methods solving this problem: 
the heterogeneous method of FEINBERG, and the heterogeneous method of auxiliary 


coupling sources. 


(Regu: le 24 février 1961.) 


Minimum-Volume Design of Discs’) 


By Tr-Cutanc Hv?) and Ricwarp THorPE SHIELD®), Providence, R. I., U.S.A. 


1. Introduction 


The problem of this paper is that of designing a disc which will carry given 
loads and which will involve minimum volume of material. It is assumed that 
the material of the disc is rigid, perfectly-plastic and the yield condition of 
TREscA is used. The case of rotationally symmetric loading has been treated 
previously [1, 2]*) and some problems have been solved for a disc composed 
of a brittle material with reinforcing rods or cords [31 


1) The results presented in this paper were obtained in the course of research sponsored by 
the Office of Ordnance Research, Department of the Army, under Contract No. DA-19-020-ORD-5124. 

*) Division of Applied Mathematics, Brown University. Now at I. B. M. Research Center, 
Yorktown Heights, N.Y. 

%) Division of Applied Mathematics, Brown University. 

4) Numbers in brackets refer to References, page 433. 
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A design is of minimum volume if it is at collapse under the given loads 
and if it admits a collapse mode for which the rate D of dissipation of energy 
per unit volume due to plastic action is constant throughout the disc [1, 2]. 
Although this provides a direct procedure for minimum-volume design, it is not 
known whether such a design exists in general. In Section 2, three theorems 
are given which supplement the direct procedure. For plane stress under 
Tresca’s yield condition, it is shown below that the constancy of the dissi- 
pation rate D leads to essentially only four types of solution, although a 
particular design may involve more than one of the four types. Several illustra- 
tive examples of minimum-volume designs are given. The designs may be 
compared with those obtained by MICHELL [4] for structures restricted to 
consist of bars in tension and struts in compression. 


2. Procedure for Minimum-Volume Design 


It is required to design a disc which is just at the point of collapse under 
given loads acting in the plane of the disc and which is of minimum volume. 
The disc is formed by distributing a variable thickness h of a given material 
at points of a plane middle-surface A7. At the boundary of the middle-surface 
A,, each component of the applied forces (acting in the plane of A7) is specified 
except where the corresponding component of displacement is prescribed to be 
zero. Loads may also be prescribed to act on specified lines in the plane Ay. 
The assumption of plane stress will be made, and the material will be assumed 
to yield under TREscaA’s maximum shearing stress criterion, Figure 1. 


Figure 1 
The Tresca yield condition for plane stress. 


It has been shown [1, 2] that if a disc can be found which is at collapse 
under the given loads and which has a collapse mode such that the rate D of 
dissipation of energy per unit volume is constant, that is 


D =0, & + 02 €2 = constant , (1) 


= | 


416 Tr-Cu1anc Hv and RicHAarp THORPE SHIELD ZAMP_ | 
N 


then the disc is of minimum volume. In condition (1), o, and a, are the principal 
stresses and ¢,, €, are the corresponding principal plastic strain-rates. It should 
be noted that although the middle-surface A of the disc may be only a part 
of the region A7, the special collapse mode must be defined throughout se | 
whole of A;. The collapse mode must satisfy condition (1) throughout the 
middle-surface A of the disc and elsewhere in A; the value of D must not | 
be greater than its value in A. 

In general, it is not known whether a design satisfying condition (1) exists, | 
and if it does exist, it may be very difficult to obtain. In order to supplement — 
the procedure for minimum-volume design, three theorems will be given. In 
some cases, the theorems can be used to obtain bounds on the value of the - 
minimum volume, while in other cases the theorems may serve to determine 
the minimum-volume design, even though the collapse mode satisfying (1) 
has not been obtained. 

A portion of the boundary of A; where all displacement components are 
prescribed to be zero (built-in support) will be referred to as a motion-free 
edge, and a portion of the boundary where the loads are prescribed to be zero 
will be called a stress-free edge. The following two theorems are concerned 
with changes in the middle-surface A; and its boundary, the loads on the 
original middle-surface and boundary remaining unchanged. Although the 
theorems may appear to be intuitively obvious, their statement and indication 
of proof seems worthwhile. 

Theorem 1: If a motion-free edge is moved outwards so that the middle-. 
surface A; is increased, the minimum volume will not decrease. 

The proof of this theorem follows directly from the fact that the minimum- 
volume design for the problem with the increased middle-surface furnishes a 
design for the problem with the original middle-surface. The material, if any, 
in the region outside the original middle-surface can only add to the volume 
of the design. The theorem is true regardless of whether loads act over the 
addition to the middle-surface or not. 

An obvious corollary is the converse: moving a motion-free edge inwards 
cannot increase the minimum volume. 

Theorem 2: If a stress-free edge is moved outwards, the minimum volume 
will not increase, provided that no loads are assumed to act over the addition 
to the middle-surface. 

For the minimum-volume design with the original middle-surface furnishes 
a design for the problem with the increased middle-surface. 

The converse of Theorem 2 is: if a stress-free edge is moved inwards so as 
to remove an unloaded portion of the middle-surface A;, the minimum volume 
will not decrease. 

Theorems 1 and 2 are the counterparts for minimum-volume design of the 
theorems on the collapse loads of bodies obtained from a known body by 
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changing motion-free or stress-free boundaries. (See [5], [6].) Although Theo- 
ems 1 and 2 (and also Theorem 3 below) are given here for the particular case 
of minimum-volume design in plane stress, they are obviously true also for 
three-dimensional structures, and they admit simple generalisations [7] to the 
more general problem of optimum design (see [8], for example). 

_ The following theorem is concerned with the continuation of designs 
‘across a motion-free boundary. 

_ Theorem 3: Let two middle-surfaces A, and A, have a portion L of their 
boundaries in common but otherwise have no points of intersection, and let 
I, h, respectively be the thickness distributions of minimum-volume designs 
for discs with middle-surfaces A,, A, under prescribed loads and with L as a 
motion-free boundary for each disc. Then, if the edge forces exerted by these 
designs at points of L are equal and opposite, the union of A, and hy is the 
‘minimum-volume design for a disc whose middle-surface is the union of A, and 
A, and which is acted upon by the same loads. 

_ This follows from the fact that any design for the problem with the com- 
“bined middle-surface provides designs for the discs with middle-surfaces A,, A, 
if the motion-free boundary L be replaced. Since no design with middle-surface 
‘A, is of less volume than the design /, and no design with middle-surface A, 
has less volume than the design /, the union of h, and h, must be the minimum- 
volume design for the combined middle-surface. 

As a simple application of Theorem 1 to obtain bounds on the value of the 
“minimum volume, we consider a square sandwich plate of core thickness H and 
‘side of length 2 R and loaded with a concentrated load P at the center of the 
‘plate. The edges of the plate are motion- free (built-in). To obtain a lower 
bound on the’ minimum volume of the face sheets, we move the motion-free 
edge inwards and consider a circular plate of radius R. The volume of the 
“face- sheet material of the minimum-volume design for the circular plate is 
“(see (91) 
5 ue (2) 


WV, = 0-333 ah 


fi 
where o, is the yield stress of the face-sheet material. The value (2) is a lower 
bound according to the converse of Theorem 1. 
To obtain an upper bound, we can move the motion-free edge outwards and 
‘consider a circular plate of radius V2 R, but the following procedure gives a 
Bioser bound. We construct a moment field in the plate which is in equilibrium 
with the load P at the center. Using polar coordinates with origin at the center, 
we assume that within a circle of radius b, where 0 <0 S R, the radial and 
’ ‘circumferential moments M, and M, are equal to the maximum moment 
M, = 0, Hh. Outside the circle, we assume that M, = —My, M,=0. At the 
“circle ry = b, the face-sheet thickness / is zero and the shear force is continuous. 


ee 


ye. 
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The volume of the face-sheet material required is found to be 


die Cols 


_ insta" 
It 


2 3 
Atle 4 2K 4oR | 


The volume V, has a minimum value when 0 is given by 


pee in tana = O74 Re 
Seis 8 


and the value of V, is then 


Pai 
V, = 0-433 (5)! 


This value is an upper bound on the face-sheet volume of the minimum-: 
volume design for the square plate. Comparison with the lower bound (2) shows. 
that the minimum-volume design has face-sheet volume 0-383 PR?/o, H+13%.. 


3. Basic Equations for Plane Stress 


It is convenient to use the lines of principal stress, Figure 2, as reference 
lines. These lines coincide with the lines of principal strain-rate for an isotropic 


Figure 2 
Principal directions. 


material. The radii of curvature of the lines of principal stress will be denoted 
by @, and @,, and we have 
1 Op 1 0p 


1 ero Q2 walla Ve (6) 


where ¢ is the angle between the first principal stress direction and the x- axis, 
and where 0/0s, and 0/ds, denote differentiation along the principal stress 
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directions. The radii 9,, 0, satisfy 


cal ae) gg pm ea 


and also we have 


0 ( 0 il eG 0 ( 0 il og 
OS, \ OSg 0, OS, OS a) Oe 0s. (8) 


In terms of the velocity components v, and v, along the principal directions, 
the principal strain-rates are given by 
O04 v2 Ov ae 20 (9) 


os + pee oo , 
i Os, Qi : OS» Q2 


The shear strain-rate ¢,, is zero so that we have 


ieee OV, V4 Dey 
ew eS a Os, Cy a oe (00) 


Equations (7) to (10) can be used to obtain 


Oe 1 o eae, 
iepra ts tt ponts > Zee) 


emia) alae) Mar al oe 


This equation is the compatibility condition on the principal strain-rates &, é>. 

The forces per unit length acting on elements of the disc perpendicular to 
the second and first principal directions are denoted by N, and N, respectively, 
and we have © 


N,=0,h, Ny=oyh. (12) 


The equations of equilibrium are 


ON, N,—N, ON. IN}, oN 
On Z b MOE 13 
Os, “s Q2 OSg A Q1 ( ) 


4. The Four Types of Solution 


The condition (1) for minimum-volume design imposes a restriction on the 
velocity field and this condition does not involve the thickness / of the disc. 
The form of the restriction depends upon the shape of the yield locus used. If 
Tresca’s yield condition is adopted, Figure 1, there are only four essentially 
different types of plastic flow which can occur. For example, plastic flow with 
stresses represented by corner A of Figure 1 differs in character from plastic 
states at D by only:a change in sign, and stress points on BC become stress 
points on EF if the subscripts 1,2 are interchanged. As typical of the four types 
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of plastic flow, we shall consider, in turn, stress states represented by corners | , 
and C, and points on the sides AB and BC, not including the end points. For 
convenience in the algebra, the constant in condition (1) will be taken through- 
out to be k oy, where k is a positive constant and gy is the yield stress. . 


Regime A: For the corner A, 0; = 0g = 0, and the principal strain-rates| 
are positive. The condition (1) then requires | 
et+e=k, a4) 

or | 
OUy OUy _ : 1 

athe Teac, (159) 


where u,, u, are the components of velocity referred to the x, y directions. 
The equations of equilibrium (13) become 


“=0, —=0. (16). 


Hence 4 = constant throughout the whole regime. Any multiply-connected 
region with uniform tension applied along its boundary admits the design of. 
a disc of constant thickness, and the result is the minimum-volume disc. A 
velocity field satisfying (14) with positive strain-rates is in this case the uniform 
dilatation 


1 1 
U,= Zk, “= Rhy (17) 
Regime C: For the corner C, o, = 0, and o, = — oy. The principal strain- 
rates must satisfy the inequalities 
—~4 2820, (18) 


and the condition of constant rate of dissipation requires 


—e=k. « (ESy 
The equilibrium equations (13) are now 
Oh, h 
~~ +—=0 
Os, in Qa ; (20) 
h 
—-—=0. 
a (21) 


From (21), @, = co when h is non-zero and the s, lines are therefore straight. 
Thus the orthogonal net of principal stress trajectories consists of families of 
straight lines and parallel curves. It is convenient at this point to introduce 
new coordinates, Figure 3. One of the curved principal stress lines is chosen 
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as a base curve and distance from this curve along the straight lines is denoted 
by «. The inclination y and the distance « are used as coordinates, and we have 


ds, =da, ds,=0,dp, O2=O+%, (22) 


Figure 3 


Coordinate system for regime C. 


where o(¢) is the radius of curvature of the base curve (for which « = 0). From 
(19) and (9) we have 


or 
v,=—kat+fy), (23) 


where f is a function of only. From (10), it now follows that 
v2 = f'(v) + e2.9(9) » (24) 


where g is a function of ¢ only. The functions f and g are not completely 
_arbitrary but must be such that inequalities (18) are satisfied. The equilibrium 
equation (20) can be written 
oh h 


Om + o+a Si) 
so that we have 
plies SAP) 
a o(p) +’ yes 


where g is a function of @ only. 
Regime AB: For stress states represented by points on the side AB of 
Figure 1, 0, = % and dy > 0, > 0. The strain-rate e, is zero and the condition 


; of constant energy dissipation rate then requires ¢, = k. 
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The equation of compatibility (11) becomes in this case 
0 (1 ofl ills Sonih 
a (a) © ae ge) eee ee 


Equations (7) and (26) give 
Ome ib 
_ O : 
OS» Eig - oj 


or 
Le Se (28); 
OS : 

Similarly we have . 
Lene (29) 
Os; 


From (28) and (29), we see that the principal stress lines form a Hencky- 

Prandtl net. 

The equations of equilibrium are equations (13) with N,=o,/ and! 
Oy hi WN 0 that is 

0 C AGe =e Oh (d9 — 0) h 

Os, (1 : i) a OS» * Q1 


=0. (30) 


Regime BC: For stress states represented by the side BC, o, — o, = op, and. 
— &, = &; = k from the normality condition and condition (1). The equation of 
compatibility (11) again reduces to (26), and it follows, as for regime AB, that 
the lines of principal stress form a Hencky-Prandtl net. This result was ob- 
tained earlier by PRAGER [3] by a different approach. 

The equations of equilibrium (13) with N, — N, = o, h become 
Ooh _ ) 


0 Oy h 
Os, (4 a Q2 a as, (08 és = Q1 


=0, (31) 


where o, — 6, = oy and 0 < o, < op. 


5. The Matching of Regimes 


As shown in the previous section, the four types of plastic flow possible 
under the Tresca yield criterion lead to four types of solution for minimum- 
volume design. However, a particular minimum-volume design may involve 
more than one of the four types of solution. We shall discuss briefly the pos- 
sible matching of solutions corresponding to two regimes. Let J” be the dividing 
line of two regimes. We choose a set of orthogonal coordinates «, B in the 
neighborhood of J" such that the line « = 0 coincides with I”. The velocity 
components in the «, 6 directions are denoted by w and », respectively. Then 
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0 0 
u v U Ou ov v (32) 


=> ——— = —— — =S —— ——> —«-— 
Tig ce aes ica Je ae ita ame) ee Saale Be 


where 9 is the radius of curvature of the curve J”. For simplicity, the coordinates 
a, B are assumed to be cartesian in the neighborhood of J" and the lines 6 = con- 
stant are assumed to be straight in deriving the expressions for ¢,, €g, and &, ». 
Tn order to avoid fracture of the material, we must have w and v continuous 
cross J’. From (32), it can be seen that ¢, must be continuous but e, and &, » 
may be discontinuous. The equilibrium conditions require the normal force N 
and the shear force T to be continuous across the dividing line. It is not necess- 
ary that the thickness / be continuous across I within the frame-work of the 
present theory. 

All possible combinations of matching two regimes can be obtained by 
considering regimes A, AB, B, BC, C, CD and D. For the purposes of this 
discussion, it is assumed that a regime consisting of points on a side of the 
yield hexagon does not include the end points of the side, e. g. regime AB does 
‘mot include the points A and B. It is not necessary to consider all possible 
‘combinations as, for example, the matching of regimes CD and D is the same 
as matching AB and A except for a change of sign. The twelve cases that it is 
‘necessary to consider are listed in Table 1 below. 


Table i 
ie a 
: Regimes Possibility | Condition which makes it impossible 
: AoAB {D2 CD) yes 
4 web 3(D20 yes 
a AD BCs) (DBC) no Equilibrium and strain-rate 
; A:C (D: B) no Equilibrium and strain-rate 
oe ACD (Dis AB) no Equilibrium 
- A:D no Equilibrium 
: AB:B (CD:C) yes 
AB: BC (CD: BC) no Equilibrium and strain-rate 
A Hip iGuen(G DB) no Equilibrium and strain-rate 
4 AB:CD no Equilibrium 
a B: BC--(C : BC} yes 
j BEC no Equilibrium and strain-rate 
si et pete) ee a 


The possibility of whether two regimes can meet or not-can be seen from 
an examination of the Mour’s circles for the forces and strain-rates for the 
in question. It is found that the condition that eg be continuous 


. 
* 
; 
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across the dividing line J" does not preclude the matching of any two of t 
regimes. If the thickness / is allowed to be zero on both sides on J’, then the 
equilibrium conditions across J" can always be satisfied and all matchings ar 
possible. However, the equilibrium conditions impose restrictions if the thick ; 
ness / is assumed to be non-zero on at least one side of J’. Table 1 indicate 
whether matching is possible with non-zero on at least one side of I’, anc 
when it is possible it can be shown that J’ must be a line of principal strain rate} 
In the case when the matching of two regimes is impossible, the table indicate 
that it is impossible because the equilibrium conditions alone cannot be satiss 
fied or that the equilibrium conditions and the continuity of strain rate e, 
cannot be satisfied simultaneously. I 


6. Plane Stress Examples 
6.1. Straight line support for a single force 


We suppose that material can be placed in the semi-infinite region x > 0 ané 
that the boundary x = 0 is a motion-free support, Figure 4. Let a single force F 
be applied anywhere in the right half plane and directed towards the right fom 


Figure 4 


Semi-infinite region, single force. 


convenience. If the force makes an angle y with the x-axis, then the minimum- 
volume design differs according as y is less than or greater than 7/4. For p <7/4, 
the minimum-volume design consists of a single bar in the direction of the 
force and connected to the supporting line as shown in Figure 4a. The bar has: 
finite cross-sectional area but has zero width. For y >a/4, the minimum- 
volume design consists of two bars from the point of application of the force 
to the support, each inclined at angles of amount z/4 with the x-axis as shown: 
in Figure 4b. For y < 2/4, the structure is in regime B. The velocity field 
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5a tisfying condition (1) is 
u,=k(1—tan*y)x, u,=—2k (tany) x, (33) 
where u, and u, are the velocity components in the x and y directions. For this 


v locity field e, = RandO S—e, Sk fory Sa/4. For y > a/4, the structure 
is in regime BC, and the velocity field satisfying condition (1) is 


“= 0, b= —2Rx. (34) 


y 


e Since the minimum-volume designs do not take the effect of buckling into 
account, reversal of the sense of the force does not change the minimum 
cture. 


a 


6.2. Point-support, two forces at an angle 


_ Here we consider the minimum-volume design for three forces, two of 
which are equal and directed to a point on the line of action of the third. 

In Figure 5a, two equal forces are applied at the points M, N and directed 
wards a point on the perpendicular bisector of MN. The forces are to be 
balanced by a vertical force at the point of support G, which also lies on the 
sisector of MN. It is assumed that the point of intersection of the lines of 


~ 


: 
: 
ac 


Figure 5a 
Two forces at an angle, point support at G. 


tion of the forces at M, N lies outside the triangle formed by the points 
SMN. The angle MGN will be denoted by 2 y. The form of the minimum- 
olume structure depends upon the value of the angle y and also on the 
irections of the forces at M, N. 

For y > z/4 and for forces inclined at less than 2/4 to the vertical, the 
design is indicated in Figure 5a, It consists of two quadrantal fans RGS and 


2 
4 
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HGJ and the four tangents MR, MH, NS, NJ. The continuous velocity fie 
of constant energy dissipation rate associated with the design is defined in th 
sector HGJ by (the velocity field for sector RGS is similar) 
3 Ze 

“=k, ug——2kr(0—F x), = 

where the polar coordinate system used is shown in the figure. For the qua te 
plane SG/, 


u,=ka +key, ly = hx — ky ,. (3¢ 


the quarter plane RGH is similar.) The restriction that the forces at M, N bt 
inclined at less than 7/4 to the vertical ensures that the bars SN and RM a 
in tension and the bars HM, JN are in compression. The arc HJ is of the samf 
thickness as the bar JN (or HM) while the arc RS is of the same thickness 1 
the bar SN (or RM). The thickness in the fans is given by 


| 
pees FE! 
| 


where P is the force in the arc HJ or RS. . 
When the forces at M, N are inclined to the vertical at an angle greates 
than x/4, the minimum-volume design is as shown in Figure 5b. The angll 


Figure 5b 
Two forces at an angle, point support at G. 


of inclination of the forces to the horizontal is denoted by y’. The minimum 


volume structure consists of two bars MH’, N J’ in compression and a fal 
H'GJ" of angle 2 y’. The velocity field of constant energy dissipation rate fo: 


| 
: 
| 
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the sector H’GJ’ is given by 


§ u,=—kr, y= —2hr (0-2), (38) 
id for the sector R’GS’ by 
u,=—kr, Hemi te br (O— 5). (39) 


ere a, is the angle of GS’ with the horizontal axis. For the quarter plane 
J’; 
a “—=kh%X+2ka,y, ts=—ky—Zkhayx, (40) 


y 
where GJ’ and GS’ are taken as the x and y axes. 


For y < 2/4, the minimum-volume design consists of the fan MGN only, 
e 5c, in the case when the forces at M, N are inclined at angles of less 


Figure 5c 
Twv forces at an angle, point support at G. 


than [2 to the lines GM, GN respectively. In this case the velocity field of 
constant energy dissipation rate for the sector MGN is again given by (38), and 


for the quarter plane NGS by 
u,=kut2kypy, u,=—2Zkypx—Aky, 
R. In the sector UGV we have 


(41) 


with a similar field for the quarter plane MG 
ies ade be up =2Akr (8-5), (42) 
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and in the sector VGS { 
u,=—Akr, uUp=—2khyr—(l1—2Ak(0—-y)7, (43) 


with a similar field for the sector RGU. The parameter A must lie between 0 
and 1. In order for the velocity field to be continuous, the parameter 4 and the 
magnitude f, of the angles VGS (or RGU) are related by 


Bi =A(n—2y) —2y. (44) 


For different values of the angle p, 4 can be adjusted so that f, is not negative. 
For 2/6 < y < 2/4, A can be taken to unity and for 0 < py < 2/6, A can be 
taken to be 1/2. Although for A = 1/2, the region UGV has an energy dissi- 
pation rate which is less than the constant value elsewhere in the plane, the 
structure is still of minimum-volume as there is no material placed in the 
region UGV (see [2]). 

The case when py < 2/4 and the directions of the forces at M, N lie between 
the perpendiculars to GM, GN, respectively, and the horizontal can be solved 
as in Figure 5b. 

The above discussion gives the solution to the problem in all cases when 
the lines of action of the three forces intersect in a point which lies outside the 
triangle MGN formed by their points of application. For structures composed 
of tiebars and struts, as considered by MICHELL [4], the solutions to the pro- 
blems considered so far in this section are similar. However, when the forces 
act so that the point of intersection of their lines of action lies within the triangle 
formed by their points of application, the solution for discs differs from that 
for tie-bars and struts. To take a specific example, in Figure 5d we consider 


eu 


OT 


Figure 5d 
Three symmetrical equal forces. 


three forces each of amount 2 T acting at the corners of an equilateral triangle 
and directed towards the center of the triangle. If the structure is restricted 
to consist of rods in tension or compression, the solution for minimum-volume 
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is any frame that consists entirely of struts in compression [4]; for example, 
the structure formed by equal rods along the sides of the triangle. The volume 
of the minimum-volume design is 6-93 T a/o9, where 2 a is the length of a side of 
the triangle. The formulation of the minimum-volume problem used here con- 
tains rods in tension and compression as special cases so that it is to be expected 
that a disc which carries the load can be found which is of volume not greater 
than 6-93 Ta/oy. Cox [10] has exhibited such a design. In Figure 5d, the 
circular arcs are each of radius 2 a and carry compression of amount T. The 
region between the arcs is of constant thickness T/2 0) 4 stressed to the yield 
; stress in compression. The total volume of this design is 6-60 T a/oy. 


6.3. Circular supporting region, single force 


Consider a rigid circular supporting region (motion-free edge) of radius a, 
Figure 6. A single force P is applied at a point A in the plane of the circle at 
a distance R from the center O of the circle. The angle of inclination of the 
force to the radial line OA is denoted by y, and for convenience we suppose that 
it is not greater than 2/2. The form of the minimum-volume design depends 
on the value of y and the ratio R/a. 

| If the force is vertically downward (y = 2/2), Figure 6a, the minimum- 
volume design is then a disc ABC bounded by two similar equiangular spirals 
of angle z/4 which intersect orthogonally at the point A of application of the 
force, as shown previously by MIcHELL [4]. The velocity field associated with 
the design is given by 


un 


u,=0, tg = —2hrIn(Z) (45) 
with polar coordinates (7, ) centered at O. The whole disc is at regime BC, 
the line of principal stress being equiangular spirals. The spirals which form 
the edges of the disc are two ribs, the upper carrying a tensile force P/ V2 and 
the other a compressive force of the same amount. The thickness of the disc 
can be found from the equilibrium equations (13) and N, — N, = oh (see 
below). Along the upper edge of the disc N, has the value — P/ (V2 7%) where 
‘ry is the radius of curvature of the upper edge, and at the lower edge N, has 
the value P|(V2 r,), where 7; is the radius of curvature of the lower edge. A 
_ graphical procedure for determining the thickness has been given by PRAGER[3]. 
However, the total volume of the disc can be obtained without determining the 
“thickness. For the rate at which work is done by the force P in the collapse 
“mode (45) is equal to the internal rate of dissipation of energy. Since the rate 
_of energy dissipation per unit volume has the constant valuek op), we have 


P2kRin(--) =kovVn (46) 
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where V,, is the volume of the disc. That is, 


2 PE in(=). 


Vie 


59 


The minimum-volume design when the direction of the force P is inclined 
at an angle less than /4 to the vertical is similar in form. With this restriction 
on the direction of the force, the upper rib is stressed in tension and the lower 
rib is in compression and the design is still associated with the collapse mode 
(45). The volume of the design is given by V,, cosy where V,,, is given by a | 
and yp is the inclination of the force to the horizontal. 

When the inclination w of the force to the horizontal is less than 7/4, the 
form of the design is as shown in Figure 6b or Figure 6c, depending on whether 
sin y is less than or greater than a|(V2 R). For siny less than a/ (V2 R), the 
line of action of the force meets the circle in a point G and the angle «, between 


siny >a//2R 


Figure 6 
Circular support, single force. 
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AG and the normal 0G produced to the circle is less than 7/4, Figure 6b. The 
minimum-volume design consists of the single bar AG in tension. It was shown 
in Section 6.1 that for the semi-infinite region to the right of the tangent GD, 
with GD as a rigid support, the minimum-volume design is the single bar AG. 
By moving the motion free edge GD outwards, it can be deformed into the 
circle and Theorem 1 of Section 2 then shows that the bar AG is the design 


for support on the circle. The volume of the design is 


e. [cosw — ( Es -- sin?y)” | 3 (48) 


0 


=) 12) 
V, = 


For siny greater than a/ (V2 R), the form of the design is shown in Figure 6c 
It consists of a bar AG in tension and a disc GBC bounded by equiangular 
spirals intersecting at right angles at the point G, where G is such that AG and 
the radial line OG intersect at an angle x/4. The upper edge of the disc is a 
rib which carries a tensile force P. The disc is at regime BC and the bar at 
regime B. The collapse mode associated with the design is given by (45) 
within the annular region bounded by the support circle and the concentric 
circle through G. The magnitude of the circumferential velocity on the circle 


through G is denoted by U and we have 


a 


U=2kbIn~ (49) 


where 0 is the length of OG, 
g b=/2 Rsiny. (50) 


a 


Between the circle through G and the tangent GD to the circle at G, the 
velocity field corresponds to a clockwise rigid rotation of amount U/®, that is 
u,=O, eae (51) 
and the strain-rates are zero in this region. Finally, to the right of the tangent 
GD, the velocity field consists of the rotation (51) with the superimposed 
velocity field 
j u,=0, Uu=—2kx, (52) 
the x-axis being directed along OG and the y-axis along GD. The existence of 
the continuous collapse mode described proves that the design is one of mini- 
mum-volume, as the dissipation rate is constant where material is placed and 
the dissipation rate does not have a larger value elsewhere (see [2]). The 


volume of the design is given by 


Via y2¥ b n(2) ip Po (53) 


Oo a Oo 
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where c is the length of AG, 
c =/2 Rsin(= —¥). 


It may be noted that Theorems 1 and 3 can also be used to prove that ing 
design in Figure 6c is the minimum-volume design. 

The minimum-volume designs obtained in this section can be shown to be 
unique by using the results of Reference 11. 

In order to determine the thickness distributions of the discs for the pro- 
blems indicated in Figures 6a and 6c, it is necessary to integrate the hyperbolic 
system of equations (13) in conjunction with the appropriate boundary con- 
ditions. When N, and N, have been determined, the thickness h of the disc is 
obtained from oy h = N, — N,. The integration is simplified by the fact that 
the principal stress lines, which are the characteristics of the equations, are 
systems of equiangular spirals. Referring to Figure 6a for definiteness, we 
introduce characteristic variables ¢, 7 by the relations 


6=4+y, r=RE", (55) 


where (7, 6) are polar coordinates, OA being the line 0 = 0. The coordinate & 
is constant along the second principal stress lines, of which AB is one, and 9 
is constant along first principal stress lines, of which AC is one. The radii of 
curvature @, and @, are given by 


e=—-V2r, o=/2r. (56) 
With equations (55) and (56), the equilibrium equations (13) can be written 
0 0 
pe" Ni) —71N2=0, whe Md tad CA eck! (57) 
It follows that 
ee 
Cer +f=0, (58) 
where 
f=orh=rN,—-r1N,. (59) 


We assume that the upper edge AB of the disc ABC is a rib which carries a 
tension T and that the lower edge AC is a rib carrying a compression S. Then 


for equilibrium, N, has the value —T///2 r along AB and Ny, has the value 


SV2r along AC, With these boundary conditions, the method of RIEMANN 


(see [12], for example) can be used to determine f throughout the disc. It is 
found that 


Oorh=flE, n) = sae 2) “Js (2 yén n) _ Vz se Ih (2 yén) 
1 p(2yén). (60) 


Le ey 
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Zusammenfassung 


Es wird die Dimensionierung auf minimales Volumen fiir Scheiben (ebene 
Spannungszustande) diskutiert, die aus starr-ideal-plastischem Material bestehen 
und der Fliessbedingung von TREsca geniigen. Dabei wird gezeigt, dass es vier 
-Typen von Lésungen gibt, und es werden mehrere Dimensionierungsbeispiele 


gegeben. 
(Received: November 23, 1960.) 


- 


ZAMP XII/28 


4 3 4 ' ZAMP 


Graphische Untersuchungen von erzwungenen Schwingungen 
mit unstetiger Geschwindigkeitsfunktion 


Von Maxtimitian MicHEeL Lawina, Ziirich?) 


1. Einleitung 


Das in dieser Arbeit behandelte System (vgl. Figur 1) besteht aus einer 
Masse m und einer Feder mit der Federkonstanten c. Die Auslenkung der Masse 


Y.~A0KRRR.WNX 


+q(t) 
Figur 1 
Schwinger mit Anschlag. 


m wird auf der linken Seite durch den Anschlag A, der sich im Abstand J von 
der Ruhelage der Masse befindet, begrenzt. Im allgemeinen Fall ist dieser Ab- 
stand g, nicht konstant, sondern eine Funktion der Zeit, die Federkonstante 
eine Funktion der Auslenkung, und an den Fiihrungsflachen zwischen Masse 
und Gehause kénnen nichtlineare Widerstandskrafte auf das System einwirken. 
In diesem allgemeinen Fall hat die Bewegungsgleichung unter einer Stérkraft 
P(t) die Form: 

mij + F(q, 4) = Pi). (1) 


Die zu dieser Gleichung gehérenden Randbedingungen bieten wesentliche 
mathematische Schwierigkeiten. Bezeichnen wir die Stosszahl zwischen der 
Masse m und dem Anschlag A nach dem n-ten Zykel mit k, (wobei k, < 1) 


1) Institut fiir allgemeine Elektrotechnik der ETH; jetzt Technische Hochschule Oberschlesien, 
Gliwice. 


Vol. XII,1961 Graphische Untersuchungen von erzwungenen Schwingungen 435 


und die Stossdauer mit f,,, so lauten die Randbedingungen fiir den m-ten Zykel: 


n—1 
fiir t= S'(T; + toi) ist g=—Go, = —FyGn-1- (2) 
i= 
ig bedeutet hier die Zeit zwischen 2 Anschlagen, ¢,; die i-te Stossdauerzeit, 
4,—1 die Geschwindigkeit, mit welcher die Masse nach n — 1-Zyklen gegen den 
Anschlag stésst. Die Gréssen T; und ¢,, sind in Figur 2 illustriert, welche den 
Zeitverlauf einer solchen Schwingung bei konstantem Wert fiir p zeigt. 


j i" 


Figur 2 
Zeitverlauf der erzwungenen Schwingung mit Anschlag. 


In der bisher vorhandenen Literatur sind dhnliche, einfachere Bewegungen 
-behandelt worden, von welchen hier zwei Beispiele angegeben werden sollen. 
KavuDERER [1, S. 426]?) behandelt einen autonomen Verband mit impuls- 
_artiger Erregung mit konstanter Energie. Er berechnet die Energie 4, des 
Systems (unter Beriicksichtigung der von der momentanen Geschwindigkeit 
_abhangenden Verluste) im m-ten Zykel und driickt diese durch die Energie des 
vorhergehenden Zykels 1 aus. So erhalt er eine Rekursionsformel, deren 
Konvergenz-Verhalten iiber die Existenz eines Grenzzykels Aufschluss gibt. 
Russakow und CHARKEWITSCH [2] lassen die Voraussetzung eines bekann- 
_ten und konstanten Energiebetrages der Erregung fallen. Die Erregerkraft 
stammt dort von einem Anschlag, der jeweils die mit der Geschwindigkeit ¢ 
anstossende Masse mit einer Geschwindigkeit — k-q abstésst. Die beiden 
Autoren gehen jedoch von anderen, vereinfachenden Voraussetzungen aus, 
_ namlich: 
a) Zwischen zwei aufeinanderfolgenden Stéssen treten keine Energieverluste auf. 
pb) Die Stérkraft ist sinusférmig. 
c) Die Stossdauer wird gleich Null gesetzt. 
Diese Arbeit behandelt im tibrigen nur Bewegungen unter der Bedingung, 
dass sie rein periodisch mit der Periode der Storkraft erfolgen. 


2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 446. 
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2. Problemstellung 


Den zwei oben besprochenen Systemen gegeniiber wollen wir dem in Figur 1 
dargestellten Schwinger etwas allgemeinere Eigenschaften zuschreiben. Von 
dem in [1] behandelten unterscheidet er sich dadurch, dass erstens eine Stér- 
kraft dauernd auf ihn einwirkt, welche erzwungene Schwingungen hervorruft. 
Die Stossimpulse, welche periodisch auftreten, haben keinen konstanten Ener- 
giebetrag, sondern ihre Energie ist jeweils proportional zur kinetischen Energie 
des Systems. Gegeniiber dem in [2] behandelten: Problem soll eine stationare 
Bewegung nicht vorausgesetzt werden, der ganze Einschwingvorgang soll hier 
untersucht werden. Dieser fiihrt nicht notwendigerweise zu stationaren Schwin- 
gungen ; vielmehr sollen auch nichtstationare Schwingungen, deren Amplitude 
jedoch in einem bestimmten Zeitabschnitt den vorgegebenen Amplituden- 
bereich nicht iiberschreitet, zugelassen sein. 

Diese Problemstellung ist auf die praktische Anwendung der Ergebnisse 
zurtickzufiihren ; das behandelte System soll némlich ein stossartig arbeitendes 
Werkzeug darstellen, etwa einen Abbauhammer, wie er im Berg- und Strassen- 
bau weiter Verwendung findet. Deshalb ist auch der Fall rein stationdrer 
Schwingungen weniger wichtig; weit gréssere Bedeutung kommt bei der 
Arbeitsweise eines solchen Werkzeugs dem Einschwingvorgang zu. Die Stér- 
kraft ist hier meistens sinusférmig, und in der Praxis haufig durch unausge- 
wuchtete rotierende Massen erzeugt. Der Abstand gq, soll als zeitlich verander- 
lich angenommen werden, wie das praktisch dann der Fall ist, wenn das Werk- 
zeug gegen das zu bearbeitende Material gedriickt wird. 

In [1] und [2] wird die Stossdauer gleich Null gesetzt. Versteht man aber 
unter Stossdauer die Zeit, wahrend der sich das Werkzeug mit dem zu bearbei- 
tenden Material in Beriihrung befindet, so macht die Stossdauer einen erheb- 
lichen Anteil der Schwingungsperiode aus und soll nicht als verschwindend 
klein angenommen werden. Physikalisch ist der Stoss ein Vorgang fiir sich 
selbst, der durch besondere dynamische Gleichungen beschrieben werden 
ratisste, wobei man eine der zahlreichen Hypothesen, die den gesamten Wider- 
stand des Stoffes in elastischen und nichtelastischen Widerstand aufteilen, 
zu Hilfe nehmen miisste. Um eine solche Komplikation zu vermeiden, betrach- 
ten wir die Stosswirkung rein phanomenologisch, indem wir annehmen, der 
Stoss wirke sich auf den Schwinger so aus, dass der Betrag seiner Geschwindig- 
keit nach dem Stoss gleich demjenigen vor dem Stoss, multipliziert mit der 
Stosszahl sei und sich das Vorzeichen der Geschwindigkeit beim Stoss umkehre. 
Die Stosswirkung lasst sich kaum durch eine einzige Konstante, die Stosszahl, 
beschreiben ; eine zweite Konstante, die Stossdauer, welche hier unentbehrlich 
zu sein scheint, muss eingefiihrt werden. In welcher Weise diese beiden Kon- 
stanten von der Art des bearbeiteten Stoffes, der Form des Werkzeuges und 
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er Stossgeschwindigkeit abhangen, kann nur durch besondere Versuche abge- 
lart werden. 

Bei einer solchen Bewegung werden die in [1] angewendeten Rekursions- 
formeln noch viel verwickelter, so dass man praktische Ergebnisse mit ihrer 
Hilfe nicht erwarten kann; die in [2] gebrauchte rechnerische Methode ist hier 
micht anwendbar, weil keine stationadre Bewegung vorausgesetzt werden soll. 
'Bevor wir eine fiir einen solchen Schwinger geeignete Methode wahlen, 
wollen wir noch kurz erwagen, zu welcher Bewegungsgruppe dieses System 
gehort. Da es in der Differentialgleichung (1) die Zeit explizit enthalt, ist es als 
heteronom zu bezeichnen. Gleichzeitig besitzt es aber auch manche autonome 
Eigenschaften; die Randbedingungen namlich, die auf die Koeffizienten der 
Bewegungsgleichung ihren Einfluss haben, sind von dem Bewegungsablaut 
nicht unabhingig, sondern hangen von der Geschwindigkeit am Ende jedes 
Zykels ab; man kénnte sagen, dass der Verband sich selbst gewissermassen das 
Gesetz des Fortschreitens seiner Bewegung in jedem Zykel von neuem vor- 
schreibe; dies ist aber eine Eigenschaft einer autonomen Bewegung. Dank 
dieser letzten Eigenschaft sollen zur Untersuchung des Systems die Methoden 
‘der Phasenebene angewendet werden. Da aber bei heteronomen Bewegungen 
in der Phasenebene eine Vieldeutigkeit auftritt, soll diese zuerst beseitigt 
werden. 


3. Einfiihrung des Winkelmasses als Zeitkoordinate 


Jede Bewegung kann bekanntlich durch eine Raumkurve im Koordinaten- 
system (g, y, ) dargestellt werden. Verzichtet man auf die ¢-Koordinate und 
betrachtet man nur das Lot der Kurve auf die (¢, g)-Ebene, so bekommt man 
eine Phasenkurve, die fiir jede qualitative Analyse der Bewegung ausreichend 
ist. Verschiedenen Anfangsbedingungen entsprechen verschiedene Phasenkur- 
ven. Bei autonomen Bewegungen entspricht nun der Gesamtheit der Anfangs- 
bedingungen eine volle Kurvenschar. Durch einen bestimmten Punkt der Pha- 
senebene (q, 7) kann nur eine einzige Phasenkurve gehen. 

- Bei nichtautnomen Bewegungen fallt diese Eindeutigkeit der Phasenebene 
weg, weil hier durch einen Punkt der Phasenebene mehrere Phasenkurven 
gehen kénnen. Dank der autonomen Eigenschaften des Systems ist es trotzdem 
zweckmassig, auch die hier behandelte Bewegung mit den Methoden der 
‘Phasenebene zu untersuchen. 

| Figur 3 zeigt den Phasenverlauf einer solchen nichtautonomen Bewegung 
‘mit einem Schnittpunkt. Die Zuordnung des Schnittpunktes A zu einer be- 
stimmten Phasenkurve a oder } kann dann eindeutig erfolgen, wenn wir in der 
Phasenebene auch die Zeitkoordinate einfiihren. Dann gehort der Punkt A mit 
fem Winkel « zur Kurve a, der Punkt A mit dem Zeitwinkel « + 2% zur 


‘Kurve 0. 


t 
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Ungedimpfte harmonische freie Schwingungen sind bekanntlich in der 
Phasenebene durch Ellipsen dargestellt. Verkleinert man den Ordinatenmass- 


Figur 3 
Phasenbild einer nichtautonomen Bewegung. 


Figur 4 
Addieren zweier Bewegungen auf der Phasenebene mittels Radiusvektoren. 
a = Bewegung 1, b = Bewegung 2, c = Bewegung 1+ 2, 
MaBstab der Ordinatenachse: 1:q. 


stab w-mal (w = Eigenkreisfrequenz), so wird die Ellipse zum Kreis. Fiir einen 
beliebigen Winkel « gilt (vgl. Figur 4): 


q _ Asin(wt+ ¢) 


paar INN) A cos(wt + 9) = tg(wi+ 9) 


und 


a=wt+o. 


Bei freien harmonischen Schwingungen ist also der Winkel «, unter Beriick- 
sichtigung des Anfangwinkels @, direkt der Zeit proportional. 
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Bei harmonisch erregten, erzwungenen Schwingungen, die nach der 
Gleichung = 
= mg+cq=Posinpt (3) 


verlaufen, setzt sich die Lésung wie folgt aus zwei Gliedern zusammen : 
q= [fi (A, B, Po) sinwt + f(A, B, Py) cosm t] + fs(Po) sinp t . (4) 


Das erste Glied, das heisst die eckige Klammer in Gleichung (4) stellt be- 
kanntlich eine freie, mit der Kreisfrequenz verlaufende Schwingung dar. Der 
zweite Bestandteil ist eine rein erzwungene Schwingung von der Kreisfrequenz 


p. Fiir das nachstehend dargestellte graphische Verfahren ist ausserdem wesent- 


lich, dass das erste Glied von den Anfangsbedingungen abhangt, das zweite 
dagegen davon unabhangig ist. Die Ausdriicke fiir die Koeffizienten f, und f, 
enthalten also allgemein die Konstanten A und B, obwohl, wie man aus 
Gleichung (5) sieht, auch noch ein von den Anfangsbedingungen formal unab- 
hangiger Ausdruck aus dem ersten Glied abgespalten werden kann. 

Beide Glieder der Gleichung (4) kénnen in der Phasenebene als Kreise dar- 
gestellt und ein Winkel « als Zeitmass angewendet werden. Soll nun derselbe 
ZeitmaBstab fiir die beiden Teilbewegungen in Gleichung (4) gelten, so mtissen 
wir jedem Punkt der Phasenkurve des ersten Gliedes einen Punkt der zweiten 
Phasenkurve zuordnen. Ein solches Verfahren wird im nachsten Abschnitt 
eingehend beschrieben. 


4. Linearer Zwischenverlauf 


Wir nehmen an, dass zwischen zwel aufeinanderfolgenden Stdssen die Be- 
wegung nach Gleichung (3) verlauft. Der Zwischenverlauf sei also linear. Soll 
die Storkraft periodisch, aber nicht sinusférmig sein, dann diirfen wir, wegen 
der Linearitat der Gleichung, die Stérkraft in eine Fourier-Reihe entwickeln, 
und die Wirkung aller sinusférmigen Glieder iiberlagern; daher geniigt es, hier 
die Erwagungen nur auf eine sinusformige Kraft, etwa nach Gleichung (3) zu 


- beschranken. 


Pere een ye VA ey 


Setzen wir in Gleichung (3) allgemeine Anfangsbedingungen ein (fiir ¢ = 0 
ist ¢ = Yo J = 4), SO wird 


q= (+ ~ PoP _ ) sinw t + dy cosw i + A sinpt. | (6) 


o w (c —m p*) c 


Zur Behandlung des in Figur 1 dargestellten Schwingers verzichten wir auf 
kontinuierliches Messen der Zeit, wie dies in den Gleichungen (2) ausgedriickt 
wurde; vielmehr wird jetzt bei jedem Zykelanfang die Zeit erneut gleich Null 
gesetzt ; der Anfangswert von q ist als — 4o gegeben (vel. Figur 1); die Anfangs- 
geschwindigkeit Go ist durch die Endgeschwindigkeit im vorhergehenden Zykel 


bestimmt. 


; : 
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Gleichung (5) stellt zwei Bewegungen dar, von welchen die erste, in eckigen 
Klammern, mit der Frequenz w, die zweite mit der Frequenz # verlauft. Wahlt 
man den MaBstab der Ordinatenachse 1: w, so kann die erste Bewegung durch 
einen Kreis dargestellt werden. Dieser andert aber in jedem Zykel seinen 
Radius, weil die Bewegungsamplitude von den Anfangsbedingungen abhiangt. 
Das Winkelmass der Zeit « = w t + o gilt aber fiir jeden Zykel, da die Frequenz 
w von den Anfangsbedingungen unabhangig ist. Wahlt man dagegen den Mass- 
stab der Ordinatenachse 1: #, so kann die zweite Bewegung durch einen Kreis 
dargestellt werden. Sein Radius und auch das Zeitmass « = pt bleiben wah- 
rend der ganzen Bewegung unverdndert. ; : 

Differenziert man Gleichung (5), so ist leicht ersichtlich, dass auch die Ge- : 
schwindigkeit die Summe der Geschwindigkeiten der beiden Teilbewegungen 
ist. Beschreiben wir die Lage eines Punktes zu einem bestimmten Zeitpunkt ¢ 
in der Phasenebene mit einem Radiusvektor r, oder rz, so gehdrt zum entspre-_ 
chenden Punkt der Summe beider Bewegungen ein Radiusvektor 


r=frfj+f,. 


Dabei muss man darauf achten, dass im allgemeinen die Ordinaten-MaBstabe 
beider Radiusvektoren verschieden sind (w + #), dass also nur eine Bewegung 
durch einen Kreis, die andere dagegen durch eine Ellipse dargestellt ist. Um 
den gesuchten Radiusvektor r, zu erhalten, muss man zuerst den zur Ellipse 
affinen Kreis aufzeichnen, auf diesem den entsprechenden Zeitwinkel pi auf- 
tragen und dann das Lot auf die Ellipse fallen, wie dies Figur 4 veranschaulicht. 

Selbstverstandlich lasst sich umgekehrt auch der Radiusvektor r, bestim- 
men, wenn die Vektoren r und r, bekannt sind. Auf Grund der bisherigen 
Uberlegungen ist es nun méglich, ein Phasenbild der durch die Gleichung (5) 
beschriebenen Bewegung aufzuzeichnen. Dabei kénnen sich die Anfangsbedin- 
gungen infolge der Stdsse bei jedem Zykel andern. Diese Konstruktion soll im 
folgenden fiir zwei verschiedene Fille, namlich w = f (Resonanz-Fall) und 
w + p (keine Resonanz), durchgefiihrt werden. 


Fall a): p + @ 


Die rechte Seite der Gleichung 5 wird fiir beliebig angenommene Anfangs- 
werte gy und q, durch einen Kreis (Bewegung 1) und eine Ellipse (Bewegung 2) 
dargestellt. Der MaBstab der Ordinatenachse ist dann 1: «. Der Kreisradius 


wird 
q yp _ 2. 
V( al err + Qo: 


da aber die Werte g, und q, fiir jeden Zykelanfang andern, werden auch die 
Kreisradien fiir jeden Zykel verschieden. Die Ellipse hingegen, die das zweite 
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Glied der Gleichung (5) darstellt, bleibt unverandert, die waagrechte Halb- 
achse wird 
mele 1 
c—mp*’ 

die senkrechte 
fe a? 


c—mp? wo 


In Figur 5 ist die Konstruktion des Phasenbildes der Bewegung g(t) veran- 
schaulicht. 


-——-99 


Figur 5 
Konstruktion der Phasenkurve fiir @ + p, a = Phasenkreis der Bewegung 1 fiir den ersten Zykel, 
b = Phasenellipse der Bewegung 2, c = der mit ihr affine Kreis, d = Phasenkurve der Bewegung 
1 + 2, e = Phasenkreis der Bewegung 1 fiir den zweiten Zykel. 


Der ausgezogene Kreis a stellt das erste Glied der Gleichung (5) dar, 
namlich denjenigen Teil der Bewegung, der von den Anfangsbedingungen 
abhingt. Die Ellipse 0 stellt das zweite Glied der Gleichung (5) dar; sie ist also 
von denAnfangsbedingungen unabhangig. 

Der Kurvenzweig von 1... 2 gibt die Konstruktion fiir einen Zeitabschnitt 
‘At an. 1’ ist der betreffende Punic der ersten Teilbewegung (« = w#,), 1” der 
entsprechende Punkt der zweiten Teilbewegung (« = #¢,). Er wird durch Pro- 
jektion des mit dem Zeitwinkel  ¢, auf dem affinen Kreis gelegenen Punktes 
=! die Ellipse bestimmt. Wahrend der Zeit At wandert der Punkt 1’ nach a4 


ane y ea Nee 


are | 
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der Punkt 1” nach 2” und der die Gesamtbewegung beschreibende Punkt 1” 
nach 2”. Fahrt man in dieser Weise fort, so kommt man schliesslich zum | 
Punkt 3” der Gesamtbewegung, also zum Anschlag. 

Den ersten Punkt 4” der Phasenkurve des zweiten Zykels findet man leichf 
durch Multiplikation von v mit — k unter Beibehaltung des Wertes von q 
(namlich g = —q). Es gilt nun, die zu 4” gehérenden Punkte der beiden 
einzelnen Teilbewegungen zu bestimmen. Der Punkt auf der Ellipse ist leicht 
zu finden, indem man auf den affinen Kreis von A aus (A = letzter Punkt 
unmittelbar vor dem Anschlag) den Zeitwinkel p ty (fp = Stossdauer) abtragt.. 


Beriicksichtigt man noch, dass die Summe der beiden Vektoren 04” und 04” 


den resultierenden Vektor 04” ergeben muss, so findet man auch den Punkt 4, 
mithin den Radius des Kreises fiir den zweiten Zykel. Von 4” aus geht dann, 
die Konstruktion in der gleichen Art weiter, wie sie in Punkt 1 begonnen wurde. 


Fall b): p=o 
Schreibt man Gleichung (5) in der Form: 


psinwt—wawsinpt P 
w (c — m p?) a2 


q= (£ sinw t + gy coset) 


so ist ersichtlich, dass bei 

c 
Pos we 

der Koeffizient von P, einen unbestimmten Wert 0/0 annimmt. Mittels der 
Regel von L’Hospirat findet man 


q= (£ sin@ t + dp COS@ i) + (. a) <inw t — ee cee cos@ t) , (6) 
o Pa j 


2m wo 2 


Die durch Gleichung (6) beschriebene Bewegung kann man sich wiederum 
aus zwei einzelnen Teilbewegungen (die beiden Ausdriicke in Klammern) zu- 
sammengesetzt denken. Figur 6 zeigt den Phasenverlauf einer solchen Bewe: 
gung. Die erste Teilbewegung ist durch einen Kreis mit dem Radius 


arom, Wed Nae 

O Ay= a V4) + 
dargestellt, die zweite durch eine Spirale a, bei welcher der Radius nach de: 
Drehung um den Winkel w At seinen Wert von 


auf 


r (a+ At) = 5% Vit @+odep 
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~90 
Figur 6 


Konstruktion der Phasenkurve fiir @ = p, 1 = Phasenkreis der Bewegung 1 fiir den ersten Zykel, 
a = Phasenkurve der Bewegung 2, b = Phasenkurve der Bewegung 1 + 2. 


andert. Wegen = w ist das Zeitmass fiir beide Bewegungen gleich, das Pha- 
senbild der resultierenden Bewegung ist durch ein rotierendes Parallelogramm 
gebildet, dassen Seitenpaar OB, und A,C, seine Lange mit dem Drehwinkel 

_ regelmassig Andert. Figur 7 zeigt, wie man die momentane Lange dieser Seiten 
graphisch ermitteln kann. 


Figur 7 
Zeichnerische Bestimmung der Radienfolge OBy, OB,, ... fiir die Phasenkurve a, (Figur 6). 


5. Nichtlinearer Zwischenverlauf 


Das oben beschriebene Verfahren kann auch fiir den Fall erweitert werden, 
dass der Bewegungsverlauf zwischen zwei Stdssen nichtlinear ist und durch 


die Gleichung 
mG +cq+mF(g,4,t) =0 (7) 
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beschrieben wird. Setzt man analog des linearen Zwischenverlaufs 


so erhalt man © 
gt+or'qt+ F¢gt) =90. 


In dieser Form eignet sich die Gleichung gut fiir die Untersuchung mittels 
der sogenannten Delta-Methode [3]. Schreibt man Gleichung (7) in der Form 


ag +orq+ FGd)=0, 

so folgt 
dj _ —otq~ Fy 4.) * 
dq q ; 8) 
Wegen der Nichtlinearitat kénnen weder die Phasenkurven der Teilbewe- — 
gungen als Kreise dargestellt werden, noch gibt es einen zur Zeit proportio- | 
nalen Winkel «. Fiir kleine Zeitabstande At kénnen diese beiden Hilfsmittel 
jedoch wieder verwendet werden. Anstatt wie im letzten Beispiel den Ordina- _ 
tenmaBstab 1: zu wahlen, soll hier eine neue Variable 


Ae ae 


@ 

eingefiihrt werden. Aus Gleichung (8) folgt dann 
gues. 
Pen (9) 

aq y 
Gleichung (9) stellt die Differentialgleichung der Phasenkurve im y-g-Koordi-_ 
natensystem dar. Sind beliebige Anfangswerte tp, 7) und y gegeben, so kénnen 
diese in Gleichung (9) eingesetzt werden, und es ergibt sich: 


i 1 Fg a 

dy 0) = — qo + Tah (Jor Vor to) Jo + or Fy VR 

dq Yo ns Vo ' (0) 
Dieser Ausdruck gibt die Neigung der Phasenkurve im Punkt A(q, Yo) an 
(vgl. Figur 8). Wir ersetzen nun den Phasenkurven-Abschnitt in der Nahe 
von A durch einen Kreisbogen mit dem Radius 7). Auf diesem Kreisbogen 
gelangt man zum Punkt A,, fiir welchen die neuen Werte ¢ = t, + At, g= 1, 
p= i einzusetzen sind, um den neuen Radius 7, zu bestimmen. Gelangt man 
mit dem ersten Kurvenzweig zum Anschlagpunkt A,, so findet man den 
ersten Punkt des folgenden Zykels auf die bekannte Weise durch Multiplikation 
des Wertes y, mit — k. Die Stossdauer wird dadurch beriicksichtigt, dass bei | 


der Berechnung des dem Punkte A, entsprechenden Abschnittes 1/w? F’ die 
Zeit t = ¢, + At, eingesetzt wird. ; 
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- 
Figur 8 
Konstruktion der Phasenkurve bei nichtlinearem Zwischenverlauf mittels der Delta-Methode. 


6. Schlussbetrachtungen 


Untersucht man mit der hier beschriebenen Methode Schwinger, deren 
Parameter konstruktiv mégliche Werte annehmen, so sieht man, dass sich das 
Phasenbild der Bewegung schon nach einigen Zykeln, unabhangig von den 
Anfangsbedingungen einem Grenzzykel oder wenigstens einem schmalen strei- 
fenférmigen Grenzbereich nahert, wie dies Figur 9 zeigt. Dies bedeutet, dass 
eine stationare (oder quasistationare) Bewegung zu erwarten ist. Wollen wir 
aber aus gewissen, zum Beispiel energetischen Griinden, ganz andere Werte 


Figur 9 
Entstehen eines Grenzzykels (d). 
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an 


der Parameter wihlen, so miissen wir uns, anhand der oben dargestellten — 
Methode, zuerst iiberzeugen, ob ein solcher Schwinger iiberhaupt stationare — 


Bewegungen ausiiben kann und welche Amplituden, Stossgeschwindigkeiten 
usw. diese aufweisen. Ein unbestrittener Nachteil der hier besprochenen 
Methode liegt darin, dass man umgekehrt handeln muss; zuerst gilt es, alle 
Parameter zu wahlen, und dann erst kann man (auf ziemlich komplizierte 
Weise) das Bestehen eines Grenzzykels untersuchen. Es ist jedoch nicht gesagt, 
dass sich nach eingehender Bearbeitung diese Methode nicht umkehren lasst und 
dann zu einem geometrisch-analytischen Verfahren wird. Dies wiirde erlauben, 
vom gewahlten Grenzzykel ausgehend, die geeignetste Eigenfrequenz des Systems 
zu bestimmen. In diesern Falle wird dann auch eine Diskussion der optimalen 
Parameter-Werte méglich. Dies iiberschreitet aber den Rahmen dieser Arbeit. 

Auf alle Falle stehen hier noch weitere Diskussionsméglichkeiten offen. Zum 
Beispiel kann bei der besprochenen Methode auch die Anderung des Abstandes 
go beriicksichtigt werden, indem man die Begrenzungsgerade bei jedem Zykel 
entsprechend verschiebt ; weiter k6nnen auch andere Parameter — wie Stosszahl 
und Stossdauer — und deren Einfluss auf die Bildung des Grenzzykels in der 
Diskussion mitberiicksichtigt werden. Auf dem rein rechnerischen Wege ware 
dies wohl kaum ausfiihrbar. 
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Summary 


An oscillator which is being struck exhibits non-linear characteristics on account 
of the instabilities in speed induced. Through a wide range of parameters the 
processes of motion brought about by a stimulating force lead to a stationary 
condition via a transient phase (limit cycle). The law of motion for each periodic 
movement between two successive stimulations is basically known. There is both 
a free and a forced oscillation, and these are superimposed. Each can be illustrated 
as a vector rotating in the phase plane (4, g), the point of which describes a phase 
curve. The resulting motion is described by a vector which is the sum of the two 
partial vectors. The author gives here a graphic method which he has developed 
from this. By establishing the connection between the angles described by the 
partial vectors and the time taken, the size and direction of the resulting vector at 
any moment is determined. The amount of the vector representing the forced 
oscillation remains constant during the whole process, while the vector describing 
the free oscillation changes its size and phase in each cycle. 

The author then extends his method to non-linear periodic movement by means 
of the so-called delta method. 


(Eingegangen: 9. Januar 1961.) 
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Constitutive Equations Involving Functional Dependence 
of One Vector on Another’) 


By Rona p S. Rrvtin, Providence, R. I., USA?) 


1. Introduction 


There are many physical situations for which the applicable constitutive 
equation involves the relation between the value of one vector at the instant of 
measurement and the values of another vector at times up to and including the 
instant of measurement. Examples of such situations are the relation between 
the magnetic intensity and magnetic induction fields in a hysteretic magnetic 
material and the relation between the electric displacement and electric fields in a 
‘lossy’ dielectric. In this paper we discuss how restrictions may be placed on the 
form of such constitutive equations as a result of the material symmetry. 

Although the method described may be employed for any type of material 
symmetry, it is applied particularly to the cases when the material is holohedral 
isotropic or hemihedral isotropic. 

The procedure used is similar to that used in previous papers [1], [2], [3]°). 
‘However, the problem discussed in the present paper is less involved from a mathe- 
matical standpoint than those discussed in the previous papers and the results 
obtained are simpler. 


2. The Basic Problem 


We suppose that a physical phenomenon taking place in a material is character- 
ized by the dependence of a vector v at time ¢ on the values taken by a vector u(r) 
at all times t in the intervall 0 <1 <¢# and on the values of u(r) and its time 
derivatives d*u/dr% (« = 1, 2,..., m) at the time ?. Let v,(t) and u,(t) be the compo- 

‘nents of v(z) and u(r) respectively in a rectangular Cartesian coordinate system ~. 
Then, we may write our initial assumption in the form 


t 
v; = &; |u,(t); | (2.1) 
t=0 
where the notation 
d*u 
ao) (C2) aie = y 

v, = 0, (2) , u, = u,(t) = 4; ne ae (hime ty 2 cr fA) (252) 
has been used and, in (2.1), a= 0,1, 2,...,4u- @, isa functional of w,(t) over the 
range 0 St S? and a single-valued function of ue (a. 5-0, 1,2, ..2 ,#)- For 


1) The work reported in this paper was carried out under a grant from the National Science 


Foundation. 
2) Brown University. 
3) Numbers in brackets refer to References, page 452. 
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t 
simplicity, we shall assume that the single-valued dependence on ul) is either | 
polynomial in character, or may be approximated with any desired degree of . 
accuracy by polynomial dependence. We also assume that the functional depend- | 
ence is continuous. | 
We now consider that the symmetry of the material is characterized by the . 
group of transformations {S}, of which S (= || S;, ||) is a generic transformation. 
Let ¥ be a rectangular Cartesian coordinate system related to x by 


%. = S: 


a ag 


a> Sez Sip = OR, (2. 3) 


and let v,(t) and u;(t) be the components in the system % of the vectors v(t) and 
u(t) respectively. Then, 


U,(t) = S45 ,;(t) and u,(t) = 5,3 u;(t) . (2.4) 


Since (2.1) is form-invariant under the transformation S, we readily see that, with _ 
abbreviations analogous to (2.2), ®; must satisfy the relation 


t 
p(t ; a) = $5 D; tn us 4 (2.5) 
TE t=0 


Let w be an arbitrary vector with components y, and y; Tespectively in the coor- 
dinate systems 7 and ¥. Then, e 
i (2.6) 
From (2.5) and (2.6), we readily obtain 


t 


t 
y,; ®; gt : ap =», ©; 
14 


u,,(T) ; “ = F(say) . (2.7) 
t=0 


=0 


Since (2.7) is valid for all transformations of the group {S}, F must be a scalar 
invariant of w, u(t) (0 <t <2) and wu (x = 0,1,..., m) under the group {S}. 
u') denotes, of course, the vector with components 4{) in the coordinate system x. 


In the next section we shall show how the restrictions imposed on @, by the rela- 
tion (2.7) may be obtained in more explicit form. 


3. Method of Solution of Problem 


We may regard ®, as a polynomial in ule ( 0, 1,..., mw), the coefficients in 
which are continuous functionals of u,(t) (0 St St). Denoting the Fourier half- 
range coefficients of u,(t) by a” (R = 0, 1,...), each of these functionals may be 
expressed as a polynomial Q') (U Na Ue ary ) which tends uniformly to the 
functional as N -> oo. ©; may therefore be expressed to any desired degree of 


approximation as a polynomial ps) in ue (R=0,1,...,N) and uy (a 0, di unis 
which tends uniformly to ®, as N +> oo. Thus, 
t 
: (a)} ~ p(N) 0) 1 N 0 1 
®; bo “p mw PS Coy : Uy anes ‘ Ul ) uy ul), sat ust), (31) 


where Pe ) tends uniformly to ©, as N + oo. 
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The Fourier half-range coefficients ue are, of course, given by 


t 
(R=. 2 Raut ih 
U; neg [yee cos dt (P>0O) and um ae Exo Cae (3.2) 
0 


Denoting the Fourier half-range coefficents of u,(t) by gs (c= OF eerere) 
we have, from (2.4), 
7(R R 
OP le sn (3. 3) 
By an argument similar to that leading to (3.1), we obtain 


©, itn wp] w PY (7, 7, 0), a, a, al). (3.4) 
Introducing (3.1) and (3.4) into (2.7), we obtain 
+B, P™ (T, 0, ..., 7), WO, a, ..., a) 
= y, P™ (u®, UM, ..., 0), uw, =F. (3.3) 


Since this relation is valid for all transformations of the group {S}, it follows that 
F is a scalar invariant under the group {S} of the N+ y+ 3 vectors wy, u 
Pai) 1,0, 2), and u®) (R= 0,1,...,.N)4). F may thus be expressed as a 
polynomial in the elements of an integrity basis for these vectors under the group 
{S}. Since F is linear in the vector y, it must be linear in those elements of the 


integrity basis which involve y. Let us denote these by J, (« = 1, 2,..., v) and the 
remaining elements which do not involve w by I, (« = 1, 2,..., %). Then, 
Fe) Ayyy: (3.6) 
y=1 

where A, is a polynomial in I, (« = 1, 2,..., %). From (3.6) and the definition of 
F, we obtain 

. = OJ, 
: | Pi a (3.7) 
We may obtain our result in terms of u,(r) rather than in terms of the Fourier 
half-range coefficients igs! (R= 0,1,..., N) by using (3.2) to substitute for oe 
‘in the expressions for J, (« = 1, 2,..., %) AIAG fe (Ole = Latter eters, Ds 


; In the next two sections, we shall carry out this procedure in detail in the 
case when the material is holohedral isotropic or hemihedral isotropic prior to t = 0. 


4. Results for Holohedral Isotropic Materials 


For a holohedral isotropic material, the transformation group {S} characterizing 
‘the symmetry of the material is the full orthogonal group. For this group an 


(R) 


4) u'®) denotes the vector with components U Ps in the coordinate system %. 
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integrity basis for the vectors YJ, u™ and U'*) is formed by 


yy, (4.1) 
yp: ul wp: ui) (4. 2) 

and 
al)» a4!) : ul®). uy : uw. v®), (4.3) 


We therefore take the expressions (4.2) as the invariants J, in (3.7) and the ex- 
pressions (4.3) as the invariants J,. With (3. 7) and (2.1), we then obtain 


ig N 
va S14, d+ SB, UM, (4.4) 
a=0 R=0 


where A, and Bp are scalar polynomial functions of the quantities (4. 3). 


Employing the expressions (3.2) for Cy we see that Equation (4.4) may be © 


re-written in the form 


t i 
“ N rs a 
x («) 2 Br Rat Bo 
Siege A u(r) cos——— dr + = / u(t) dr, (4.5) 
a=0 R=1 0 
where Ania = 0)1, 2. ,,) aud Be t=, 1, ... -W) are polynomial functions of | 


the quantities (4.3). Hence, it may be written in the form 


t 
“ 2 
o= Aa ielter) wt, 2) u(x) de, (4.6) | 
a=0 
0 


where p(t, tT) is a single-valued function of its indicated arguments and a polynomial 
function of the invariants (4.3). Employing (3.2), we obtain 


R Coane Ra ies 
uU®). up — a | | u(t,) * u(t.) cos se cos —* dt, dt, (4.7) 
where C is 1, 2 or 4 depending on whether both, one or neither of the quantities 
Rand T are zero. Similarly, 


t 
G » 


ul). yl) — ra ul) - u(r) cog Wt, (4.8) | 


0 


where C= 1 if R=0 and C= 2 if R + 0. Thus, the coefficients A, and g(t 
. - > , 2 Tt 
in (4.6) must be expressible as polynomials in expressions of the form . eee, 


Bc te! 
yl) . i y(t, vt) u(r) dr and / ¥ O(t, tT, T)) u(t,) + u(t,) dt, dt,, (4.9) 
0 0 0 


where y(é, Tt) and A(t, T,, T,) are single-valued functions of their indicated arguments. | 
In the particular case when the functionals ®, in (2.1) are hereditary functionals, . 


the kernels g(t, t) and y(t, t) must depend on ¢ and t through ¢ — t and O(t, T,, T.) 
must depend on #, t, and t, through ¢ — 1, and ¢ — 75. 
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We note also the particular case when the functional ®, is not particularly 
dependent on the values of u(r) and its time derivatives at time ¢. Then, in (4.6), 


we may take u\) — 0. We obtain the result 


t 


v = |v tT) u(t) dt, (4.10) 


0 


where y(t, t) is a polynomial in invariants of the form 


bt 
[| [% %, Ts) W(t) > U(t,) At, dt, - (4.11) 
0 0 


5. Results for Hemihedral Isotropic Materials 


_ For a hemihedral isotropic material, the appropriate transformation group {S} 
is the proper orthogonal group. For this group an integrity basis for the vectors 
w, u and u'®) is formed by the invariants (4.1), (4.2) and (4.3), together with 


and 
Tu, ul), u)] Z fu), wu), u'))] : [u®), u's), ul] , fu), u's), u?)] . (5.2) 


We therefore take the expressions (4.2) and (5.1) as the invariants J, in (3.7) and 
the expressions (4.3) and (5.2) as the invariants I,. Then, following a procedure 
similar to that employed in the previous section, we obtain 


t Ft 
ia rn 
v= S74, ul + | olt,2) ule) ar | [altri wl) x ule) den dr 
a= 0 6 0 0 
Me 


t 
a7] 
+ yf ate x) ul x u(r) dr + SY” Aap ul x ai?) = (523) 
=05 


a,B=0 


where Ag, Aag, p(t, t), x(t, Tr, T2) and p(t, tT) are polynomials in invariants of the 


forms (4.9) and 
t 
Cu, (6), 0?) i A(t, c) [u!, u!?), w(x)] de, 
0 
est 


[ [ow TT.) [w', u(t), u(t») ] dr, dre, (5.4) 
| 0 0 
and 


Fe ated: 
it a 4 oll, Tus tye Ze) [e0(ty), t2(ca), (c9)] dry dre ts - 
02677:0 


r attention to the case when the functionals are of the 


We again draw particula 
he case when they do not depend particularly on the 


hereditary type and to t 


~—ae 
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values of u (« = 0,1,...) at r=+#. The modifications which must be made to 
(5.3) and (5.4) in these cases are obvious. 
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Zusammenfassung 


Die vorliegende Arbeit befasst sich mit einer Klasse von Stoffgleichungen, die > 
so beschaffen ist, dass der momentane Wert eines Vektors vom gleichzeitigen und 
von allen vorhergehenden Werten eines zweiten Vektors abhangt. Es wird gezeigt, . 
wie man diejenigen Einschrankungen der Form solcher Stoffgleichungen, die von 
materiellen Symmetrieeigenschaften stammen, herleiten kann. Die Herleitung 
wird v6llig durchgefiihrt fiir den Sonderfall eines isotropen Materials mit oder> 
ohne einen Symmetriemittelpunkt. 


(Received: Mars 4, 1961.) 


A Solution of the Steady State Thermoelastic Equations 


By Witt1am E. Witiams, Liverpool, Great Britain?) 


Introduction 


It is well known [1]?) that certain classes of solution of the elastostatic equations 
may be derived from scalar potential functions. In the present note it is shown that 
there exists a particular solution of the equations of thermoelastic equilibrium, with 
nonuniform temperature variation, which can also be derived from a scalar potential 
function. This particular solution can be combined with the known solutions of the 
elastostatic equations to solve thermoelastic boundary value problems when the? 
shear stresses vanish on a plane. Boundary value problems of this type have been 
considered by McDoweELt and STERNBERG [2], SNEDDON and Lockett [3] andl 
SNEDDON and OLEs1ak [4]. The first named authors employ a GREEN’s function) 
approach whilst the other authors employ transform techniques. 

The solution obtained here is applied to the boundary value problems treated | 
by the above authors. It is found that the solutions may be expressed in a simpl 
form. It is shown, without elaborate calculation, that the stress field in a thick: 
slab with stress free boundaries is plane and parallel to the boundary [3]. 

The problem treated by SNEDDON and OLESIAK [4] of a penny shaped crack: 
on the surface of a semi-infinite solid is also considered. It is shown that the normal 
stress is of the form of a constant C multiplying a function which is independentt 


1) Department of Applied Mathematics, The University of Liverpool. 
*) Numbers in brackets refer to Reference, page 455. 
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of the prescribed values of temperature on the surface of the solid. The constant C 
is proportional to the total heat flux over the crack. This result was established 
by SNEDDON and OLEsIAk assuming the heat flux in the crack possessed a Fourier 
Bessel expansion. The present result is independent of any such assumption. 


Basic Equations 


We consider an elastic body in which there is established a temperature field 
T + &(x,y, 2) where T is the temperature in a state of zero stress and strain. 
%, y, z are the Cartesian coordinates of an arbitrary point within the solid and k 
will be taken to denote the unit vector along the z axis. It will be assumed that 
there are no body forces or heat sources within the solid. The equations of thermo- 
elastic equilibrium thus become [5] 


(1 — 2) V2u+4+ grad[d4 —2a(14+ 7) B)=0, (1) 
y270=0. (2) 


u is the displacement vector, A = div u, 4 is Poisson’s ratio and « represents the 
coefficient of expansion of the solid. 
The stress vector Z (cf. [6]) across an element of surface whose normal is 
parallel to k is given by 
Z [2nA4—2a(1+%) 9] 


aoe (1 — 27) 


"where pu is the modulus of rigidity. 
: It is well known that [1] a particular solution of equations (1) and (2) is 


Ou 
aU a + grad (u,k), (3) 


u=t,=gradF, #=0, (4) 
where F is a potential function. The corresponding stress vector is 
OF 
Ay 2 grad— : 


It may be verified by direct substitution that another solution of the above 
equations is 


0G 
u=u,=2 grad. ; (5) 
— 0?G 
i 3— 47 (6) 


2a(1+%m) 02?’ 


where G is also a potential function. The corresponding stress vector is 


0G 0?G 0?G 
Z= bb {grad je + 2 z grad = 2, Oz? Rl 
From equations (4) and (5) we see that a solution of the equilibrium equations 
with zero shear stress on z = 0 is 

: 


0G 1 
= = it : | 
u =U, = z grad Fy. 3 gradG (7) 


f @ is still given by equation (6). 
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: 


A more general solution with zero shear stress on z= 0 can be obtained by 
adding to wu, the solution uw, [1] where 


Ox Ox j 
tu, = (1 — 2m) grady + z grad] — (3-4) Gk | 
where x is a potential function. The temperature # associated with u, is zero. It i 
convenient for mani-pulative purposes to write G as [4 (1 — m) 4]/(3 — 4 m) and 47 
as — ¢/(3 — 4) — » where ¢ and y are potential functions. We then obtain as a1 


solution of equations (1) and (2). 


ad ad 
u = —grad¢ + 2 grad-5 +k ris 


wh 
0z 


+(3-4n 2k, (9) 


— (1— 27) grady — z grad 


(1— 7) 0° 
OT Oe ene (10)) 


z = 0 is independent of ¢ on z = 0. The stress vector corresponding to the display 
cement and temperature of equations (9) and (10) is 


The solution given by equations (9) and (10) has the advantage that (u, k) on 


Z=2u[z grad, (6—v) —k 55 -y). (11)) 


Application to Particular Problems 


The above solution can now be applied to the solution of particular boundary? 
value problems. We consider first the problem of the semi-infinite elastic solid] 
occupying the region z = 0 with the plane z = 0 stress free. It will be assumed that: 
either # or 09/0z is known on z= 0. & and the components of uw are required to 
tend to zero as z-> oo, 

The boundary value problem for # (and hence ¢) can be solved by Fourier or-> 
Hankel transform methods or by a GREEN’s function approach. The boundary’ 
condition on z= 90 gives 0?/dz*(¢— yp) =0 on z= 0, hence since ¢ and y are} 
harmonic functions 0?/dz? (¢ — py) = 0 and hence we obtain STERNBERG and Mc Do-- 
WELL’s result that Z = 0. In order that the conditions at infinity are satisfied we: 
must in fact have ¢ = wy and hence 


u = 2 (yn — 1) gradd + 4 (1 — ») eee: (12)) 


If ¢ is expressed as a double Fourier transform it can be shown that equation (12)) 
agrees with the displacement given by SNEDDON and Lockett [3]. We now con-- 
sider a similar problem for the thick elastic slab bounded by the planes z = 0, h,| 
these planes being stress free. Equation (12) once more represents the solution of! 
equations (1) and (2) which is stress free on z = 0, h provided that ¢ is now deter-- 
mined from the known values of # or 09/dz on the boundary surfaces. We thus: 
have once more that Z = 0. 

We can also employ equations (9) and (10) to give a simple solution to the: 
problem of a penny shaped crack on the surface of a semi-infinite solid. Axial! 
symmetry will be assumed and the solid will once more be assumed to occupy the: 


region z = 0. We shall consider the particular problem when the applied normal! 
stress on the crack is zero. . 


y 
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_ The boundary value problem for # will depend on whether the temperature or 


the heat flux is specified on the crack. In either case the boundary value problem 
_is of a type which is amenable to solution by established methods [4], [7]. If the 


radius of the crack is a and r is a radial coordinate in the plane z = 0 then the 


elastic boundary conditions give 
va @—y)=0, 087 Sa, (13) 


Oe tog V 2G (14) 


_ Equation (13) may be rewritten 


Ory (1 + 7) 

= <r< 

ye & (Lacs a) (Eos OS ete aX (15) 
The potential boundary value problem equations (14) and (15) is of a type whose 
solution is well known [7]. By direct manipulation of this known solution it can be 


_ shown that az Ree (+n) 
Vo San n n 
xD Spee fey erie eR (16) 
where 
° 
g(r, 2) =| cos u Jo“) e “du. 
0 


1 
Cis f 0 Q(e) de, where Q (r/a) is the flux within the crack. Thus 
0 


02 
= edt) Gio E gradg(r, z) —R sn | é (17) 
Equation (17) agrees with the result of OLESIAK and SNEDDON [4] who obtained it 
by assuming that Q possessed a Fourier-Bessel expansion. u can be expressed in 
terms of @ and g from equations (9) and (15). 
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Résumé 


Une solution particuliére des equations d’équilibrium thermoélastique est 
derivée d’une fonction potentielle. Cette solution est combinée avec les solutions 
connues des equations élastostatiques et appliquée a la solution d’une classe parti- 
culiére de problémes aux limites. 
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Réunion de Printemps de la Société Suisse de Physique, 
a Rolle les 6 et 7 Mai 1961 


Compte rendu des communications de physique et de mathématiques appliquées 


Abtastregelung fiir ein hochkonstantes Magnetfeld. Von B. WALTHARD, 
Fribourg?). 

Die hier beschriebene Stabilisierung wurde fiir ein doppelt fokussierendes Eisen- 
spektrometer fiir B-Teilchen von 20 keV bis 2,5 MeV gebaut [1], [2]2). Seine maxi- 
male Energieauflésung soll 10-4 betragen. Das fokussierende Magnetfeld H bzw. die 
Induktion B in der Vakuumkammer muss demnach beim Normalradius des Spek- 
trometers von 30cm von 30 bis 350 Gauss variiert werden kénnen, und zwar in 
Schritten von 2- 10-5. Ausserdem muss sein Betrag auf 2- 10-5 konstant sein. 


+300V 


Sampler 


Figur 1 


Figur 1 zeigt das Prinzipschema der dazu gebauten Regelung. Der Strom I zur 
Erregung der Feldspulen wird von einem stabilisierten Netzgerat geliefert. Er kann 
mittels 12 parallelgeschalteter Tetroden der Type 807 kontinuierlich von 0 bis 1 A 
varuert werden. Zwei Regelkreise sorgen nun fiir die Stabilitat des Magnetfeldes. 
Einerseits wird durch eine starke Wechselstromgegenkopplung dafiir gesorgt, dass 


1) Physik-Institut der Universitat. 
*) Die Zahlen in eckigen Klammern beziehen sich auf das Literaturverzeichnis, Seite 458. 
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der Strom I durch die Feldspulen nur sehr langsam von einem jeweils eingestellten 
Wert abweicht, namlich héchstens um 2- 10-5 pro 15 Sekunden. Er wird zu diesem 
Zweck iiber einem Shunt von 100 Q gemessen; das gewonnene Signal wird kapazitiv 
an den Eingang eines Gleichspannungsverstarkers gekoppelt. Verstarkt um den 
Faktor 10° werden damit wieder die Langsréhren ausgesteuert. Diese Anordnung 
stellt eine nahezu ideale Stromquelle dar. Sie ist ausserdem mit der unter dem 
Namen Miller-Integrator bekannten Schaltung identisch. Tatsdchlich ist der 
Strom I mit grosser Genauigkeit dem Integral einer Steuerspannung U, proportio- 
nal, die man allenfalls zusatzlich an den Eingang des Gleichspannungsverstarkers 
anlegt. 

Ein zweiter Regelkreis hat nun die bei einem Eisenmagneten unvermeidlichen 
‘Anderungen des Zusammenhanges zwischen Strom und Magnetfeld auszuregulieren 
sowie auch den Nullpunktfehler des Integrators. Das dazu notwendige MeBsignal 
fiir das Feld wird durch eine im Spektrometer rotierende Spule H, erzeugt. Als 
Referenz fiir die Regelung dient das Feld eines Dauermagneten, in welchem eine 
weitere, auf derselben Antriebswelle rotierende Spule Hp angebracht ist. Die Welle 
wird von einem Synchronmotor mit 4500 Umdrehungen pro Minute angetrieben. 
Die beiden induzierten Wechselspannungen werden nun gleichgerichtet und geglat- 
tet. Es erfolgt eine Spannungsteilung des Referenzsignals, dem gewiinschten Wert 
des Feldes entsprechend, und schliesslich Bildung der Differenz durch Serieschal- 
tung. Zur Gewdhrleistung einer hinreichenden Langzeitstabilitat wird die ganze 
Mess- und Vergleichsanordnung auf 0,1°C thermostiert. Das Differenzsignal AU 
betragt bei der zulassigen Abweichung des Feldes von seinem Sollwert mindestens 
100 wV. Es wird nun iiber einen Abtaster (Sampler) an den Eingang des inte grie- 
trenden Stromreglers gefiihrt. Der Abtaster ist ein mechanischer Schalter, der alle 

6s periodisch den Regelkreis fiir die Dauer von 1 s schliesst. 


; 


rat = 


: Figur 2 


In Figur 2 ist ein méglicher Verlauf des Fehlersignals AU aufgetragen. Die 
_Abtastimpulse (schraffiert) bewirken jedesmal eine kleine bleibende Anderung des 


ce ek ee veh 


| 
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Stromes I und damit auch des Magnetfeldes H. Sie ist dem zeitlichen Integral des” 


entsprechenden Spannungsimpulses von U, proportional. 


Sampler I-Regler Spulen 


Feldmessung 


Figur 3 


Aus dem Signalflu8schema (Figur 3) gewinnt man einen Uberblick iiber die 
Stabilitat des 4usseren Regelkreises. Die einzelnen Regelglieder sind durch Recht- 
ecke symbolisiert. Die darin eingezeichneten Kurven bedeuten deren Ubergangs- 
funktion, das heisst ihre Antwort auf die Schrittfunktion. Man erkennt daraus, 
dass schon bei tiefen Frequenzen Phasenverschiebungen von z und mehr auftreten. 
Der Abtaster tibertragt aber nur Frequenzen, die kleiner sind als die halbe Abtast- 
frequenz. Wahlt man diese klein genug, bewirkt das auf einfache Weise ein stabiles 
Verhalten dieses Regelkreises. Die eingangs aufgestellten Forderungen hinsichtlich 
Stabilitat wurden damit erfiillt. 

Wir moéchten Herrn Prof. O. HuBER fiir das Interesse, das er dieser Arbeit ent- 


gegengebracht hat, danken sowie dem Schweizerischen Nationalfonds fiir den 
finanziellen Beitrag. 


LITERATURVERZEICHNIS 


[1] H. Witp und O. Huser, Helv. Phys. Acta 30, 3 (1957). 
[2] O. Huper, L. SCHELLENBERG und H. Wip, Helv. Phys. Acta 33, 535 (1960). 


Etude d’un systéme a 4 bobines pour la production d’un champ 
magnétique homogéne. Par P. Cornaz, Lausanne) 


1. Introduction 


Un systéme a 4 bobines permet d’obtenir des champs trés homogénes avec des 
dimensions réduites et par conséquent une faible dissipation électrique. Par rapport 
au solénoide corrigé on profite encore de l’avantage de l’accés latéral. 

Nous avons calculé et construit la bobine représentée schématiquement sur la 
figure 1. Elle fait partie d’un spectrograph de résonance paramagnétique a 1000MHz. 
Elle permet d’obtenir un champ de 350 Gauss avec un courant de 10,5 A et une 


| 


| 
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puissance de 90 W (refroidissement par ventilation). L’homogénéité a été mesurée 
avec une sonde Hatr 4 l'aide d’un montage d’opposition; elle est meilleure que 
AH/H = 3-10-® sur une distance de 22 mm sur l’axe‘). 


2. Détails constructifs 


Le corps de bobine est en anticorodal. Les bobines en sont isolées par du mica. 
Chaque bobine est constituée de deux flasques bobinés en ruban de cuivre, 10- 0,3 mm 
avec isolation en prespan 0,1 mm et de la laque. Les deux flasques sont reliés au 
centre par une spire épaisse commune; ainsi l’entrée et la sortie du courant sont a 
la périphérie. Le bobinage en ruban donne des bobines solides et précises. D’autre 
part la conductibilité thermique axiale du bobinage est celle du cuivre massif de 
sorte que le refroidissement se réalise trés facilement par les flancs du bobinage. La 
température intérieure des bobines n’est guére plus élevée que celle des flancs. En 
plongeant cette bobine dans un liquide agité on pourrait réaliser des champs de 
Yordre de 3000 Gauss avec une puissance de 6,5 kW et l’homogénéité indiquée plus 
haut. Pour augmenter le refroidissement on pourrait remplacer l’isolant. entre les 
deux flasques de chaque bobine par de |’aluminium oxydé avec une circulation 
d’eau a la périphérie. La construction doit prévoir un déplacement axial de + 2mm 
sur les bobines centrales pour ajuster expérimentalement l’homogénéité a sa valeur 
optimum. 


|. 100 


io 

2 

a 10; 70.\ | 

rt: Figure 1 

; Plan schématique de la bobine. 

— i Pew c Nome 
- 4) Nous remercions M. A, CHATELAIN et M. J. Burret de ce laboratoire de la réalisation soignée 


de ces mesures. 
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3. Calcul 
Une méthode de calcul de tels bobinages épais sera publiée ultérieurement;; 
nous nous bornons ici A un résumé succinct des pas a suivre (figure 1): | 


a) Le champ en P(&) d’une spire de courant I, de rayon 1, située a l’abscisse +! 
peut étre mis sous la forme d’une série ; 


mf Sx (Ef 


9 n=0 


ou 7, est un rayon moyen de référence et Y,, sont les coefficients; ils sont fonction: 
de ~ et 7/7: 


Y 
yay, “2 

b) On considére d’abord les bobines comme infiniment minces. Avec %, 3, 
7,/%x, N,/N, (N = nombre de spires) on a quatre paramétres a varier, ce qui perme 
d’annuler Y, et Y, 4 l’aide des courbes Y, (*, 7/7), tout en satisfaisant a quelque: 
exigences constructives. Les coefficients impairs sont nuls par symétrie. 

-Les données ainsi trouvées permettent de faire un avant-projet des bobin 
complétes. On garde 7, et x, comme variables (figure 1). 


c) Par intégration numérique ou graphique on détermine la contribution 
Y, et Y, des bobines entiéres pour quelques valeurs de x. On trace les courbes 
Y,(¥) et Y,(¥) et on détermine les valeurs x, et ¥2 tel que Yoiot = Yatot = 0. Vu le 
imprécisions inévitables dans le bobinage il est prudent de prévoir une certaine 
latitude de réglage des valeurs ¥, et 7%. 


Performances d’un amplificateur différentiel 4 transistors. Par E. Vit- 
toz et R. DEssouLavy, Lausanne‘). 
1. Introduction 
But et définitions (figure 1a) 


Tensions d’entrée: 


U, jen UR SET ih ie ede a. 
’ 2 ’ 

ou 
UU. = CLE, Cee ae SS 


E doit étre amplifié, U doit étre rejeté. 
Le signal de sortie peut étre l’une des 3 tensions: 


U, 2A(E+", 0), (1) 
U, QA (— Eee, U), (2? 
Uy, =U —-U SAGE PY 7) Ue aA Geo, Ue (3) 


5) Ecole Polytechnique de Il’ Université, 
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On définit ainsi 3 rapports de réjection 71, 72, 74, et un facteur d’amplification 


| différentielle 


AD U; D —U, a) U; — 0; 
Ss at. er ne BD Ie 
On veut: a) A> 1 
b) 7, et 7. <1 (sortie asymétrique) 
ou (et) 74. <1 (sortie symétrique) 


(pou Ui 0)" 


Figure la 
Schéma de l’amplificateur différentiel. 


Effet de Ry 
On peut montrer que 
Ren 1 
rg pitt pa (=2* __ ) 
cite 2, aay eee 2 


pour deux transistors 7, et 7; parfaitement identiques. Nous avons mesuré cette 


courbe (figure 2) en compensant les dissymétries entre T, et T, (c.f. 2). 


ee eee ba tare hn aes 


05-10°* 


0 
5 70 R, (Ma 


-10"* 
Figure 2 
Effet de Ro. 
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; 
Pour obtenir 7, et 7, < 1, il faut que la valeur de R, soit voisine de | 


Iiare n A 
hove — heey 


2. Montage pratique (figure 1a) 
Equilibrage et compensations 


Pratiquement les 2 transistors T, et T, ne sont jamais identiques. 

a) Pour Vgr1 = Veee, Ic: # Ice (tensions et courants continus). 
Le potentiométre P, permet d’équilibrer les courants continus au repos dans les 
2 branches de l’amplificateur. 

b) Pour U,, = U,, = U, Ig, F Ln et U2 FO 

(tensions et courants alternatifs BF.) 

On compense cette dissymétrie en répartissant judicieusement dans les charges R, 
un courant proportionnel 4 U, ceci au moyen des éléments P, et R. 

c) En HF, la différence de phase entre J,, et [,, est compensée en ajustant le 
condensateur différentiel Cj. 


Réalisation de R, 


On a vu (figure 2) qu’il était intéressant d’avoir une valeur de R, élevée. Si 
cela est impossible avec les courants de repos choisis et la tension positive V, a 
disposition, on peut recourir au montage donné a la figure 1b. 


Figure 1b 


Source de courant a résistance interne élevée. 


On obtient Ry > 1/hg., pour R; > hy, 


3. Mesures 


Mesure de A 
On applique une tension 2 E entre les 2 entrées. 
AD U, a) —U, oa Uy. 
yg eS ae a ai 3 


Mesures de 11, rg et 115 


On applique une tension U simultanément aux 2 entrées (bases reliées ensemble), 
Les équations (1), (2) et (3) permettent de déterminer RAY IE 
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Valeurs numériques 
Les valeurs numériques utilisées pour nos mesures sont les suivantes: 


R, =10KQ, 


P, = 150 Q (une valeur plus faible peut convenir et présente l’avantage d’une aug- 
mentation de A) , 
P, = 250 KQ, 


R = 0,5 a 2 MQ (suivant le type des transistors utilisés) , 
. C, = 100 pF différentiel. 
Point de fonctionnement des transistors au repos et a 20°C 
—Io 2 250 pA, —Vcor 25V. 


Types de transistors: OC 71 (Germanium) 
BCZ 11 (Silicium) 


Nous avons effectué chaque série de mesures avec les 6 combinaisons de 4 tran- 


sistors de méme type. 
Facteur d’amplification différentielle: 


A & 50 (varie peu d’une paire a l’autre). 


Rapports de véjection 


| Nous avons fait subir 4 chaque paire de transistors un cycle de températures 
| 20°, 40°, 60° et retour A 20°C en relevant chaque fois les valeurs de 71, 7, et "2. 

i La figure 3 montre un exemple des résultats obtenus. On remarque une dispa- 
| rition de l’effet des compensations aprés un cycle de températures (transformations 


Drie ble BF O20 20° = 40° = 60° 20°C 


Ih max 360 C 


hac’ as 


15-10" 


aA. 4 


Figure 3 


Variation des rapports de rejection au cours d’un cycle de températures pour 6 paires de transisto 
Diagramme de gauche: sans compensation. Diagramme de droite: avec compensation BF a 20°C, 
T, et Ty: OC 71 (Ge), Ro= 2MQ, fmesure = 150 Hz. Courants continus équilibrés 4 20°C. 


rs. 


Cee ee i ee Ae eee en 


at ete tek De 908 
’ 
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irréversibles des transistors). Ceci n’apparait pas a ces températures avec des 
transistors au Si. . 

Nous donnons (tableau 1) les valeurs de (| 7, | 0u | 72 |) max €t | %12 | max Mesurées. 

Du point de vue des rapports de réjection, on ne constate pas de différences — 
appréciables entre les transistors au Ge et au Si. I 

Par contre la dérive équivalente max a l’entrée pour T variant de 20° a 60°C _ 
qui est de 10 4 20 mV (selon la dissymétrie des courants inverses) avec des tran-_ 
sistors au Ge n’est que de 1 mV environ avec des transistors au Si. On aura done 
avantage a utiliser ces derniers dans les amplificateurs travaillant a faible niveau 
et soumis a de fortes variations de température. 


Tableau 1: Rapports de réjection max mesurés. 


| Transistors OC 71 (Ge) | Transistors BCZ 11 (Si) 
ifcaregiee 150 Hz 5 Iz 150 Hz 5 KHz 
Ry 12KOQ| ~2MQ | Z~21MQ | ~2MQ | Z,~1MQ 


| 
Compens. BF a 20°C | non | non | oui oui oui 


Compens. HF a 20°C | non | non | non | non | oul 


Re) eens Peay eT oe 
10°| 7% |max | 40°C pee te te ta 
ou 8 Seer 
Fa be Mtb athe eas nein CA he 
20°C xk Pa ee hee 
20°C | 73 me le Spat 
aed ong ea cc 
eet 60°C, 1 96 18 13 13 14 
20°C |, 72 | It |) > eel oe pole) oF ee 


4. Conclusion 


Les mesures effectuées confirment lavantage d’une valeur élevée de R, si on 
veut que | 7, | et | 7, |.soient faibles. 

Un ajustage du dispositif de compensation permet d’obtenir facilement 
| “12 | < 10-5 a température constante. 


Avec des variations de température de 20° A 60°C on peut encore compter sur 
une valeur de |7,.| < 1,5- 10-4. 


Compteur d’impulsions bidirectionnel A tores magnétiques. Par R. L 
Favre, Lausanne‘). 
Indépendamment de Vapplication des ferrites aux mémoires matricielles dont 
il serait vain de discuter l'intérét, les compteurs d’impulsions a tores magnétiques 
sont loin, semble-t-il, d’avoir donné la limite de leurs possibilités. 


6) Laboratoire de Recherches Nucléaires. 
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Les circuits actuellement parvenus a notre connaissance [1]7) reposent, avec 
quelques variantes, sur le schéma de la figure 1. Un seul des tores de ferrite est 


Figure 1 


polarisé dans le sens inverse 4 celui correspondant au courant du transistor, lequel 
est bloqué a l’état de repos. Supposons qu’il en soit ainsi du tore 0: Un bref déblo- 
cage du transistor rectifie la polarisation de ce tore, d’ot charge du condensateur 
tampon Co, a travers la diode de blocage do. La décharge différée de Co dans le 
bobinage du tore 1, provoque l’inversion magnétique de ce dernier jusqu’a l’impul- 
sion d’entrée suivante, et ainsi de suite. 

L’intérét du comptage bidirectionnel et du relevé de la mémoire, dans un 
probléme d’automation en cours d’étude, nous a conduit au développement d’un 
circuit apte a aiguiller, de maniére sélective, la décharge du condensateur tampon 
dans le tore prédéterminé par l’opération requise. La décharge.d’un condensateur 

' tampon quelconque, par exemple C/ (figure 2) est asservie au déblocage de l’une 
| ou l’autre des voies de commutation contrélées par les transistors Td et Tr. Les 
| 


Figure 2 


7) Le chiffre entre crochets renvoit a la Bibliographie, page 466. 
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' 
impulsions d’entrée sont requises de double polarité, la premiére alternance dé 
bloquant le transistor To et la seconde, l’un ou l’autre des transistors de commu- 
tation Td ou Tr selon que l’impulsion est appliquée a l’entrée directe Ed ou rétro- 
grade Ev. Dans le premier cas, c’est le tore suivant qui s’inverse et dans le second, 
le tore précédent. ‘| 

Le procédé est sans autre applicable 4 une décade (10 tores en série). Il est 
possible d’augmenter le taux d’information spécifique en constituant deux cycles 
par série. C’est ainsi qu’une série de 5 tores, divisée en un cycle de 2 et l’autre de 
3 tores, permet 6 états distincts au lieu de 5. 

La généralité de la méthode d’aiguillage sélectif permet son application au 
relevé de la mémoire par affichage lumineux, impression automatique, etc. I] suffit 
d’introduire une voie de commutation supplémentaire, réinversant toujours le 
méme tore pour qu’apparaisse sur le condensateur tampon de ce dernier une 
succession d’impulsions propres 4 permettre son identification. Ces impulsions sont 
entre autre applicables au déblocage d’un tube d’affichage du type DM 160. 

Les contingences d’une communication qui se doit courte, nous contraignent a 
renvoyer le complément de nos informations a un prochain article. 

Il nous reste & exprimer nos sentiments de gratitude 4 Monsieur le Professeur 
HaEnny, a notre technicien Monsieur BAIcHE et a la Commission Suisse pour la 
Science Atomique, en reconnaissance de leur contribution respective a ce travail. 


BIBLIOGRAPHIE 
[1] J. C. Barton, Simple core scaling circuits, Nuclear Inst. 5, 332-334 (1959). 


Contribution des transistors aux performances de circuits 4 tubes 
electroniques. Par R. L. Favre, Lausanne’). 

Les quelques circuits élémentaires proposés ci-dessous s’ajoutent, dans le méme 
esprit, aux communications [1] et [2]9). 

a) La figure 1 schématise un limiteur de courant pour stabilisateur électronique 
de tension, artifice dont nous apprécions journellement les services dans le dévelop- 


Figure 1 


Circuit limiteur de courant. 


8) Laboratoire de Recherches Nucléaires. 
®) Les chiffres entre crochets renvoient A la Bibliographie, page 468. 
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_pement et la mise au point de circuits. Par son action sur la tension d’écran du 
tube rhéostat EL 86, la résistance anode-base de polarisation complémentaire d’un 
transistor tampon, permet de limiter le courant de court-circuit 4 une valeur égale, 
-voir inférieure, au courant utile maximum. Dans le premier cas, la caractéristique 
horizontale tension-courant se termine par une chute verticale. Le courant utile 
maximum est alors indépendant de la tension stabilisée, d’ot la possibilité de calibrer 
le potentiométre R, par exemple entre 0 et 100 mA, a + 10%, pour une tension 
stabilisée variable de 0 4 1000 V, tel que nous |’avons expérimenté. 

b) D’un emploi généralisé par celui des photomultiplicateurs, les hautes tensions 
stabilisées d’un débit de plusieurs mA s’accommodent mal de la résistance interne 
élevée des oscillateurs 4 haute fréquence, par ailleurs trés avantageux pour la caté- 

_ gorie de tensions intéressée. La figure 2 illustre le principe appliqué avec succés pour 


300 V 


Figure 2 


Circuit compenseur de débit. 


-annuler la résistance interne d’une tension redressée dans ces conditions. Le courant 
de retour du bobinage HT est partiellement dérivé vers la base d’un transistor 
couplé de maniére a ce que la polarisation de grille du tube oscillateur EL 86 soit 
d’autant moins négative que le courant absorbé est plus élevé. Il est aisé d’annuler 
sensiblement la résistance interne pour une tension de 2000 V, par exemple, tandis 
que la compensation est encore trés efficace aux valeurs extrémes de 500 et 2500 V. 

c) La figure 3 est relative 4 un procédé de mise en forme des impulsions d’un 

' photomultiplicateur, desquelles il était requis une amplitude maximum et une 


220 
ASZ21 
1K 
4uH 
Ea 
Figure 3 


Circuit de mise en forme. 
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: . . x . | 
durée trés réduite, en vue de leur application 4 un processus de commutation aussi 


rapide que possible. Le front négatif de 2 a 30 V est généré sous haute impédance > 
résistive, tandis que le front de retour fait intervenir la brusque conduction dun 
transistor «de rappel» qui rétablit, en quelques m s, la tension de repos. Les per-_ 
formances peuvent étre améliorées aux dépens de la simplicité. | 

Nous remercions Monsieur le Professeur HAENNyY, nos techniciens Messieurs — 
Jomint, HotzEeR et QUINCHE, ainsi que la Commission Suisse pour la Science — 
Atomique, de leur précieuse contribution respective a ce travail. 
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[2} R. L. Favre, Eléments de compensation d’un stabilisateur de tension, Helv. Phys. 
Acta, 33, 598-599 (1960). 


Ein- und Ausdiffusion von Dotierungssubstanzen in Silizium. Von 
H. Kunz, Baden??). 

Zur Erreichung bestimmter Konzentrationsprofile mit kleinen Randkonzentra- 
tionen bei der Eindiffusion von Fremdsubstanzen in Silizium wurde das folgende 
Diffusionsverfahren gewahlt: 

1. Eindiffusion der Fremdsubstanz in Siliziumscheiben aus der Gasphase. 

2. Nachtemperung der vordiffundierten Siliziumscheiben ohne weitere Zufiih- 
rung von Dotierungsmaterial. 

Der erste Teil des Verfahrens, Eindiffusion der Dotierungssubstanz aus der 
Gasphase, fiihrt im eindimensionalen Fall auf eine Konzentrationsverteilung der 


Form: 
te 
C(x pd) Go.) 1-8 | Ih, 
: ‘ ( 2 VDt, ) 


C, = Randkonzentration ; 2 ae = Fehlerintegral. 


2 Di 


Im zweiten Teil des Verfahrens, Nachtemperung ohne weitere Zufiihrung von 
Dotierungsmaterial, findet einerseits eine weitere Eindiffusion, anderseits aber auch 
eine Ausdiffusion der Fremdsubstanz aus dem Silizium statt. 

Pur dieses Verfahren der Diffusion in 2 Schritten wurden die Konzentrations- 
profile von WarInG [1]") unter folgender, vereinfachender Annahme mit einem 
Analog Computer bestimmt: Die Randkonzentration der Si-Probe soll wahrend der 
Nachtemperung dauernd auf 0 oder konstant gehalten werden. 

_ Fur diese Falle kann aber die Diffusionsgleichung gelost werden. Zur Erreichung 
einer wahrend der Nachtemperung konstant gehaltenen Randkonzentration machen 
wir fur die Anfangsbedingungen folgende Annahmen: 


Case, l be Neral fiir 0 <x <0o, 
2/Dt, 


10) Physiklaboratorium Brown Boveri & Cie. 
le ap ee F 3 . : : 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 471. 
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C=C|—(1- A) 9 i ) fir 0>%>— 0, 


2) Di 


-wobei f/2 C, die konstante Randkonzentration bedeutet. 
Als Lésung der Diffusionsgleichung erhalten wir dann eine Konzentrationsver- 
teilung 
Canc: | -!\o gy Sow ac spect |e aa 
2 : 2 / a a 2 VDt 


a+ 1 Vides 
‘t, = Zeit der Nachtemperung , 

Fats Zeit der Nachtemperung, 

t, Zeit der Diffusion aus der Gasphase ’ 


M, +t,=1t. 


gy = 


Machen wir nun noch die Annahme f/2 = 0, so erhalten wir auch die Konzen- 
trationsverteilung fiir den Spezialfall, dass die Randkonzentration wahrend der 


— | Dt 
a+l 


i ———S= SS a a sae @ : aad @ — 
; Verlauf der Funktion C 97 2 2/ ‘3 . 2\/Dt 


fiir verschiedene Werte von uw iM a = konstant = 10. 


Paw ae 


1 
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Nachtemperung dauernd auf 0 gehalten wird. 


Clr.) = C, |e x 2 ae (2) 
a yDi fo | 


t 


Die Funktion C/C, fiir « = 10 und verschiedene Werte von #/2 zeigt Figur 1. Aus 

dieser Darstellung geht hervor, dass die Gleichung (2) fiir alle Falle mit grosser — 

Ausdiffusionsrate [8/2 = klein] eine gute Naherung ist. 
Die Funktion C/C, fiir 8/2 = 0 und verschiedene Werte von « zeigt Figur 2. 


; 


Figur 2 


Verlauf der Funktion ee =D Z 


(ape) 


) fiir verschiedene Werte von «, 


Von uns durchgefiihrte Versuche zur Eindiffusion von Ga in Siliziumscheiben 
haben gezeigt, dass der gemessene Konzentrationsverlauf sehr gut mit dem nach 
Gleichung (2) berechneten Verlauf iibereinstimmt. Figur 3 zeigt die gemessenen 
Werte der Konzentrationsverteilung und den nach Gleichung (2) berechneten 
Konzentrationsverlauf fiir Ga in Si mit « = 95. Fiir die gerechnete Kurve wurde 
die von FULLER und DitzENBERGER [2] angegebene Diffusionskonstante verwendet. 
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Figur 3 


Gemessener und nach Formel (2) berechneter Konzentrationsverlauf fiir Ga in Si « = 95 
berechneter Konzentrationsverlauf. 


LITERATURVERZEICHNIS 


1] W. Warrne, J. Electrochem. Soc. 705, 12, 702 (1958). 
] C. S. Futrer, J. A. DitzenBerceEr, J. Appl. Phys. 27, 544 (1956). 


: Couplage de deux milieux hétérogénes multiplicateurs de neutrons. 
Par Cu. Manprin, P. Gavin et B. Vittoz"). 
A paraitre prochainement dans ZAMP. 


Kompressibilitatseinfliisse bei Gaslagern. Von V. STINGELIN, Ziirich"*). 

Die Differentialgleichung der hydrodynamisch geschmierten Gaslager lasst sich 
bekanntlich, selbst fiir das unendlich lange Lager (ebenes Problem), nicht geschlos- 
sen integrieren. Man war daher lange Zeit auf Naherungslésungen mit Hilfe von 
Potenzreihenentwicklungen angewiesen [1], [2]*), deren Giiltigkeitsbereich natur- 
gemass beschrankt war. Um diesen zu erweitern, haben wir die Differentialgleichung 
fiir das unendlich lange Lager auf der Analogie-Integrieranlage [A-58, der Firma 
Contraves berechnen lassen. 

In Figur 1 sind die Druckverteilungen fiir ein festes « in Abhangigkeit von c 
‘aufgetragen. c ist die Ahnlichkeitszahl der Differentialgleichung und direkt ein 
Mass fiir die Starke der Kompressibilitatseinfliisse 1°). e = relative Exzentrizitat. Die 
‘Druckverteilungen sind offensichtlich ausserst empfindlich auf Variationen von ¢. 
Obwohl wir in der Praxis das unendlich lange Lager nicht exakt realisieren 
‘kénnen, durften wir doch bei den Messungen am endlich breiten Lager ein qualitativ 
12) Lab. de Génie Atomique, EPUL, Lausanne. 


) 

13) Institut fir Aerodynamik, BH. 

14) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, 
) 


c = 0: inkompressibel (SOMMERFELD); ¢ > 0: kompressibel. 


Seite 473. 
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ahnliches Verhalten der Druckverteilungen erwarten. Was wir bis jetzt veel 
konnten, waren Rechnungen!*) und Messungen fiir ein L/D = 1 (L = Lagerlang' 
D = Lagerdurchmesser), wobei Ciyeor = 0 UNA Copp = 0,079 betrug. Siehe Figur 


oa 


‘| 
a 


: 


a 


i 


0.5 


‘ - tm ot 
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Figur 1 
Berechnete Druckverteilungen des unendlich langen Gaslagers. 


Die Ubereinstimmung Messung—Rechnung war fiir ¢ = 0,5 recht gut, hingegen 
machten sich bei e = 0,7 die erwarteten Kompressibilitatseinfliisse bereits bemerk- 
bar. Bei den relativ kleinen Oberflachengeschwindigkeiten durften wir diese Ab- 
weichungen von der inkompressiblen Lésung mit Sicherheit auf die bereits vari- 
ierende Dichte zuriickfiihren und nicht auf Sekundireffekte innerer Temperatur- 
schwankungen. 

Die nachtragliche energetische Betrachtung hat uns diese Vermutung bestatigt 
und zudem ein Mittel gegeben, die Grenzen der Theorie, welche isotherme Zustands- 
anderungen im Gas voraussetzt, von Fall zu Fall abzukliren. 

Aus dem Energiesatz in dimensionsloser Darstellung ziehen wir folgende 
Schlussfolgerungen: 

Solange das relative Lagerspiel y eine oder mehrere Potenzen von Zehn kleiner 
ist als der reziproke Wert der Pe-Zahl (Pe = Res Pr, Res = Up 6/y), tiberwiegt die 
Warmeleitung senkrecht zur Stromungsrichtung diejenige in Richtung der Stré- 

16) Die Berechnungen des endlich breiten Lagers, c = 0, wurden auf der ERMETH ausgefiihrt. 
Die Rechnungen fiir c > 0 sind zur Zeit noch im Gange. 
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ung und den Warmetransport. Unter diesen Voraussetzungen wird in erster 
Naherung die Temperaturverteilung eine Funktion vom Druckgradienten und der 
AShe z allein. Demzufolge lasst sich der Energiesatz integrieren, und wir erhalten 


Schnitt = 140° 


Mittelebene 


* 
Theorie inkompr. (C = 0), = , pr= ee a, OSs 


0 
—5=—0,7, C=0,0788, Urp=6,3 m/sec, po= 1 bar. 
Or e=0,5, C=0,0785, Ug=6,3 m/sec, py = 1 bar. 


Figur 2 
Druckmessungen am endlich breiten Gaslager, L/D = 1. 


die Temperaturverteilung als Funktion von z und dem Druckgradienten Pp als 
Parameter. ait 

- Die Auswertung dieser Berechnung ist in Figur 3 wiedergegeben, und zusammen 
mit dem kleinen numerischen Beispiel, welches die Bezugstemperatur mw UR/A fiir 
einen typischen Lagerfall angibt, erhalten wir einen Begriff von den zu erwartenden 
Temperaturveranderungen. Die Groéssenordnung betragt hier rund 2 Grad, womit 
die Brauchbarkeit der isothermen Theorie fiir dieses Beispiel nachgewiesen ware. 


a LITERATURVERZEICHNIS 


[1] J. S. Ausman, The Fluid Dynamic Theory of Gas-Lubricated Bearings. Trans- 
action of the ASME 79, 1218-1224 (August 1957). 
[2] Y. Katto and N. Sopa, Theory of Lubrication by Compressible Fluid, with 
_ Special Reference to Air Bearing. Proc. of the 2nd Japan Nat. Congr. for App. 
Mech. 267-270 (1952). 
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Beispiel: Luft bei300°K, Urp=50m/sec, 6=30u, 
sec? 
cm2 


+ = 6,93- 10-9] 


; Grd.]_, f UL = 173°, Pew. 


| = 
Figur 3 

Geschwindigkeitsprofile und Temperaturprofile im Gaslager in Abhangigkeit vom Druckgradienten 
Pg. Voraussetzung ist: p< 1/Pe. Randbedingungen: Twana = "weltle = To = konst. 


Eigenschaften und AnwendungsmO6S$lichkeiten des Omegatrons. Von 
R. H. Bieri und W. RutscuH!’). 

Erscheint demnachst als ausfiihrliche Arbeit, voraussichtlich in Review of Sci- 
entific Instruments. 


Uber einen Einfluss des AgJ auf den atmospharischen Vereisungs- 
prozess. Von RoLanp List, Weissfluhjoch-Davos }8). 

Grdéssere atmospharische Eispartikeln wie Graupeln und Hagelkérner wachsen 
in der Atmosphare zur Hauptsache durch Akkreszenz kleinerer Wolkenpartikeln. 
Handelt es sich hierbei um fliissige Wassertropfen, so bleiben dieselben beim Auf- 
prall auf ein grésseres Hindernis infolge Benetzung haften. Fallt jedoch ein Hagel- 
korn in einer Wolke, die nur aus Eisteilchen — seien es gefrorene Trépfchen oder 
Schneekristalle — besteht, so ist seine Wachstumsrate klein, da ein grosser Teil der 


1”) Physikalisches Institut der Universitat Bern. 
18) Eidg. Institut fiir Schnee- und Lawinenforschung. 
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auftreffenden Partikeln wegprallt und nur wenige haften bleiben. Eine gewisse 
Klebekraft zwischen sich beriihrenden Eisteilchen besteht; sie ist jedoch gering und 
immt mit sinkender Temperatur ab [1]1). 


0 -10 -20 -80 


Figur 1 
Spezifischer Massenzuwachs dm/dt einer vereisenden Kugel in Funktion der Lufttemperatur t;. 
Relative Luftgeschwindigkeit v = 12 m/s, Durchmesser der vereisenden Partikeln d = 2 cm, 
Wassergehalt der Luft (in Form von Tropfen); wy = 10,5 g/m’. 
Kurve I: Massenzuwachs ohne Verwendung von Ag J-Partikeln. 
Kurve II: Massenzuwachs unter Verwendung von Ag J-Partikeln. 


Wie Figur 1 zeigt, wird diese Vorstellung durch Experimente im Hagelversuchs- 
kanal [2] bestatigt. Der Massenzuwachs pro Zeiteinheit wurde bei konstanten 
Werten der relativen Luftgeschwindigkeit vereisendes Objekt — Partikelwolke, des 
freien Wassergehaltes, bzw. Eisgehaltes der Luft sowie des Durchmessers des ver- 
-eisenden kugeligen Objektes in Funktion der Lufttemperatur gemessen. Kurve I 
zeigt dabei das spezifische Wachstum einer Kugel, das sich innerhalb einer Tropfen- 
_wolke einstellt. Es bleibt konstant bis zu Werten von —33°C und fallt dann rasch 
ab, um bei —40°C einen Betrag anzunehmen, der 5% des urspriinglichen Wertes 
_betragt. Dieser Abfall riihrt davon her, dass bei diesen Temperaturen ein mehr oder 

weniger spontanes Gefrieren der Wassertropfen eintritt, sei es auf Grund der maxi- 
-malen Unterkiihlbarkeit des verwendeten Tropfengréssen-Spektrums (Durchmesser 
des mittleren Volumens 50 y) oder aus apparativen Griinden. 

~ Wiederholen wir die gleiche Experimentreihe bei gleichen Versuchsbedingungen, 
jedoch unter Beimengung von AgJ-Partikeln zum Tropfennebel vor der Mefstelle, 
wo der Vereisungsversuch stattfindet, so ergeben sich Massenzuwachsraten gemass 
-Kurve II. Obwohl die als Eisbildungskerne wirkenden Ag J-Teilchen einen Schwel- 
lenwert von —4°C aufweisen und daher bereits unterhalb dieser Temperatur eis- 
pildend wirken, so lisst sich erst bei —10°C ein Einfluss feststellen. Dies ist nicht 
verwunderlich, da das Verhaltnis zwischen den wirksamen Kernen zu den total vor- 
handenen exponentiell mit fallender Temperatur zunimmt (gemass [3] und [4] von 
4° bis —10°C um einen Faktor 10? bis 10%). Das heisst, dass bei Temperaturen 
hoher als — 10°C ein wesentlicher Teil der Wolkentrépfchen (ca. 160/cm* bei @ =50u 
und einem totalen freien Wassergehalt von w = 10,5 g/m’) fliissig bleibt und dadurch 
“eine nasse Oberflachenschicht auf der vereisenden Kugel verursacht, in der auch die 
‘auftreffenden Eispartikeln kleben bleiben. Der Abfall zwischen —10° und —15°C 


19) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 476. 
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ist also dem Gefrieren aller Wolkentropfen zuzuschreiben, wahrenddem der f 
gende, leicht sinkende Ast der Kurve II bei noch tieferen Temperaturen auf d 
oben erwahnten Abfall der Haftung zwischen reinen Eispartikeln zuriickzufiihren 


ist. 


Die hier beschriebenen Experimente wurden unter Verwendung sehr hoher Kern- 
konzentrationen durchgefiihrt (Gréssenordnung 10°/cm* der bei —10°C wirksamet 
AgJ-Partikeln). Die Zeit, die dem Vereisen der Tropfenwolke zur Verfiigung stand, 
betrug 0,7 s, gegeben durch die Stromungsgeschwindigkeit und den Abstane 
zwischen Einspritzstelle und Mefstelle. 

Zusammenfassend kann also festgestellt werden, dass die Grundidee einer Hage 
bekampfung durch Gefrieren samtlicher unterkiihlter Wolkenpartikeln prinzipie! 
richtig ist. Die angefiihrten Resultate erlauben jedoch noch keine allgemeiner 
Schlussfolgerungen fiir die Praxis, da speziell die Fragen der Kernwirksamkeit i 
Funktion der Zeit, des Kernwachstums durch Diffusion sowie der Kerndosierun 
noch einer Lésung harren. 


. 
LITERATURVERZEICHNIS : 


[1] R. E. Hatrcren und C. L. Hosier, Geoph. Monograph Nr. 5, Amer. Geoph. 
Union, 257—263 (1960). | 
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[3] H. Wercxmann, in H. Byers, Thunderstorm Electricity, University Press, Chi- 
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Uber die Erzeugung eines prazisen 1,5-Mill.-Volt-Elektronenstrahles 
fiir Elektronenmikroskopie. Von E. B. Bas, Ziirich?*). 

In der Durchstrahlungs-Elektronenmikroskopie werden normalerweise Elek 
tronenenergien bis zu 100kV angewendet. Es bestehen verschiedene Beweg: 
griinde, die Elektronenenergie weiter zu erhéhen. Einer davon ist zum Beispiel dic 
MOglichkeit, dickere Objekte zu durchstrahlen. Man kann beispielsweise die leben 
den Bakterien in feuchtem Milieu zwischen zwei durchstrahlbaren Folien vakuum 
dicht einbetten und durchstrahlen, wobei eine hohe Energie auch die Abtétungsrats 
vermindern kann. Die Hauptschwierigkeit einer Elektronenmikroskopie mit seht 
hohen Spannungen besteht in der Erzeugung des hochkonstanten Elektronen- 
strahles. Ein fiir diese Zwecke geeigneter Elektronenbeschleuniger wurde von det 
Firma Haefely in Basel in Zusammenarbeit mit der Abteilung fiir industrielle 
Forschung (AFIF) an der ETH gebaut und steht im Institute d’Electronique (Prof 
Dupovy) in Toulouse fiir elektronenmikroskopische Untersuchungen zur Verfiigung 
Im folgenden wird kurz iiber das Beschleunigungsrohr, welches an der AFIF ent. 
wickelt wurde, berichtet. 

Der Beschleuniger ist fiir 1500 kV und 200 wA ausgelegt. Die Spannung wir 
aut 10~® und der Strom auf 10-? stabilisiert. Die ganze Anlage ist frei aufgestellt 
Die hochstabilisierte Spannung wird von einem symmetrischen Kaskadengleich 
richter erzeugt {1]*4). Um die geforderte hohe Stabilisierung der Spannung z 
gewahrleisten, musste die ganze Anlage und auch das Beschleunigungsrohr seh 
sorgfaltig dimensioniert und ausgefiihrt werden. Die Beschleunigungsstrecke ist it 


20) Abteilung fiir industrielle Forschung, ETH. 
*1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 481. 
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0 Stufen unterteilt. Die Elektronen werden also pro Stufe um 150 kV beschleunigt 
Diese kleine Stufenspannung ermdglicht die Unterdriickung der Korona-Entladung 
usserhalb des Rohres und eine sichere Beherrschung des elektrischen Durch- 
chlages zwischen den Beschleunigungselektroden. Die Stufenhéhe betragt 45 cm. 
igur 1 zeigt schematisch den Aufbau der Beschleunigungsstrecke. Aus einem 


HSp. 
Kaskade 


605 cm ++, 
405 
as 
inn Mish Rs, Bg OO ee 


200 


> Figur 1 
Prinzipschema des Beschleunigers: A Injektionssystem, P Elektronenquelle, Q Bild der Elektronen- 
; quelle. 


- 


EinschuBsystem A werden die Elektroden mit variabler Energie in die Beschleuni- 
gungsstrecke eingeschossen. Der EinschuBstrahl kann in einem beliebigen Punkt P 
in der Nahe der Achse vor der ersten Beschleunigungsstufe fokussiert und um 
diesen Punkt P gekippt werden. Dadurch ist ein genauer Einschuss gewahrleistet. 
In Figur 2 ist das Einschu8system schematisch dargestellt. Es besteht aus einer 
‘Fernfokus-Elektronenkanone und einer magnetischen Kondensorlinse. Die Elek- 
‘tronenkanone fokussiert den Strahl auf eine Crossoverblende von 0,2mm @, die 
Kondensorlinse bildet diese Blende im Punkt P ab, wobei der Offnungswinkel des 
‘Strahles durch eine Aperturblende auf 1,25- 10-2 rad beschrankt wird. Das ganze 


rd 
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Erste Beschleunigungslinse 
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Figur 2 
Prinzip des Injektionssystems. 


System kann in zwei Richtungen x, y verschoben werden, wodurch der Strahl in. 
einem beliebigen Punkt vor der 1. Beschleunigungslinse fokussiert werden kann. 
Ausserdem lasst sich die Elektronenkanone samt der Aperturblende in einer Kugel-. 
kalotte mit der Crossoverblende als Zentrum kippen. Der Strahl bleibt dadurch im: 
Punkt P fokussiert; nur seine Neigung zur Achse wird verandert. Die im Durch- 
messer gross dimensionierten Polschuhe der Kondensorlinse erlauben einen schiefen 
Durchgang des Strahles, ohne dass sich die spharischen Fehler des Linsenfeldes aus- 
wirken. Der Grossteil des von der Kathode emittierten Stromes wird durch die 
beiden Blenden ausgeblendet; zum Beispiel bei einem EinschuBstrom von 100uA 
bendtigt man einen Emissionsstrom von 500 HA. Die Beschleunigungsspannung in 
der Elektronenkanone U, kann zwischen 20 und 60 kV variiert werden. 

In Figur 3 ist die konstruktive Gestaltung des Einschu8systems und der Be-, 
schleunigungsstrecke ersichtlich. Das Einschu8system kann von der Beschleuni- 
gungsstrecke mittels eines Schieberventils abgetrennt werden, so dass fiir den 
Kathodenwechsel nicht das ganze Rohr beliiftet werden muss. Bei der Konstruktion 
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Figur 3 
Konstruktiver Aufbau des Injektionssystems und der Beschleunigungslinsen: 1 Kathodenpatrone 
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2, Fernfokus-Elektronenkanone, 2 Justierung des Elektronenstrahles auf die 
des Injektionsstrahles, 5 Servomotoren fiir «,$-Kippung, 6 Kon- 
densorlinse, 7 x,y-Verschiebung des Injektionsstrahles, 8 Servomotoren fiir x,y-Verschiebung, 
§ Lufteinlass- und Vorevakuierungsventil, 10 Absperr-Schieberventil, 11 Doppel-Messkafig fir 
Strahlstrommessung und -justierung, 12 Servomotor fiir Messkafigantrieb, 13 Beschleunigungslinse, 
14 Abschirmelektroden, 15 Porzellanisolator, 16 Justierbare Halterung der Linsenelektrode (drei 
am Umfang), 17 Stufenblende, 18 Korona-Abschirmung, 19 Halterungsnasen (drei am Umfang), 
° 20 Anschlussflansch fiir eine Vakuum-Messzelle. 


mit Bolzenkathode, 
Crossoverblende, 4 «,$-Kippung 
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leunigungsstrecke waren drei Punkte von Bedeutung: prazise Justieru 
ie SS tee Aan REG Vee Sicherheit in der Spannungsfestigkeit und gute 
Abschirmung des magnetischen Erdfeldes. Die Beschleunigungselektroden a 
verschiebbar und kippbar angeordnet. Sie werden bei der Montage durch einen 
mit hochempfindlichen Libellen versehenen Zentrierdorn sehr genau aufeinander 
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Figur 4 


Profil des Strahles auf der Beschleunigungsstrecke und die Ablenkung des Hauptstrahles im 
magnetischen Erdfeld. 


ausgerichtet. Um das Erdfeld gut abzuschirmen, sind sowohl Beschleunigungs- 
elektroden wie die um diese angeordneten Abschirmelektroden aus hochpermeablem 
Material (Mumetall)hergestellt. In Figur 4 ist das gerechnete Strahlprofil [2] und die 
Ablenkung des Zentralstrahles im Erdfeld beim Fehlen einer Abschirmung graphisch 
dargestellt, und zwar fiir eine Beschleunigungsspannung von 1500 kV. Ware eine 
magnetische Abschirmung nicht vorhanden, so wiirde demnach der Strahl in det 
10. Stufe um etwa 4cm und am Eingang des Mikroskopes um 7,3 cm durch das 
Erdfeld abgelenkt sein, und diese Ablenkung wiirde stark durch die Veranderung 
der Beschleunigungsspannung variieren. Dank der guten Abschirmung kann nun 
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die Spannung von 200 kV aufwarts variiert werden ohne eine beobachtbare Strahl- 
ablenkung am Ausgang. Dadurch dass der Strahl genau auf der Achse gefiihrt wird, 
k6nnen in den Beschleunigungselektroden dem Strahlprofilangepasste Stufenblenden 
(17 in Figur 3) vorgesehen werden. Die Aufgabe dieser Blenden ist die Versperrung 
des Weges fiir die auf der Beschleunigungsstrecke durch Zusammenstoss der Strahl- 
elektronen mit den Restgasmolekiilen entstandenen Ionen und Elektronen. Dadurch 
ist eine Verbesserung der Spannungsfestigkeit des Rohres und auch die Verminderung 
‘der Réntgenstrahlung zu erwarten, welch letztere auf die Ausbildung von Korona- 
Entladung in der ganzen Anlage von Bedeutung sein kann. 

Der berechnete Strahldurchmesser am Ausgang des Beschleunigungsrohres 
stimmt sehr gut mit den gemessenen Werten iiberein; man erhalt zum Beispiel bei 
einer Einschu8spannung von 40 kV einer Beschleunigungsspannung von 1000 kV 
fiir den Strahldurchmesser am Ausgang 0,78 mm gerechnet und etwa 0,7 mm 
gemessen. : 


| LITERATURVERZEICHNIS 


‘{1) H. Apter, R. MINKNER, G. REINHOLD und J. Seitz, Delft Conference on Elec- 
| tvon Microscopy, 1960 (im Druck). 
[2] E. B. Bas, L. Preuss und W. SCHNEIDER, Z. angew. Math. Phys. 70, 533 (1959). 


<f Uber die Ausmessung von progressiv-periodischen Potentialfeldern 
mit einem Widerstandsnetzwerk. Von L. Preuss, Ziirich™). 
Betrachten wir das Feld eines elektrostatischen Beschleunigers mit einer hohen 
Anzahl geometrisch gleicher Elektroden (Figur 1). Dieses Feld wollen wir im Grenz- 
fall, wo diese Zahl gegen unendlich strebt, und also gilt: 


g(r, 2)=9(2+1I)—K, (1) 


_als progressiv-periodisch mit einer raiumlichen Pseudoperiode / und einem Zuwachs 
des elektrischen Potentials pro Periode K bezeichnen. (Das gesamte System ist 
‘rotationssymmetrisch um die z-Achse.) Wir wollen zeigen, dass man einen Aus- 
‘schnitt der Lange n-/ (m = 1, 2, ...) eines solchen Feldes auf einem in z-Richtung 
‘endlichen, aus diskreten Widerstanden aufgebauten Netzwerk darstellen kann. Die 
an anderen Orten behandelte Beriicksichtigung der endlichen Maschenweite wollen 
wir hier voraussetzen. 

Trennt man den zwischen z, und z, + ” / liegenden Teil eines unendlichen Netzes 
heraus und verbindet entsprechende Punkte (das heisst solche mit gleichem 7) 
der neu entstandenen Grenzen paarweise iiber galvanisch unabhiangige Spannungs- 
quellen des inneren Widerstandes Null und der Spannung  K, so bleibt ein vor- 
-handenes, progressiv-periodisches Feld ungestort. Da man leicht zeigen kann, 
dass es bis auf eine additive Konstante nur eine mégliche Potentialverteilung gibt, 
“geniigt der Beweis, dass die urspriinglich vorhandene Potentialverteilung auch nach 
‘dem beschriebenen Heraustrennen einen mdglichen Gleichgewichtszustand dar- 
“stellt. Letzteres geht daraus hervor, dass erstens das Zuschalten der Quellen allein 
an der schon erfiillten Bedingung (1) offensichtlich nichts andert. Zweitens wird 
durch die gegenseitige Unabhangigkeit der Quellen auch nach dem Heraustrennen 
“die Gleichheit der durch entsprechende Punktepaare der Trennebenen z = 2) und 
z= 2,+ 0! hindurchfliessenden Stréme gewdahrleistet. Dies gilt allerdings nur, 
Ss — 


22) Abteilung fiir industrielle Forschung (AFIF), ETH. 
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Figur 1 
Achspotential eines Linearbeschleunigers (4 MeV) . 


wenn die zurtickgebliebenen Teile des unendlichen Netzwerkes dergestalt beschaffen 
sind, dass sie, langs der Trennebenen wieder zusammengesetzt, ein normales Netz- 
werk ergeben. Dies kann zum Beispiel durch Berandung des herausgegriffenen 
Stiickes mit doppelten Widerstandswerten erreicht werden. Die beschriebene An- 
ordnung ist in Figur 2 fiir m = 1 schematisch dargestellt, wobei die doppelt ausge- 
zogenen Linien die doppelwertigen Randwiderstinde andeuten sollen. Die abgebil- 
dete Elektrode wird elektrisch «freischwebend» belassen, kann aber zur Festlegung 
des Nullpunktes des nur bis auf eine additive Konstante bestimmten Potentials 
gebraucht werden. Zur Kontrolle durch die Bedingung (1) kénnen auch mehrere 
raumliche Pseudoperioden abgebildet werden. Figur 3 entspricht einer zur Priifung 
des Feldes eines elektrostatischen Beschleunigers durchgefiihrten Messung. Es 
sollte bei gegebener Elektrodenform untersucht werden, ob (und eventuell bis zu 
welchem Wert von 7) das achsennahe Feld innerhalb einer gewissen Toleranz 
homogen sei. 

Wahrend beim Abstecken von bis zu fiinf in iiblicher Weise auf die Potentiale V, 
bis Y, + 4K gebrachten Elektroden die Bedingung (1) wegen unsauberer Rand- 
bedingun gen nicht genau genug eingehalten werden konnte, erlaubte die hier vor- 
geschlagene Anordnung die Bestimmung des Potentials auf drei Stellen genau. 
D abei wurden die Spannungen K behelfsmassig von 30 durch Potentiometer abge- 
glichene 4,5-V-Trockenbatterien geliefert. Ein Teil der Messungsergebnisse ist in 
der Tabelle 1 zusammengefasst, welche deutlich zeigt, dass sich hier trotz der 
asymmetrischen Elektrodenform bis zu einer Achsentfernu ng von rund einem 
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Figur 2 
Darstellung eines progressiv-periodischen Feldes. 


Prittel des Elektrodenradius ein bis auf die dritte Dezimalstelle homogenes Feld 
ausbildet. 


Tabelle 1 


Es ist K = 1000 und V (Elektrode) = 0 
oe Se ea 
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Figur 3 
Potentialverteilung im Beschleunigungsrohr des 4-MeV-Beschleunigers. 
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Kiinstliche Erdsatelliten. 2. Sonderband der Fortschritte der Physik. Uber- 
setzung aus dem Russischen (Akademie-Verlag, Berlin 1959). 357 S., 146 Abb., 
gt Tab.; DM 38.—. 

Zu Beginn des Internationalen Geophysikalischen Jahres und zwar kurz vor 
dem Abschuss des ersten sowjetischen Erdsatelliten, namlich im September 1957 
erschien eine zweiteilige Sondernummer zum Band 63 der «Fortschritte der Physik» 
(Uspechi fiziéeskich nauk), die sich mit einigen der wichtigsten mit den kiinstlichen 
Erdsatelliten im Zusammenhang stehenden Problemen befasst und gleichzeitig das 
ausfiihrliche Programm darlegt, das sich die russischen Forscher auf diesem Gebiet 
fiir das IGJ vorgenommen hatten. Diese Sondernummer ist nun auch in deutscher 
Sprache zugdnglich (nachdem friiher schon eine amerikanische Ubertragung er- 
schienen war), vermehrt um einige zu diesem Thema gehdrende Arbeiten aus 
spadtern Faszikeln. — Vier Arbeiten befassen sich mit dem giinstigsten Abschuss, der 
Lebensdauer und den Bahnstérungen von kiinstlichen Erdsatelliten, eine weitere 
sehr einlasslich mit dem Mondflug. Die iibrigen Arbeiten sind im wesentlichen den 
durchzufiihrenden Experimenten gewidmet. Es werden die folgenden Problemstel- 
lungen behandelt: Untersuchung der Zusammensetzung der primaren kosmischen 
Strahlung und ihre Variation; die ultraviolette und korpuskulare Strahlung der 
Sonne; Bestimmung der Zusammensetzung, des Druckes und der Dichte der hohen 
Atmosphare (Exosphare); Untersuchung der Ionosphare; Messung der Haufigkeit 

| der Mikrometeoriten; geomagnetisches Feld und elektrostatisches Feld. Im grossen 
und ganzen sind dies gerade die Experimente, die vom Mai 1958 an mit Hilfe des 
Sputnik III (19586) realisiert wurden. Daneben finden sich noch Untersuchungen 
iiber die Siliziumsonnenbatterien, die optische Beobachtung kiinstlicher Erdsatelliten 
| und die Méglichkeit der Priifung der allgemeinen Relativitatstheorie. - Wenn diese 
-sorgfaltige und gut ausgestattete Ubersetzung leider auch etwas spat erscheint, so 
ist sie doch noch aufschlussreich genug (zu bemangeln ist allerdings die falsche 
-Riicktranskription vieler nichtrussischer Namen und auch von Zeichen: aus 
‘THIELE wurde TILLE, aus W. HOHMANN ein V. HocHMANN, aus F. R. MouLTON 
ein F. P. Mut’Ton usw.). E. RotH-DESMEULES 


Probability and Statistics. Herausgegeben von U. GRENANDER zum 65. Ge- 
_burtstag von HaraLp CRAMER. John Wiley & Sons, New York 1959, 434 Sire u2-005 
Dem Ansehen und der Beliebtheit des schwedischen Mathematikers HARALD 
CRAMER ist es zu verdanken, dass ein derart reichhaltiger Band mit Arbeiten der 
bekanntesten Forscher Westeuropas und Amerikas entstehen konnte. 
Die Ubersicht iiber den Inhalt dieses Buches zeigt die Tendenzen der Entwick- 
ung der Forschung in den letzten Jahren an. Ein grosser Teil der Arbeiten befasst 
sich mit stochastischen Prozessen, Construction du processus de W. Feller et H. P. 
McKean en partant du mouvement Brownien (Paut Lévy), Cramér’s theorem on mono- 
tone matrixvalued functions and the Wold decomposition (P. Masan1), Some non linear 
problems in probability theory (ULF GRENANDER), Some remarks on stable processes with 
independant increments (M. Kac), Statistical Analysis of stochastic processes with 
' stationary individuals (M. ROsENBLATT), Oder mit ihren Anwendungen, An introduc- 
tion to the measurement of spectra (JOHN W. TuKEy), Non linear prediction (P. MASANI 
and NorBERT WIENER). Diskrete Prozesse behandeln Unitary dilatation of Markov 
transition operators and the corresponding integral representations for transition- 
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q 
probability matrices (Davip G. KENDALL) und A Markov Chain Theorem (J. L. Doos). 
In The impact of stochastic process theory on statistics gibt M. S. BaRTLEIT eine 
Ubersicht iiber all die statistischen Probleme, die eigentlich nur mit Hilfe von 
stochastischen Prozessen umschrieben und korrekt gelést werden kénnen. } 

On combinatorial methods in fluctuation theory (W. FELLER) bringt einige neue 
Beweise im Zusammenhang mit dem arc sin Gesetz. Mit den Schatzungsmethoden > 
und Priifverfahren der modernen Statistik setzen sich JERZY NEYMAN in Optimal 
asymptotic tests of composite statistical hypotheses, HERBERT ROBBINS in Sequential 
estimation of the mean of a normal population, T. W. ANDERSON in Some scaling 
models and estimation procedures in the latent class model, EVELYN Fix, J. L. HopGEs 
and E. L. LEqmMann in The restricted chi-square test auseinander. Sehr willkommen 
ist auch die neue Ubersicht iiber den Stand der Forschung in Nonparametric statisti- 
cal inference von S. S. WiLKs und iiber die Forschungen zur Erreichung optimaler 
Versuchsanlagen, Design of linear experiments, von G. ELFYING. 

Schliesslich fiihren uns A survey of results in the collective theory of risks (CARL- 
Otto SEGERDAHL) und Ends and means in econometric model building (HERMAN 
O. A. WOLD) in weitere angewandte Gebiete, in denen H. CRAMER ebenfalls tatig 
wat. P. ScHMID 


La science et la théorie de l’information. ;Von LEon BriLLoutn (Masson 
& Cie, Paris 1959). 302 S., 69 Fig.; 48 NF. 

Die franzésische Ubersetzung des 1956 vom gleichen Autor erschienenen 
Buches «Science and Information Theory», welches in der ZAMP 8 (1957), H. 5, 
S. 428 besprochen wurde, liegt vor. Der Text weicht nicht ab, wurde aber da und 
dort erganzt, insbesondere wurde ein neuer Abschnitt beigefiigt iiber das Problem 
der divergenten Integrale in der Quantenelektronik. Die Ubersetzung dieses wert- 
vollen Buches ist sehr zu begriissen. Trotz des etwas zu klein gewahlten Druckes 
kann es bestens empfohlen werden. H. WEBER 


Orbit Theory. Proceedings of Symposia in Applied Mathematics of the Am. 
Math. Soc., Vol. 9. Herausgegeben von G. BrrKHOFF und R. E. LANGER. (Am. Math. 
Soc., Providence, Rhode Island, 1959). 195 S., 61 Bigs; $1220: 

Die Theorie der Bahnen eines Teilchens, die vor allem am Ende des letzten Jahr- 
hunderts in der Himmelsmechanik einen gewissen Hohepunkt erreicht hatte, ist heu- 
te wieder in den verschiedensten Gebieten von besonderer Bedeutung geworden. 
Aus diesem Grunde befasste sich das durch die American Mathematical Society im 
April 1957 veranstaltete neunte Symposium iiber angewandte Mathematik aus- 
schliesslich mit der Bahntheorie, um die Aufmerksamkeit der Mathematiker ver- 
mehrt auf dieses interessante und wichtige Gebiet zu lenken. Die zahlreichen neuen 
und alten noch ungelésten Probleme (periodische und unperiodische Bahnen, Stabili- 
tat, Extremalprobleme, Bahnstérungen) vermégen ihm Anlass zu neuen Entwick- 
lungen zu geben, und der Einsatz moderner Rechenautomaten lisst zudem in ganz 
anderer Weise numerische Untersuchungen durchfiihren als dies friiher méglich 
war. Die in den Proceedings wiedergegebenen 10 aufschlussreichen und bedeutungs- 
vollen Darstellungen bekannter Forscher umfassen einen grossen Teil des Gebietes 
der Bahntheorie: Teilchenbeschleuniger, geladene Teilchen in kosmischen M. agnet- 
feldern (kosmische Strahlung), Abschussbahn eines kiinstlichen Evdsatelliten (bezeich- 
nenderweise beansprucht die diesbeziigliche Arbeit von Srry mehr als ein Drittel 
des Bandes), Bahn eines Satelliten um einen abgeplatteten Zentralkérper, eingeschrénk- 
tes Dreikorperproblem, allgemeine Theorie der Planetenbahnen, Bahnen in relativisti- 
schen Gravitationsfeldern. Diese Proceedings vermitteln einen recht guten Uberblick 
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liber manche der aktuellen Probleme und: Methoden und sind fiir Mathematiker, 
Physiker, Astronomen und Ballistiker in gleicher Weise von Interesse. 
: E. RotH-DESMEULES 


Combinatorial Analysis. Vol. 10 ([Proceedings of Symposia in Applied 
Mathematics] American Mathematical Society, Providence 1960). 309 S.; $ 7.70. 
Dieser Band enthalt die Mitteilungen, die 1958 im Symposium tiber angewand- 
te Mathematik an der Columbia University gemacht wurden. Die Verfasser sind 
ausschliesslich in den USA tatige Mathematiker. Entsprechend der damaligen 
Organisation des Symposiums enthalt das Buch die folgenden Kapitel: 
1. Existence and construction of combinatorial designs ; 
2. Combinatorial analysis of discrete extremal problems; 
3. Problems of communications, transportation and logistics; 
4. Numerical analysis of discrete problems. 

Die Titel beweisen, dass die behandelten Fragen den verschiedensten Gebieten 
der Mathematik angeho6ren. In der Regel gehen die Autoren von einem praktischen 
Problem aus, um bei dessen Lésung in vielen Fallen Methoden der abstrakten 
Mathematik zu benutzen. Dieses Vorgehen zeigt, dass man die reine und ange- 
wandte Mathematik nicht voneinander trennen kann. 

Einzelne Beitrage sind besonders interessant, wie z. B. der Zusammenhang 
zwischen der Theorie der linearen Ungleichungen und dem Pauli-Prinzip (Verfas- 
ser: H. W. Kuun) oder die Abhandlung von R. BELLMANN iiber Combinatorial 
Processes and Dynamic Programming. W. SAXER 


Tafeln héherer Funktionen. Von JAHNKE-EmpDE-Léscu. 6. Auflage, bear- 
beitet von F. Léscu (B. G. Teubner, Stuttgart 1960). 318 S., 189 Fig.; DM 58.80. 
Die 6. Auflage dieses weltbekannten Werkes macht einen sehr vertrauten Ein- 
‘druck, indem die meisten der schénen und niitzlichen Reliefzeichnungen der hoheren 
Funktionen in drucktechnisch verbesserter Form wieder aufgenommen sind. Bei 
/naherer Durchsicht zeigt es sich aber, dass gegentiber der 5. Auflage ganz ent- 
scheidende Verbesserungen und Fortschritte erzielt wurden. Zundchst sind die be- 
handelten Funktionen iibersichtlicher in Klassen zusammengefasst; speziell wurde 
ein Abschnitt iiber Orthogonalpolynome aufgenommen. Vor allem aber sind die 
Tabellen neu iiberarbeitet worden und stehen jetzt leicht auffindbar hinter den 
Formelsammlungen und nicht mehr verstreut im Text. Gegeniiber der 5. Auflage 
sind erste Differenzen oder (wenn lineare Interpolation nicht gentigt) zweite Dif- 
ferenzen in allen Tafeln aufgenommen worden, so dass deren Gebrauch bequemer 
“wird. Einige neue Tabellen wurden berechnet: Insbesondere findet man nun bei 
den Bessel-Funktionen das Zumutbare; auch die Funktion J,,;, steht nun da, die 
- friiher fehlte. Dem Physiker werden die neuen Tafeln der Planck-Einsteinschen und 
- Debyeschen Funktionen willkommen sein. Die Bezeichnungen sind durchwegs mit 
denjenigen der modernen Standardwerke iiber spezielle Funktionen in Einklang 
gebracht worden; ein ausfiihrliches Literaturverzeichnis entspricht ebenfalls dem 
heutigen Stand der Forschung. Der «Jahnke-Emde» war schon immer unentbehr- 
lich, der « Jahnke-Emde-Lésch» wird noch unentbehrlicher sein. E. STIEFEL 


Handbuch der Physik - Encyclopedia of Physics. Bd, 48: Geophysik II. 
_Herausgegeben von S. FLUGGE (Springer-Verlag, Berlin 1957). 1045 S., 497 Fig.; 
DM 198.-. 
Der zweite Band der Geophysik, der von J. BARTELS redigiert wurde, enthalt 
nicht weniger als 17 Beitrage. Das einleitende und grésste Kapitel von A. ELIASSEN 


ain 


488 Buchbesprechungen — Book Reviews — — bibliographiques — ZAM) 


. ty 


und E. KLeinscumipt jr. ist der Ped logie, vor allem in theoreti- | 
scher Hinsicht, gewidmet. Die folgenden fiinf Beitrage umfassen das weite Gebiet | 
der Strahlung. F. M6rLER arbeitet in sehr klarer Weise das Gesamtproblem de | 
Strahlung in der unteren Atmosphare heraus; die Messmethoden werden jedoch 
nur gestreift. W. E. K. Mippreron befasst sich mit den optischen Eigenschaften_ 
der Atmosphare, wahrend Z. SEKERA unser Wissen iiber die Polarisation des_ 
Himmelslichtes zusammenstellt. Die Probleme, die mit der Diffusion der Strahlung 
an in der Atmosphare suspendierten Wassertrépfchen zusammenhangen, werden 
von J. Bricarp behandelt. Es ware wiinschenswert gewesen, wenn man hier auch 
die entsprechenden Forschungen auf dem Gebiet der Radarmeteorologie herange- 
zogen hatte. Sehr instruktiv ist die Abhandlung Ozon der Atmosphare von H. K. 
PAETZOLD und E. REGENER +, gibt sie doch einen sehr guten Uberblick iiber den | 
heutigen Stand des Ozonproblems. Nach kurzen Darlegungen iiber die geophysi-_ 
kalischen Aspekte von Meteoren von A. C. B. Lovet und iiber die Schallausbrei- 
tung in der Luft von E. F. Cox, folgt eine Einfiihrung in die Wolkenphysik von 
F. H. Lupiam und B. J. Mason. Es ist sehr erfreulich, dass dieser jiingste Zweig 
der Meteorologie in dem vorliegenden geophysikalischen Werk Aufnahme gefunden 
hat. Dann gibt R. MUHLEISEN einen ausgezeichneten Uberblick iiber die atmo- 
spharische Elektrizitat. Die folgenden Beitrage von H. U. SvERpDRuP, H. U. Rott, 
J. BarTtets, R. TomAscHek, A. DEFANT und W. KErRtTz haben zum Gegenstand die 
Ozeanographie und die Gezeitenkrafte. Der mehr als tausend Seiten umfassende 
Band schliesst mit einem Beitrag iiber physikalische Vulkanologie von S. SAKUMA 
und T. NaGata. 4 
In Anbetracht der immer umfangreicher werdenden und weit zerstreuten Spe- 
zialliteratur fullt diese mit grosser Sorgfalt vorgenommene Zusammenfassung 
unseres geophysikalischen Wissens eine schon seit langem bestehende Liicke aus. 
J.C. THams 


Differentialgleichungen fiir Ingenieure (2. Auflage). Von L. CoLratz 
(Verlag B. G. Teubner, Leipzig 1960). 197 S., 115 Abb.; DM 21.60. 

Der Teubner-Verlag stellt dieses Buch an die Spitze einer geplanten Folge von 
Lettfaden der angewandten Mathematik und Mechanik, die in erster Linie fiir die 


-Studierenden unserer Hochschulen gedacht sind. CoLLatz gibt eine meisterhafte 


Einfiihrung in das Gebiet der Differentialgleichungen, in der sowohl die prazise 
mathematische Theorie als auch numerische Methoden und Anwendungen gegen- 
einander ausgewogen sind. Nach einer Schilderung der elementaren expliziten 
Integrationsmethoden findet man eine in ihrer Einfachheit und Strenge bemerkens- 
werte Darstellung der Existenzsatze und der Poincaré-Bendixsonschen Klassifikation 
der Singularitaten. Der zweite Teil des Buches ist einer breitangelegten Theorie der 
linearen Differentialgleichungen gewidmet, die reich durch Beispiele aus der 
Schwingungslehre illustriert ist. Es folgen die Rand- und Eigenwertaufgaben bei 
gewohnlichen Differentialgleichungen, iiberleitend zu einer knappen, aber metho- 
disch interessanten Theorie der Orthogonalfunktionen und speziell der Kugel- und 
Zylinderfunktionen. In einer Ergdnzung behandelt das Buch erste Grundlagen fiir 
die Lésung partieller Differentialgleichungen und die fiir den Ingenieur wichtigen 
numerischen Methoden fiir gewéhnliche Differentialgleichungen. Wie in allen 
Biichern von Corrartz fallt die iibersichtliche Zusammenfassung der Ergebnisse in 
Tabellenform auf. 

Das Buch ist vorziiglich fiir die Verwendung als Begleittext von Vorlesungen 
fiir Physiker und Ingenieure geeignet. Der Studierende der Mathematik wird es als 
Grundlage fiir seine spateren vertieften Studien nehmen. E. STIEFEL 


“Lehrbuch : | 
‘ates righ baetss 3 


‘Direktor des Schweiz. Vereins 
fiir Schweisstechnik 


yon Dr.C.G.KEEL isa ; Schweisstechnik : 


Band I: 


Privatdozent an der Eidgenéssischen 


‘Technischen Hochschule, Ziirich Lichtbogenschweissverfahren ie yee 


J 


Band 22 Sammlung «Lehr- und H. andbiicher der I ngenteuy-Wissenschaften » 


(1961), 392 Seiten mit 448 Figuren und 113 Tabellen. Ganzleinen Fr./DM 64.— 


Inhalt: 
I. Allgemeines. Schweisstechnische Verfahren und Definitionen — Die Schweisstechnik — 


Definitionen — Kurze Kennzeichnung der schweisstechnischen Verbindungsverfahren — 


-Grundsatze fiir die Gruppierung der Verfahren — Die Warmequellen — der Grad der Erwar- 
mung — Gruppierung nach dem Schutz des Schweissprozesses — Der Begriff «Verbindungs- 


der Verbindungsmethoden — Normen: und Literatur. 
“II. Die elektrische Lichtbogenschweissung. Ubersicht iiber die Verfahren: Schematische 


-Gruppierung — Gruppierungsméglichkeiten der Verfahren — Ubersicht iiber die einzelnen — 


Verfahren — Der Lichtbogen — Allgemeines — Temperatur und Temperaturverteilung — Der 
_Werkstoffiibergang im Metallichtbogen — Der Spannungsabfall im Lichtbogen — Die Licht- 
_ bogencharakteristiken — Die Charakteristiken bei Wechselstrom — Die Stabilitat des Licht- 
_bogens — Die magnetische Blaswirkung des Gleichstromlichtbogens — Die Stromdichte bei 
der Lichtbogenschweissung — Die Lichtbogenschweissmaschinen — Allgemeines — Schweiss- 
-umformer — Schweissgleichrichter — Schweisstransformatoren — Arbeitstechnik, Einrich- 
-tungen, Hilfsmittel der Lichtbogenschweissung — Der Aufbau der Schweissnahte — Spezial- 
‘falle der Lichtbogenschweissung — Schweissfehler und deren Behebung — Die Schweisser- 


7: 


“schaften — Das Schweissgut — Elektrodenfabrikation, Lagerung und Handelsformen — Ab- 


 schmelzgréssen und Wirtschaftlichkeit — Die Wahl von SchweiSstromstarke und Elektrode — 


— Literaturverzeichnis — Die Schweissautomaten — Allgemeines — Halbautomaten — 
-Schweissautomaten mit sichtbarem Lichtbogen und Schlackenabdeckung des Bades — 
 Schweissautomaten mit verdecktem Lichtbogen — Mehrzweckautomaten — Spezialschweiss- 
_verfahren — Die Regelung der Schweissautomaten — Schweissmaschinen fiir Automaten — 
- Anwendungsgebiete und Zusammenfassung — Literaturverzeichnis. 


nt 
“Dem Verfasser ist es gelungen, den nur in Fachzeitschriften zerstreut publizierten Stoff, 


den er seit 1946 in seinen Vorlesungen behandelt, in einem leichtverstandlich dargestell- 
ten und iibersichtlich geordneten mehrbandigen Lehrbuch zusammenzufassen. 
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+ hoe ACZEL 
- Vorlesungen tiber Funktionalgleichungen 
‘ “und ihre Anwendungen 


Von Prof. Dr. J: AczéL, Professor an der Universitat Debrecen (Ungarn). 
331 Seiten mit 14 Figuren. Preis gebunden Fr. 37.50 (DM 37.50) (1961) 
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Inhalt: Inhaltsverzeichnis. Einleitung. — Erstey Teil: Gleichungen fiir Funktionen 4 


von einer Veranderlichen: 1. Gleichungen, in denen Veranderliche nicht nur unter 
Funktionenzeichen vorkommen, 2. Gleichungen mit lauter inneren Veranderlichen. 
3. Gleichungen mit mehreren unbekannten Funktionen. 4. Zuriickfiihrung auf Diffe- 
rential- und Integralgleichungen. — Zweiter Teil: Gleichungen fiir Funktionen meh- 
rerer Verdnderlicher: 5. Einfache Gleichungen. 6. Zusammengesetzte Gleichungen. 
7. Gleichungen mit mehreren unbekannten Funktionen mehrerer Veranderlicher. 
Zuriickfihrung auf partielle Differentialgleichungen. 8. Vektoren- und Matrizen- 


gleichungen.—Schlussbemerkungen. Einige ungeléste Probleme. Literaturverzeichnis. — 


Das vorliegende Buch ist das erste, das ein grosses Gebiet der Theorie der Funktional- 
gleichungen zusammenfassend bearbeitet. Die grossen Anwendungsgebiete, ins- 
besondere Wabhrscheinlichkeitstheorie, (nichteuklidische) Geometrie und Mechanik, 
werden auch eingehend behandelt. Das Literaturverzeichnis, das mehr als 1000 Ar- 
beiten aufzahlt, will gewissermassen auch als Bibliographie dieses Gebietes dienen. 
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